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SECTION  I 


INTRODUCTION 

The  general  problem  to  which  this  work  was  directed  was 
that  of  generating  a  usable  amount  of  energy  with  piezoelectric 
or  ferroelectric  materials.  The  input  to  the  active  material 
of  interest  was  to  be  a  forcing  function  on  one  face  of  the 
crystal,  and  the  desired  output  was  an  electrical  current  which 
would  be  dissipated  in  a  fixed  resistor. 

The  forcing  functions  which  would  be  applied  to  the  mate¬ 
rial  were  investigated  in  the  work  reported  on  in  Volume  I  of 
the  final  reports  of  the  contract  supporting  this  work.  The 
problem  of  particular  interest  was  to  maximize  the  energy  avail¬ 
able  for  dissipation  in  the  fixed  resistor. 

Important  factors  which  influence  the  amount  of  energy 
available  for  use  are  certain  material  properties  of  the  piezo¬ 
electric  or  ferroelectric  material,  the  geometry  of  the  material, 
the  pulse  signature  of  the  input  forcing  function,  the  fixity  or 
support  of  the  material,  and  the  circuits  used  with  the  crystal. 
The  work  done  for  this  study  included  consideration  of  each  of 
these  factors. 


SECTION  II 


LITERATURE  REVIEW 

Initial  studies  on  this  problem  included  a  search  of  the 
literature  to  determine  the  current  state  of  the  art  of  energy 
generation  from  piezoelectric  and  ferroelectric  crystals,  as 
defined  by  the  published  works  on  the  subject.  Because  of  the 
broad  nature  of  this  problem  and  to  present  the  information 
as  efficiently  as  possible,  the  list  of  papers  considered  as 
pertinent  are  given  as  the  references,  and  the  author's  ab¬ 
stract  of  each  of  the  more  significant  works  are  given  in  Ap¬ 
pendix  A.  The  first  four  references  give  information  of  a 
general,  introductory  nature,  and  do  not  have  abstracts  given 
for  this  reason. 
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SECTION  III 


REVIEW  OF  BASIC  RELATIONSHIPS  AND  THEORY 


Since  the  primary  purpose  of  this  study  is  the  determi¬ 
nation  of  the  energy  that  a  piezoelectric  or  ferroelectric 
device  can  deliver  to  some  electrical  load,  a  review  of  the 
necessary  computations  and  theory  are  in  order, 
a.  Charge  generated  by  the  crystal. 


Q, 


=  d  0  A 
em  m 


(1) 


The  above  relation  implies  a  linear  relation  between  charge 
and  force.  In  general,  this  is  not  true,  particularly  over 
a  wide  stress  range  for  quartz,  and  for  ferroelectric  mater¬ 
ials. 

b.  Capacitance. 


C  = 


K  Ae^ 
o  11 


(2) 


The  above  relation  implies  that  is  constant.  It  has  been 

shown  that  this  may  change  under  strong  shock  conditions. 

c.  Voltage  between  faces  of  cry 5 tel  for  open  circuit 
conditions  (RT  -*■  00 ) 


v  =  Qe/C 


(3) 


d. 

crystal 


Current  through  a  short  circuit  across  faces  of  a 
during  initial  passage  of  a  shock  wave. 


i  (t)  = 
sc 


d  a  (t)Ac 
em  m 


for  t<- 
c 


(4) 
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Reference  9  shows  that  this  relation  is  good  for  quartz  for 
the  time  for  passage  of  the  first  wave  through  the  crystal. 

The  subscript  "sc"  indicates  this  is  for  a  short  circuit 
from  one  face  of  the  crystal  to  another.  It  neglects  the 
effect  of  any  crystal  capacitance,  lead  capacitance,  or  line 
resistance,  which  cannot  be  zero. 

Since  Equation  (4)  is  the  basis  of  the  work  done  report¬ 
ed  herein,  the  limitations  of  this  equation  as  related  to 
piezoelectric  and  ferroelectric  materials  are  discussed  below. 

The  equation  is  based  upon  a  stress  induced  polarization 
which  is  a  linear  function  of  the  stress  in  the  direction 
normal  to  the  surface  upon  which  charge  is  being  collected. 
This  is  true  only  in  limited  cases.  One  of  these  is  a  state 
of  uniaxial  strain.  Another  is  uniaxial  stress,  where  the 
strain  in  other  directions  be  written  as  a  function  of  the 
stress  in  the  normal  direction  by  use  of  appropriate  consti¬ 
tutive  relations.  Similarly,  it  is  also  necessary  that  the 
electric  field  generated  in  the  crystal  be  a  linear  function 
of  the  normal  stress. 

It  is  also  necessary  that  stress  is  a  linear  function 
of  strain.  This  allows  the  relation  for  the  current  output 
to  be  written  as  a  linear  function  of  stress,  even  though  it 
is  the  strain  in  the  crystal  which  produces  the  current.  In 
order  to  keep  the  relation  one-dimensional,  it  is  assumed 
that  che  stress  on  the  charge  collecting  face  is  uniform  on 
that  face . 

Other  necessary  conditions  related  to  the  properties  of 
the  piezoelectric  material  are 

1.  Shock  velocity  is  independent  of  stress  amptitude. 

2.  Shock  velocity  is  steady. 

3.  Resistivity  of  the  material  is  infinite.  This  is 
particularly  true  in  short  circuit  mode. 

4.  Dielectric  constant  is  unchanged.  This  is  normally 
the  case,  except  for  high  strain  rates,  as  when  a 
shock  wave  passes. 


•i 


These  are  reasonably  valid  for  most  piezoelectric  or  ferro¬ 
electric  materials. 

Necessary  conditions  for  Equation  (4)  to  be  applicable 
which  are  valid  for  quartz  but  are  in  general  not  valid  for 
ferroelec tries  are  the  condition  of  linear  polarization  with 
stress  and  an  electromechanical  coupling  effect  which  is 
negligible.  For  ferroelectrics ,  the  polarization  is  approx¬ 
imately  linear  with  stress  for  very  low  stress  levels,  but 
for  high  stress  levels  it  is  highly  nonlinear.  The  results 
of  later  experimental  tests  make  this  quite  apparent. 

The  effect  of  a  large  electromechanical  coupling  coeffi¬ 
cient  is  to  convert  some  of  the  electrical  energy  generated 
in  the  crystal  back  into  mechanical  energy  through  induced 
stress.  For  a  material  such  as  quartz  where  the  coupling  co¬ 
efficient  is  small  (K  =  .1) ,  this  effect  is  negligible.  But 
in  a  ferroelectric  where  the  coupling  coefficient  can  be 
large  (K  =  .3  to  .6),  the  effects  are  certainly  not  negligi¬ 
ble.  They  are,  however,  difficult  to  calculate  and  account 
for. 

Even  though  this  equation  does  have  significant  limita¬ 
tion  on  when  it  is  valid,  it  will  be  used  as  a  basis  for  a 
first  approximation  in  the  studies  for  determining  energy 
outputs  from  dynamically  loaded  piezoelectric  and  ferro¬ 
electric  materials. 

e.  Current  through  a  short  circuit  across  faces  of  a 
crystal  during  passage  of  a  pulse  whose  duration  is  long 
relative  to  h/c. 


i 

sc 


(t) 


d  A  o  ( t)  Ac 
em  m 


(5) 


in  the  above  expression,  Ao(t)  is  the  stress  difference  be¬ 
tween  the  front  and  rear  electroded  surfaces  of  the  generator. 
This,  equation  is  applicable  for  any  degree  of  fixity  at  each 
face  of  the  crystal  and  is  based  upon  an  extension  of  Equa¬ 
tion  (4) ,  which  represents  the  output  of  the  crystal  as  ex¬ 
plained  on  the  basis  of  piezoelectric  theory  and  observed 
experiments . 
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f.  Energy  available  under  open  circuit  conditions 
(Rl  =  ->. 


E 


(6) 


In  this  expression,  E  is  in  joules  wnen  Q  is  in  coulombs  and 
C  is  in  farads.  The  energy  value  determined  from  this  expres¬ 
sion  is  not  a  realizeable  value  under  actual  conditions  of 
use,  i.e.,  when  the  generator  is  connected  to  a  fir Lte  load. 
Rather,  it  gives  an  upper  limit  on  the  amount  of  energy  avail¬ 
able. 

In  the  practical  application  of  a  piezoelectric  or  ferro¬ 
electric  as  an  energy  generator,  the  problem  is  to  find  the 
energy  delivered  to  the  external  load.  The  generator,  circuit, 
and  external  load  are  shown  in  Figure  1.  The  inductive  effect 
of  the  elements  has  been  assumed  negligible.  If  the  following 
further  assumptions  are  made 

^  ”  rl 


RDist  <<:  ^ 

(7) 

ZBond  °  ' 


» 


♦ 


cT  .  <•<  c  +r  . 

Load  p  Dist 

then  the  circuit  can  be  approximated  as  shown  in  Figure  2. 
If  it  is  assumed  that  the  energy  generated  by  the  crystal 
in  the  real  circuit  of  Figure  2  is  the  same  as  when  in  a 
short  circuit,  the  following  expression  is  obtained: 


i 
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Figure  1.  Generator,  Circuit  and  External  Load  Schematic 


Figure 


for  t  <  h/c 


(8) 


Sc(t)  "  iq<t)  " 


de^<t)Ac 


This  expression  is  subject  to  all  the  assumptions  previously 
discussed  in  connection  with  Equation  (4) . 

The  current  through  the  load  resistor  may  be  found  from 
circuit  theory 


i  + 
c 


i 


R 


C 


de 

Total  dt 


where  e  is  the  voltage  drop  across  C 
above  equations  may  be  written  as 


Total ' 


The  first  of  the 


de 
T  dt 


1R 


(9) 


This  equation  may  be  solved  for  the  current  delivered  to  the 
external  load  if  the  input,  i  (t)  ,  is  known. 

g.  Energy  delivered  to  the  external  load. 


(10) 


It  must  be  remembered  that  iR  in  the  above  equation  has  been 

obtained  utilizing  i  .  which  is  based  upon  the  assumption  that 

Si 

the  load  resistance  in  the  circuit  is  zero.  However,  the 
energy  available  is  now  a  function  of  RT .  Experiment  has  shown 
that  for  R^  small,  say  less  than  100  ohms,  the  equation  may 
still  be  used  for  quartz  (Ref.  9) . 

As  before,  in  many  practical  applications  the  rise  time 
of  the  stress  pulse  will  not  be  less  than  the  transit  time 
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of  the  pulse  through  the  piezoelectric  generator.  It  is 
necessary  to  extend  Equation  (8)  to  cover  these  cases.  To 
do  this.  Equation  (8)  is  reinterpreted  to  the  form: 


iq(t) 


d  Ao  (t)Ac 
em  m 

- E - 


(11) 


where  Ao  (t)  is  the  stress  difference  between  the  front  and 
m 

rear  surface  of  the  crystal . 

The  above  discussion  is  based  upon  the  energy  generation 
from  the  piezoelectric  effect.  Ferroelectric  materials  also 
generate  charge  in  a  mode  different  from  the  piezoelectric 
one,  i.e.,  by  depoling  from  passage  of  a  shock  wave.  Refer¬ 
ence  6  gives  a  discussion  of  phenomenon. 

The  maximum  charge  release  possible  under  these  condi¬ 
tions  is  determined  from  the  expression 


Q  =  P  A 
Tnax  r 


(12) 


Actual  tests  show  that  all  of  this  amount  cannot  be  obtained 
(Ref.  10) .  The  percentage  which  can  be  obtained  is  dependent 
on  such  factors  as  the  intensity  of  the  shock  wave  and  the 
planeness  of  the  shock  front. 

In  general,  a  greater  charge  may  be  released  from  the 
remanent  polarization  than  with  the  piezoelectric  effect. 

This  phenomenon  has  been  given  the  name  "charge  bulge." 
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SECTION  IV 


SUMMARY  OF  IMPORTANT  PROPERTIES  CF  COMMON  MATERIALS 


There  are  many  materials  which  exhibit  a  piezoelectric 
or  ferroelectric  effect,  and  consideration  must  be  given  to 
the  ones  most  suitable  for  use  in  energy  generation.  This 
section  gives  information  on  the  more  important  properties 
of  materials  currently  listed  in  catalogs  of  certain  manu¬ 
facturers.  The  balues  have  been  obtained  primarily  from  the 
manufacturer's  literature,  but  some  has  come  from  technical 
journals.  The  values  are  summarized  here  for  convenience. 
Information  is  given  on  the  following  materials: 


Quartz 

HD- 11 

HS-21 

G-53 

HDT-31 

G-1408 

G-1500 

HST-41 

G-1512 

G-2000 


Manufactured  by 
Gulton  Industries,  Inc. 


Manufactured  by 
Clevite  Corporation 


Blanks  on  these  sheets  indicate  that  the  information  was 
either  not  readily  available  or  the  quantity  was  not  appli¬ 
cable  for  that  material. 
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REMANENT  POLARIZATION: 

*  An  Introduction  to  Piezoelectric  Transducers,  Valpey  Corp. , 
May  1965. 

**  Piezoelectricity,  by  Hans  Jaffe,  Reprinted  from  the  Ency¬ 
clopedia  Britannica,  Copyright  1961. 

#  A  Primer  on  Ferroelectricity  and  Piezoelectric  Ceramics, 
Bernard  Jaffe,  Electronic  Research  Division,  Clevite  Corp., 
Engr.  Memo  #60-14,  1960. 
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MATERIAL:  HD- 11 

(Eariuro  Titanate) 

DENSITY: 

SPECIFIC  WEIGHT: 

ELASTIC  MODULI: 


ONE-DIMENSIONAL 
WAVE  VELOCITY: 

CURIE  TEMPERATURE: 

FREE  DIELECTRIC 
CONSTANTS : 


PIEZOELECTRIC  CHARGE 
COEFFICIENTS : 


PIEZOELECTRIC  VOLTAGE 
COEFFICIENTS: 


COUPLING  COEFFICIENTS: 


REMANENT  POLARIZATION: 


MANUFACTURER:  Gulton  Industries,  Inc. 

3 

=  >5.6  gm/cm 
=0.2  lb/in3 


,E 

33 

=  13.0 

X 

1010 

newtons 

2 

m 

or 

1.89 

X 

107 

psi 

,E 

11 

=  13.7 

X 

1010 

newtons 

2 

m 

or 

1.99 

X 

107 

psi 

E 

55 

=  5.7 

X 

1010 

newtons 

2 

or 

C .  83 

X 

107 

psi 

m 

c  =  4.8  x  105  cm/sec  or  1.89  x  105  in/sec 
°C  >135 


00 


pcb  pcb 

newton  or  3  82  lb 

sUfes  or-i33  !ir 
sSisir  -  sss 

-3 

x  10  volt  meters/newton 

-3 

x  10  volt  meters/newton 

_3 

x  10  volt  meters/newton 


ul/cm3 
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MATERIAL:  HS-21 

(Barium  Titanate) 
DENSITY: 

SPECIFIC  WEIGHT: 
ELASTIC  MODULI: 


MANUFACTURER:  Gulton  Industries.  Inc. 

=  >5.6  gm/cm^ 

=  0.2  lb/in^ 

vE  _  n  t  „  1n10  newtons  . 

=  11*1  x  10  - ^ -  or  1*61  x  10  psi 

m 

VE  _  in  c  „  ...10  newtons  .  ...7 

Y..  =  11.5  x  10  - = -  or  1.67  x  10  psi 


ONE-DIMENSIONAL 
WAVE  VELOCITY: 


vE  a  .  .  A10  newtons  ..  .. .  ,.7 

Yj.,.  =  4.4  x  10  - ^ -  or  0.64  x  10  psi 

in 

5  5 

c  =  4.45  x  10  cm/sec  or  1.75  x  10  in/sec 


CURIE  TEMPERATURE: 


*C  >125 


FREE  DIELECTRIC 
CONSTANTS : 


=  1150 


=  1050 


PIEZOELECTRIC  CHARGE 
COEFFICIENTS : 


PIEZOELECTRIC  VOLTAGE 
COEFFICIENTS : 


148  -E£b..-  or  658  E£b 

newton  lb 

-50  -ESh—  or  -220  Egb 

newton  1J3 

225  -E£$_  or  1000  E£b 

newton  lb 


=16  x  10  volt  meters/newton 


=  -5.2 x  10  volt  meters/newton 


=25  x  10  volt  meters/newton 


COUPLING  COEFFICIENTS : 


-  0.51 


=  0.18 


=  0.48 


REMANENT  POLARIZATION: 


=  50  y  ccul/cii,‘ 
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MATERIAL:  G-53  MANUFACTURER:  Gulton  Industries,  Inc. 

(Lead  Zirconate  Titanate) 

DENSITY:  =  >7.6  gm/citi 

SPECIFIC  WEIGHT:  =  0.27  lb/in3 


ELASTIC  MODULI:  Y^  =  6.5 

X 

1010 

newtons 

5 

or 

0.94 

X 

107 

psi 

m 

’'ll  •  8-! 

X 

1010 

newtons 

2 

or 

1.18 

X 

107 

psi 

in 

VE  _  3.8 
*55  ~ 

X 

1010 

newtons 

2 

or 

0.55 

X 

107 

psi 

m 


ONE- DIMENSIONAL 

WAVE  VELOCITY:  c  =  2.92  x  1C'*  cm/sec  or  1.15  x  10 5  in/sec 

CURIE  TEMPERATURE:  °C  >330 


FREE  DIELECTRIC 


CONSTANTS : 

K3 

=  720 
\ 

K1 

=  960 

PIEZOELECTRIC  CHARGE 
COEFFICIENTS : 

d33 

=  190 

pcb 

newton 

or 

845 

pcb 

nr~ 

d31 

=  -84 

pcb 

newton 

or 

-370 

pcb 

lE~ 

d15 

=  300 

pcb 

newton 

or 

1330 

lb 

PIEZOELECTRIC  VOLTAGE 

COEFFICIENTS:  g^  =  30  x  10  volt  meters/newton 

g^  =  “13  x  10  3  volt  meters/newton 


COUPLING  COEFFICIENTS: 


_  3 

36  x  10  volt  meters/newton 
0.60 
0.29 
0.64 


REMANENT  POLARIZATION:  Pr  =26  Pcoul/cm2 


MATERIAL:  HDT-31  MANUFACTURER:  Gulton  Industries,  Inc. 

(Lead  Zirc^nate  Titanate)  o 


DENSITY: 


SPECIFIC  WEIGHT: 


ELASTIC  MODULI: 


=  >7.6  gm/cm 
=0.27  lb/in3 

,,E  ^  .  A10  newtons  .  -.7 

Y33  =  6.7  x  10  - 2 -  or  0*97  x  10  psi 

in 


„E  .  ,  . -10  newtons  -  ,  _  . 

Y11  =  S.l  x  10  - 2 -  or  1*18  x  10  psi 

rn 

„E  _  . -10  newtons  „  ,  „7 

Y^-  =  3.2  x  10  - = -  or  0.46  x  10  psr 


ONE- DIMENSIONAL 
WAVE  VELOCITY: 


5  5 

c  =  2.98  x  10  cm/sec  or  1.17  x  10  in/sec 


CURIE  TEMPERATURE: 


>330 


FREE  DIELECTRIC 
CONSTANTS : 


K3  =  1300 


K 1  =  1350 


PIEZOELECTRIC  CHARGE 
COEFFICIENTS: 


L,  =  280  —  or  1243 

33  newton  lb 


d^  12  0 


newton 


or  -534 


d15  =  360 


newton 


or  1600 


PIEZOELECTRIC  VOLTAGE 
COEFFICIENTS : 


^33  =  23  x  10  volt  meters/newton 


=  -11  x  10  volt  meters/newton 


g. _  =  29  x  10  volt  meters/newton 
15 


COUPLING  COEFFICIENTS:  k33  =  0.66 

k31  =  0.35 

k, r  =  0.59 
Id 


REMANENT  POLARIZATION:  Pr  =27  ycoul/cm 
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MATERIAL:  G-1408 


MANUFACTURER:  Gulton  Industries,  Inc. 


DENSITY: 

= 

>7.5 

gm/cm3 

SPECIFIC  WEIGHT: 

= 

0.27 

lb/in3 

ELASTIC  MODULI: 

ye 

*33 

= 

8.2 

iglO  newtons 
m^ 

or 

1.19 

x 

10 

ye 

*11 

= 

9.0 

.-,10  newtons 

X  ±u  - = - 

or 

1.30 

X 

1C 

in 

ye 

*55 

= 

4.0 

newtons 

2 

or 

0.58 

X 

1C 

m 

ONE-DIMENSIONAL 

WAVE  VELOCITY: 

c 

— 

3.3  x 

105  cm/sec  or  !• 

3  x 

105 

CURIE  TEMPERATURE: 

FREE  DIELECTRIC 
CONSTANTS : 


PIEZOELECTRIC  CHARGE 
COEFFICIENTS: 


>300 

1000 

1250 


PIEZOELECTRIC  VOLTAGE 
COEFFICIENTS: 


33 


31 


15 


*33 


=  20C 

pcb 

newton 

or 

890 

pcb 

lb 

=  -80 

pcb 

newton 

or 

356 

pcb 

lb 

=  315 

pcb 

newton 

or 

1401 

pcb 

lb 

COUPLING  COEFFICIENTS: 


g3i 


*15 


‘33 


‘31 


‘15 


REMANENT  POLARIZATION:  P  =22  ycoul/cm 


=22  x  10  3  volt  meters/newton 
-3 

=  -9. Ox  10  volt  meters/newton 

-3 

=29  x  10  volt  meters/newton 
=  0.60 
=  0.26 
=  0.60 

2 


7 

psi 


7 

psi 


7 

psi 


in/sec 
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MATERIAL:  G-1500  MANUFACTURER:  Gulton  Industries,  Inc. 

(Lead  Zirccnate  Titanate) 


DENSITY: 

SPECIFIC  WEIGHT: 
ELASTIC  MODULI: 


ONE- DIMENSIONAL 
WAVE  VELOCITY: 

CURIE  TEMPERATURE: 


=■  >  7 . 6  gm/cm~ 


=  0.27 

lb/in3 

.E 

33 

=  4.9 

newtons 

or 

0.71 

X 

io7 

psi 

• 

m3 

,E 

11 

=  6.3 

^glO  newtons 
m3 

or 

0.91 

X 

107 

psi 

* 

.E 

55 

=  2.2 

,rt10  newtons 
x  10  2 

or 

0.32 

X 

iO7 

psi 

c 

°C 


m 


=  2.54  x  103  cm/sec  or  1  x  103 
>360 


in/sec 


FREE  DIELECTRIC 
CONSTANTS : 


K3  =  1700 
=  1700 


PIEZOELECTRIC  CHARGE 
COEFFICIENTS : 


PIEZOELECTRIC  VOLTAGE 
COEFFICIENTS: 


d33  -  370 


d  3  =  -166 


d15  =  540 


_E£b  1642  Egb 

newton  lb 

Pcb  or  -739  Pcb 
newton  lb 

or  2400  Egfe 
newton  lb 


.-3 


g33  =25  x  10  volt  meters/newton 


*31 


=  -11  x  10  volt  meters/newton 


COUPLING  COEFFICIENTS: 


-3 

g15  =  36  x  10  volt  meters/newton 
=  0.69 
k31  =  0.34 
k^^  =  0.65 


REMANENT  POLARIZATION: 


P  =  24  yccul/cm 
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MATERIAL:  HST-41  MANUFACTURER:  Gultcr.  Industries,  Inc. 

(Lead  Zirconate  Titanate)  . 


DENSITY: 

SPECIFIC  WEIGHT: 
ELASTIC  MODULI: 


=  >7.6  gm/cm 
=0.27  lb/in3 


ye 

*33 

=  5.9 

X 

1010 

newtons 

2 

or 

0.86 

X 

1C7 

psi 

in 

ye 

*11 

=  7.0 

X 

1010 

newtons 

2 

or 

1.01 

X 

1C7 

psi 

in 

ye 

*55 

=  1.8 

X 

1010 

newtons 

2 

m 

or 

0.26 

X 

107 

psi 

ONE-DIMENSIONAL 
WAVE  VELOCITY: 


=2.79  x  10  cm/sec  or  1.10  x  10"  in/sec 


CURIE  TEMPERATURE: 


>270 


FREE  DIELECTRIC 
CONSTANTS : 


K3  =  1800 
K.  =  2100 


PIEZOELECTRIC  CHARGE 
COEFFICIENTS : 


PIEZOELECTRIC  VOLTAGE 
COEFFICIENTS : 


=  325 


h  =  -157 
31 


newton 


newton 


=  625  P£*L 


newton 


1445 

lb 

-698  Pcb 
lb 

2780  E£b 


g33  =22  x  10  volt  meters/newton 


g31  =  -11  x  10  volt  meters/newton 

_3 

g^3  =37  x  10  volt  meters/newton 


COUPLING  COEFFICIENTS:  k33  =0.66 

*31  =  0.35 
*15  *  °-69 


REMANENT  POLARIZATION:  Pr  =  23  tcoul/cm^ 
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MATERIAL:  G-1512  MANUFACTURER:  Gultcn  Industries,  Inc. 

(Lead  Zirconate  Titanate) 


DENSITY: 

=  >7.4 

SPECIFIC  WEIGHT: 

=  0.26 

ELASTIC  MODULI: 

vf3  =  5.4 

Y11  *  6'3 

is  • 2-6 

gm/cm 

lb/in3 


X 

1010 

newtons 

2 

or 

0.78 

X 

107 

psi 

in 

X 

1010 

newtons 

2 

or 

0.91 

X 

107 

psi 

ni 

X 

O 

f— 1 

o 

f— 1 

newtons 

2 

or 

0.38 

X 

1C7 

psi 

m 


ONE- DIMENSIONAL 
WAVE  VELOCITY: 


5  c 

=  2.69  x  10  cm/sec  or  1.06  x  10  xn/sec 


CURIE  TEMPERATURE:  °C  >240 


FREE  DIELECTRIC 

CONSTANTS:  K3  =2600 

Kx  =  2200 


PIEZOELECTRIC  CHARGE 
COEFFICIENTS: 


PIEZOELECTRIC  VOLTAGE 
COEFFICIENTS : 


COUPLING  COEFFICIENTS: 


d 

d 

d 

g 

g 

g 

k 

k 

k 


33 

= 

500 

pcb 

newton 

or 

2220 

pcb 

lE~ 

31 

= 

-232 

pcb 

newton 

or 

-1030 

pcb 

lb 

680 

pcb 

or 

3100 

pcb 

15 

newton 

lb 

33 

= 

20  x 

10~3 

volt 

meters/newton 

31 

= 

-9. 3x 

10~3 

volt 

me  te  r s/newton 

15 

= 

35  x 

10-3 

volt 

meters/newton 

33 

= 

0.72 

31 

= 

0.37 

15 

= 

0.78 

REMANENT  POLARIZATION:  Pr 


=  23 


ycoul/cm 


2 
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MATERIAL:  G-2000 
(Lead  Ketaniobate) 
DENSITY: 


SPECIFIC  WEIGHT: 
ELASTIC  MODULI: 


ONE-DIMENS IONAL 
WAVE  VELOCITY: 

CURIE  TEMPERATURE : 

FREE  DIELECTRIC 
CONSTANTS : 


PIEZOELECTRIC  CHARGE 
COEFFICIENTS : 


PIEZOELECTRIC  VOLTAGE 
COEFFICIENTS : 


COUPLING  COEFFICIENTS: 


REMANENT  POLARIZATION: 


MANUFACTURER :  Suitor  Industries,  Inc. 
=  >5.8  gm/cm3 
=0.21  lb/in3 


ye  = 
33 

4.7 

X 

1010 

newtons 

2 

or 

0.68 

X 

10 

in 

ye  = 
11 

4.0 

X 

1010 

newtons 

2 

or 

0.58 

X 

10 

m 

ye  = 
55 

1.8 

X 

1010 

newtons 

2 

or 

0.26 

X 

10 

m 


psi 

psi 

psi 


c  -  2.85  x  10 5  cm/sec  or  1.12  x  l05in/sec 
°C  >400 


K3  =  250 
Kx  =  250 


d33  =  80 

pcb 

or 

356 

pcb 

I3BT 

newton 

d31  =  -10 

pcb 

newton 

or 

44.5 

pcb 

lb~ 

d15  -  115 

pcb 

newton 

or 

511 

pcb 

lb 

g33 

36  x 

io-3 

volt 

meters /new ton 

g31  = 

-4 . 5x 

10-3 

volt 

meters/newton 

g15  = 

50  x 

io-3 

volt 

meters/newton 

k33  = 

0.38 

k31  " 

0.04 

k15  = 

0.33 

P 

r 

2 

-  Mcoul/cm 

MATERIAL:  Ceramic  "B' 

DENSITY: 

SPECIFIC  WEIGHT: 
ELASTIC  MODULI: 


MA>UFACTURER:  Clevite  Corp. 

3 

=  5.55  gm/cm 
=0.2  lb/in3 

,,E  tit  .  ~  10  newtons  .  , . 

Y^^  =  11*1  x  10  - ^ -  or  1.61  x  10  psi 

in 

Y^  =  11.6  x  10  - 2 -  or  1*68  x  10  psi 

nr 


5  . 


ye 

55 

s 

ONE-DIMENSIONAL 

WAVE  VELOCITY: 

c 

= 

4.48  x 

:  105 

cm/sec  or 

1.65  x  10 

CURIE  TEMPERATURE: 

°C 

115° 

FREE  DIELECTRIC 
CONSTANTS : 

K3 

= 

1200 

K1 

= 

1300 

PIEZOELECTRIC  CHARGE 
COEFFICIENTS: 

d33 

= 

149 

pcb 

■  —  or 

newton 

662 

pcb 

lb 

d31 

= 

-58 

-EES—  or 
newton 

258 

pcb 

lb 

d15 

= 

242 

-E£6_  or 
newton 

1070 

pcb 

lb 

PIEZOELECTRIC  VOLTAGE 
COEFFICIENTS : 

g33 

= 

14. lx 

10-3 

volt 

meters/newton 

g31 

= 

-5.5x 

10'3 

volt 

meters /newton 

g15 

= 

21  x 

10-3 

volt 

meters/newton 

COUPLING  COEFFICIENTS: 

k33 

= 

0.48 

k31 

= 

-0.194 

k15 

= 

0.48 

REMANENT  POLARIZATION: 

P 

r 

= 

-  coul/cm3 
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MATERIAL:  PZT-2 


DENSITY: 

SPECIFIC  WEIGHT: 
ELASTIC  MODULI: 


ONE-DIMENSIONAL 
WAVE  VELOCITY: 

CURIE  TEMPERATURE: 


MANUFACTURER:  Clevite  Corp. 

2 

=  7.6  gm/cm 
=  0*27  lb/in^ 


c  = 

°C  370 


FREE  DIELECTRIC 
CONSTANTS : 


PIEZOELECTRIC  CHARGE 
COEFFICIENTS: 


PIEZOELECTRIC  VOLTAGE 
COEFFICIENTS : 


COUPLING  COEFFICIENTS: 


REMANENT  POLARIZATION: 


-3 

g33  =  38.1*  10  volt  meters/newton 
-3 

g31  =-15. lx  10  volt  meters/newton 

_3 

=  50. 3x  10  volt  meters/newton 

^33  =  ®  •  8^6 
k31  =  -0.28 

k15  =  0.701 

2 

Pr  = - coul/cm 


MATERIAL:  PZT-4 


MANUFACTURER:  Clevite  Corp. 


DENSITY: 

-  7.5 

gm/cnt 

SPECIFIC  WEIGHT: 

=  0.27 

lb/in^ 

ELASTIC  MODULI: 

ye 

33 

=  6.6 

x  1010  newfns 
in 

or  0.96 

x  107 

psi. 

*11 

=  8.2 

„  1Q10  newtons 

or  1.19 

x  107 

psi 

ONE- DIMENSIONAL 
WAVE  VELOCITY : 

CURIE  TEMPERATURE: 


c  =  2.97  x  10"  cm/sec  or  1.17  x  lO^in/sec 
°C  328 


FREE  DIELECTRIC 
CONSTANTS : 


1300 

1475 


PIEZOELECTRIC  CHARGE 
COEFFICIENTS: 


PIEZOELECTRIC  VOLTAGE 
COEFFICIENTS: 


289  —  or  1240  ^ 

newton  lb 

-123  —  or  -546 

newton  lb 

496  —  or  2201 

newton  lb 


=  26.  lx  10  volt  meters/newton 


--H.ix  10  volt  meters/newton 


COUPLING  COEFFICIENTS: 


g-^  =  39. 4 x  10  volt  meters/newton 
k33  =  0.70 
k31  =  0.334 
k, c  =  0.71 


REMANENT  POLARIZATION: 


=  -  coul/cm 
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MATERIAL :  PZT5A 
DENSITY: 

SPECIFIC  WEIGHT: 
ELASTIC  MODULI: 


MANUFACTURER:  Clevite  Corp. 
=  7.  75  gm/cm3 
=0.28  lb/in3 


ONE-DIMENSIONAL 
WAVE  VELOCITY: 

CURIE  TEMPERATURE: 

FREE  dielectric 
CONSTANTS : 


PIEZOELECTRIC  CHARGE 
COEFFICIENTS : 


PIEZOELECTRIC  VOLTAGE 
COEFFICIENTS : 


COUPLING  COEFFICIENTS : 


c  = 

°C  365 


d33  -  374 

pcb 

newton 

or 

1660 

pcb 

IE~ 

d31  *  -171 

pcb 

newton 

or 

-760 

pcb 

lb 

d15  =  594 

pcb 

newton 

or 

2598 

pcb 

Ib~ 

933  24. 8x  10  3  volt  meters/newton 

931  ~  11,4x  10  volt  meters/newton 

g15  =  38*2x  10  3  volt  meters/newton 

k33  =  °-705 

k31  =  ~0*344 
kl5  =  0.685 


REMANENT  POLARIZATION:  Pr  =  _  coul/cm2 
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MATERIAL:  PZT-5H 


MANUFACTURER:  Clevite  Corp 


DENSITY. 

SPECIFIC  WEIGHT: 
ELASTIC  MODULI : 


ONE- DIMENSIONAL 

WAVE  VELOCITY:  c 

CURIE  TEMPERATURE:  °C 

FREE  DIELECTRIC 
CONSTANTS:  K3 

K1 

PIEZOELECTRIC  CHARGE 
COEFFICIENTS:  d33 


PIEZOELECTRIC  VOLTAGE 
COEFFICIENTS:  g33 

g31 

*15 

COUPLING  COEFFICIENTS:  k33 

1- 

'31 

k15 

REMANENT  POLARIZATION:  P 


=  7.5  gm/cm3 
=  0.27  lb/in3 


193 


=  593 

pcb  or 

2639  pcb 

newton 

lb 

-  -274 

pcb  or 

-1220  pcb 

newton 

Ib^ 

=  741 

pcb 

— -  or 

3299  pcb 

newton 

lb 

=  19. 7x 

10-3  volt 

meters /newton 

=-9.11x 

10-3  volt 

meters /newton 

=  26. 8x 

10  3  volt 

meters/newton 

=  0.752 
=  -0.388 
=  0.675 

= - coul/cm3 


SECTION  V 


COMPUTATIONS  UTILIZING  BASIC  RELATIONSHIPS  AND  CIRCUITS 

EXAMPLES  OF  COMPUTATIONS  USING  BASIC  RELATIONSHIPS 

In  the  previous  two  sections  of  this  report,  the  basic  re¬ 
lationships  between  pertinent  quantities  were  given,  as  well 
as  values  for  material  properties  of  several  materials .  In 
this  section,  sample  calculations  are  given  to  illustrate  the 
use  of  the  equations  and  material  properties,  and  to  give  in¬ 
dications  of  actual  numerical  values  for  the  various  parameters 
of  interes,  such  as  currents  and  voltages.  These  sample  cal¬ 
culations  will  be  based  upon  the  following  two  crystals: 

1.  Quartz  disk,  1/2  in.  dia.  x  0.10  in.  thick 

(identified  by  "g") 

2.  HST-41  disk,  1/2  in.  dia.  x  0.10  in.  thick 

(identified  by  "h") 

a.  Crystal  capacitance 


KoAen 


(8. 85xl0-12)  J(i)  (0.115) 
<3=  c - 07TCT - 


=2.0  pfd 


h:  C  = 


( 8 . 85xl0-12)  (45.6) 

OTTO 


C  =  795  pfd 


b.  Voltage  across  capacitor  for  open  circuit  condi¬ 
tions  : 
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V  =  Q/C 


q:  V  =  (10.2  pcb/lb)/2.0  pfd 

V  =  5.1  volts/lb  at  applied  force  (for  one- 

dimensional  stress) 

h:  V  =  (1440  pcb/lb) /795  pfd 

V  =  1.82  volts/lb  of  applied  force  (for  one¬ 

dimensional  stress) 

c.  Short  circuit  current  across  faces  of  a  loaded 
crystal: 


d  o  ( t)  Ac 

isc  (t)  =  - ^ -  f  for  t  <  h/c. 


,_2, 


(10.2  pcb/lb) (.196  in  ) (216,000  in/sec) 
q:  Sc  =  -  - OTTO - 


4.3y  amp/psi,  for  t  <  .  46y  sec  (for  one¬ 
dimensional  stress) 


h:  i 


sc 


(1440  pcb/lb) (.196  in2 ) (2 . 77xl05  in/sec) 

0.10 


=  79 2y  amp/psi,  for  t  <  .  36y  sec  (for  one¬ 
dimensional  stress) 

d.  Energy  stored  in  capacitor: 


E  =  o! 

h  2C 
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(10.2  pcb/lb) 
2  x  2.0  pfd 


-12 

E  =  26.2  x  10  joules/lb  of  applied  force 


t..  -  _  (1440  pcb/lb)2 

h-  E  "  unffpid- 


-9 

E  =  1.31  x  10  joules/lb  of  applied  force 

Equation  (11)  will  be  used  to  illustrate  the  current  and 
energy  output  which  would  be  expected  from  a  piezoelectric 
generator  in  an  application  where  the  input  pulse  rise  time 
is  long  compared  to  the  transit  time  through  the  pulse.  The 
case  considered  will  be  one  for  which  the  input  pulse  is  a 
ramp  type,  and  the  impedance  of  all  elements  of  the  system 
are  matched.  Figure  3a  shows  the  pulse  and  system  components, 
and  3b  shows  the  expected  current  time  pulse.  Figure  4  shows 
the  effect  of  various  rise  times  of  input  pulse.  The  charge 
generated  (area  under  i(t)  curve)  is  the  same  in  each  case, 
but  since  energy  involves  current  squared,  the  energy  released 
is  increased  as  the  rise  time  of  the  pulse  is  decreased. 


Stress  wave  Propagating 
With  Velocity  c  and 


Note:  The  mechanical  impedance  of  elements  F,G,  and  R 
are  assumed  equal  for  this  example. 

(a) 


Current 

(i) 


Time  (t)  in  Units  of  h/c 
(b) 


Figure  3.  Er.ample  for  Finite  Rise  Time  Pulse 


Current 

(i) 


Time  (t) 

(b) 


Figure  4.  Effect  of  Rise  Time  on  Crystal  Output 
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STUDY  OF  BEHAVIOR  OF  IDEAL  CIRCUIT 

Equation  (9)  gives  the  differential  equation  relating  the 
current  generated,  i  or  igc#  to  the  current  passing  through 
the  load  resistor  in  which  the  energy  will  be  dissipated.  For 
any  experiment  conducted  on  a  crystal,  i  may  be  measured,  but 

X\ 

i  is  the  one  predicted  by  theory.  Equation  (8) .  This  section 
gives  the  results  of  a  study  made  of  the  effect  of  the  circuit 
parameter  R^  and  CT  on  iR  for  assumed  forms  of  i^. 

Equation  (9)  wqs  programmed  for  evaluation  on  a  digital  com¬ 
puter.  This  solutions  was  evaluated  for  two  forms  of  i^(t): 
a  ramp  function  and  a  haversine  pulse.  These  two  functions 
were  selected  since  they  approximate  the  excitation  applied  to 
the  crystal  for  two  of  the  types  of  tests  planned. 

The  parameters  describing  these  pulses  were  defined  as 
follows : 

t  =  Rise  time  for  ramp  function  and  pulse  duration 

IV 

for  haversine. 

RC  =  Time  constant  of  the  circuit. 

XK  =  t  /RC 

K 

Runs  were  made  for  values  of  XK  from  10.0  to  10,000  for  the 
ramp  function,  and  100  to  50,000  for  the  haversine.  It  was 
believed  that  these  ranges  would  more  than  amply  cover  that 
encountered  in  the  experimental  work. 

Figure  5  shows  examples  of  the  results  for  each  of  the 
inputs.  The  curves  plotted  represent  those  for  the  parameter 
values  resulting  in  the  greatest  deviation  in  the  two  currents. 

Table  I  gives  the  computer  output  obtained  for  the  runs 
made  for  the  ramp  input.  Similar  results  are  given  in  Table 
II  for  the  haversine.  The  one  undefined  term  in  these  tables 
is  "DFCRP . "  This  is  the  percent  difference  between  i  and  iR, 
based  upon  the  maximum  i^.  These  tables  may  be  used  to  esti¬ 
mate  the  difference  in  i^  and  iR  for  experimental  data  involv¬ 
ing  similar  inputs.  A  study  of  these  tables  shows  that,  in 
general,  the  difference  between  the  two  quantities  is  small 
and  that  significant  difference  occur  only  from  t  /RC  less 

x\ 

than  500 .  This  condition  could  be  easily  satisfied  for  the 
test  conditions  utilized  in  th_s  work. 


32 


7\|  'Yl 


TABLE  t.  Results  for  Ramp  Function  Input  to  Ideal  Circuit 


VA1UE  OE  RC  *  4.0000E-08 


VAU'E  OF  XK  =  in. 00  ’0 

TIME  lMr«E'-*E\iT  CURRENTIQ 
2„ 000o000r-f,9  R.OOOOOOnE-O? 
4.000o000F-nfl  1  .  OOOOOOOE-O 1 

6-OOOoOOJE-C3  1.5000000E-01 
P.OCOOOOOF-oR  2.0000000E-01 
1.  OOOoOOOF-n7  ?.5000000E-0l 
1.2000000F-O7  3. OOOOOOOE-O 1 

1-4000000F.-0?  3.5000000E-01 

1.6000000F-07  4.0000000E-01 

L.g00n00OF-o7  4.5000000E-01 

,OOOoOOOF-fl7  ^.0000000E-01 
,2000000E-07  S.5000000E-01 

2,40000n0E-(l7  A,  OOOOOOOE-O  l 
2.6000000F-07  4.SOO0000E-01 

2,800O000F-07  7.0000000E-01 

3.0 00O00OE-O7  7.5000000E-01 

3.2000000F-07  rt.OOOOOOOE-01 
3.400000Uc-o7  8.5000000E-0 i 

3. 6000000E-n  7  9. OOOOOOOE-O 1 

3.800OO00F-07  9.5000000E-01 

4  .  OOOqOOoE  —  07  1 ,O0OOOOOE*0O 


RISE  TIME  s  4.0000E-0? 

CUPREMTTR  OFcRP 

1 .0653068E-02  3,934Aq34E*00 

3.A7R7944E-02  A. 321 2056E*00 
7,231301 AE-02  7,7A86q84E*00 

1.1353353E-ni  8.64664t2F*00 

1 .5820850E-01  9.1791500F+00 

2.0497871E-O1  9,5021293E*00 

?.S3ol974E-f'l  9.A980?62F*00 

3.0183156E-01  9.8168436E*00 

3,5111 090E-O 1  9.8889100F+00 

4.00A7379E-01  9 ,9326?05E*00 

4.5n4086RE-Ol  9,9S913H3F*00 

5.0024788E-01  9,9752l25E*00 

S.5ni5034E-Ol  9.9E a9a56E*00 

6.0009119E-01  9.9908Bl2E*Oo 

A.5005531E-01  9.9944A92E*00 

7,0O03355E-01  9,996A454E*Oo 

7.5002035E-O1  9.9979653E*0n 

8, 0  00 1 234E-0 1  9,9987A59E*00 

R.5O00749E-O1  9,9992S15E*00 

9, 0000454F-0 1  9.9995460F400 


VALUE  OF 
VAt.  JF  qE 
TIM-  iMCRE'-E’-lT 
Z.OOOOOOOE-OS 
4„OOOnOOOK-n8 
6.000000 0^-09 

H,  0000000^-09 

I .  0000000F-O7 
I . 200000OF-O7 
i . 4000000F-07 
I .6000000E-07 
I. flOOOOnOF-07 
2.O00o000f;-07 
2.2000000F-o7 
2.4OOO0OOF-O7 
2.600f)00  )E-«7 
2.8o0o000r-o7 
3, OOOOOOOE-07 
3. 20 On  0 00F-07 
3.4000000E-07 
3.6000000F-07 
3. 8O0OO00F-O7 
4. OOOoOOOE-07 


PC  =  8.0000E-09 

XK  =  50,00  in 

C'JRREMTI  0 
S. OOOOOOOE-O? 

1 #000000nE-0l 
1 .5000000E-01 
2, OOOOOOOE-O 1 
?,50000nor.O] 

3  •  OOOOOOOE-O l 

3.5000000E-0] 

4.0000000E-01 

4.5000000E-01 

5.0000000E-01 

5.5000000E-01 

A.OOOOOOOE-Ol 

A. 5000000E-01 
7, 0  0  00  0  OOE-0 1 
7 .5000000E-0 1 

B. OOOOOOOE-Ol 
3.5000000E-01 
9.0000000E-01 
9 • 500  00  0  OE-O 1 
1 .OOOOOOOE+OO 


RISE  TIME  *  4.0000E-07 


CUPREMTIR 
3,1641 700E-02 
R.0134759E-O2 
1 .3001106E-01 
1 ,8000091 E-01 
2.3000007E-01 
2.8000001E-01 
3.3000000E-01 
3.R0O0000E-O1 
4.3000000E-01 
4,8000000E-01 
5, 30000 00 E-01 
5,8000000E-ni 
6.3000000E-01 
6.8000000E-01 
7.3000000E-01 
7,8000000E-01 
«.3000000E-01 
B.8000000E-01 
9.3000000E-01 
9.8000000E-01 


OEcRp 

1,B35R300E*Oo 
1,9865?41F*00 
1 ,998R93HE*Oo 
1,9999o92E*0o 
1 ,9999925E*00 
1 ,9999994F*Oo 
1 ,9999999E*Oo 
2, 0000Q00F  40o 
2.0000000E+00 

2,oooooooe*o'o 

2,OOOOOOOE*On 

2,0000000E*0n 

2,0000000F*0q 

2,OOOOOOOF>Oo 

2,oooooooe*oo 

2,ooooooof46o 

2.0000000F *00 
2, OOOOOOOE*O0 
2,OOOOOOOF*00 
2,OOOOOOOE*On 
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TABLE  I 


(continued) 


VALUE  OF" 
VALv'E  OF 

time  Increment 

2-0000000F-O8 
4.0000000F-OR 
6.  OOOOOOOE-08 

r.ooooooof-ob 

i  .000n00<iF-n7 
t  ,2oOnOOOF-n7 
1  .4000000C-07 
l  .600000OE-A7 

1  tSOOOOO«iF-o7 

2  .OOOnOOOf -o7 
2,2oonnooE-o7 
2,400000«F-07 
2,600o0noc-07 
2,800nO0UF-o7 
3.0000000F-07 
3,200o000F-o7 
3 .40000OOE-07 
3.6000000F-O7 
3.8000000F-07 
4.0000000E-07 

VALUE  OF 
VALUF  OF 

T  *  me  Increment 

2.0000000F-P8 

4.00000O0E-O8 

6.0000000E-O8 

8.oooooooe-oR 

I  .OOOOOOOF-07 
I  ,200O00UF>n7 
l  .4000000^-07 

I «6000000E-n7 
I  .8000000E-07 
2.000O000E-07 
2 ,200o000E-07 
2,400oO00E-p7 
2.&00O000E-07 
2.80000n0E-o7 
3.0000000E-07 
3.2000000F-07 
3.4000000E-07 
3.600000OE-07 
3.8000000F-07 
4.0000000F-07 


RC  *  4.0000E-09 

XK  a  100,00  0 
CURRENT  10 

R. OOOOOOOE-O? 

1  .OOOOOOOE-Ol 
1 .500000OE-O1 
?”«  OOOOOOOE-Ol 
2.5000000E-0 1 
3. OOOOOOOE-Ol 
3.8000000E-0 1 
4. OOOOOOOE-Ol 
4.5000000E-01 
ri.OC00003E-01 
5.5000000E-01 
f>. OOOOOOOE-Ol 
4.SOOOOOOE-01 
7. OOOOOOOE-Ol 
7.50000C9E-0 1 
M.OOOOOnOE-Ol 
8.5000000E-0 1 
•9. OOOOOOOE-Ol 
9.500O000E-01 

i  .nooooooE*oo 

RC  a  8.0000E-10 
XK  =  500,00('0 

CURRF  4TT(J 

S. OOOOOOOE-O? 

1  .OOOOOOOE-Ol 
1.5000000E-01 
?. OOOOOOOE-Ol 
2.8000000E-01 
3.0000000E-0 1 
3.SOOOOOOE-Q1 
4.00000OOE-01 
4.5000000E-01 
5.0000000E-01 
5.5000000E-01 
ft.OOOOOOOE-Ol 
^,5000000E-01 
7. OOOOOOOE-Ol 
7.5000000E-01 
A, OOOOOOOE-Ol 
8,50O0QO0E-0, 

9. OOOOOOOE-Ol 
9.5000000E-01 
1 .OOOOOOOE^OO 


RISE  TIME 
CURRENT  I R 
4.0Q67379E-O? 

R. 0000454E-A? 

1  .4000000F-01 
1  .9000000E-A1 
2.4O0O000E-O1 
2.9000000E-01 
3.4000000E-01 
3.9000000E-01 
4.4000000E-01 
4.9000000E-01 

S. 4000000E-01 
8.9OO00O0E-01 
f».4000000E-01 
4.9000000F-ni 
7.4000000E-01 
7.9000000E-01 
M.4000000E-A1 
8.9OO0000E-01 
9.4000000F-01 
9.9000000E-O1 


RISE  TIME 
CURRENTIR 
4.8O0000QE-02 
9.8000000E-02 
1.4800000E-01 
1.9800000E-01 
2.4ROOOQOE-0 1 
2.9800000E-01 
3.4ROOOOOF-0) 
3.98OO000E-ni 
4 . 4WO0000E-A 1 
4.9800000F-01 
5.4800000E-O1 
5.9AOOOOOE-01 
A.4ROOOOOE-A1 
4.9ROOOOOE-nl 
7 , 4ROOOOOE-n 1 
7.9R00000E-ni 
A.4»oOOOOE-ni 
M.9ROO0O0E-01 
9.4R0000OE-01 
9.9300000E-01 


a  4, rtOOOE-07 

OFcRR 

9.93?fi?05F-Ol 
9.9995460E-01 
9.99999&9E-01 
l,OOOOr>OOF*Oo 
1.0000000E*00 
1 ,0000000r+0n 
1.0000000F*00 
1.0000000E*00 
1.0000000F*00 
1.0000000F*00 
1 ,OOOOOOOF*Oo 
1 ,OOOOOOOF*On 
1.0OOO000F*0o 
1.0000000F*00 
1 . 00  000  00E*O0 
1 . OOOOOOOF*Oo 
1.0000000E*00 
1 ,OOOOOOOE*Oo 
1.0090000F*no 
1 ,OOOOOOOE*00 


4.O000E-O7 
nFcRP 
2, OOOOOOOE-Ol 
2.OOOO00OF-01 
2, OOOOOOOE-Ol 
2, OOOOOOOE-Ol 
2.0000000E-01 
2.000000PE-01 
2 . 00  OOOOOF-O i 
2.00OO000F-O1 
2.0000000^-01 
2.O000000E-01 
2.000AOOOE-0] 

2 , 00  OOOOOF-O 1 
2.0000QOOE-01 
2.0000000F-01 
2.0000000F-01 
2.0000000F-01 
2. OOOOoOOE-O 1 
2,OOOuOOOF-Ol 
2.0000000E-01 
2.0000000F-01 
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TABLE  I 


(continued) 


VALUE  OF 
v\L:>F  OF 
TIME  ImorE^EmT 
2, 00000 Oof -o« 
4.OOOOO0OF-08 
6.000000GF-08 
8.000noa<'^-oP 

I  .0000000^-07 
l  .2000000F-07 
l  .4000000*  -07 
I .6000000F-07 

1. eoooooof -o7 

2.  OOOOOOOE-0 7 
2.200O00OF-O7 
2.400OO0:'F-o7 
2.6000000F-07 
2.800OOO0F-O7 
3. OOOnOOOF -o7 
3.2o0onoor-o7 
3.400OOOOE-O7 
3.600O000E-O7 
3.8000000F.-o7 
4.0O0O00OF-O7 

VALUE  OF 

VAU'F  OF 
TIME  INCREMENT 
2.000000UF-08 
4.000OC0UF-08 
6.0000000F-08 

a.oooooooE-oft 

1.0000000F-O7 
1.200O00OF-07 
1 .4000000F-07 
1.6000000F-*17 
1.8000000E-07 
2.O000000F-O7 
2.2OOOOOOE-07 
2. 40000  OOF-07 
2.6000000E-07 
2.8000000F-07 
3.000O000F-O7 
3.2000000E-07 
3.4000000E-07 

3. LOOOOOOF-07 
3.9000000E-07 
4.0000000E-07 


RC  *  4.0000E-10 


XK  =  1000,00 

CURRENTIQ 
5.O000O00E-02 
T.nOOOOOOF-Ol 
1.5000090E-01 
2. OOOOOOOE-0  1 
2.5000000E-01 
3. OOOOOOOE-0 1 
T.5000000E-01 
4, OOOOOOOE-0 1 
4.5000000E-01 
5. OOOOOOOE-Ol 
5.5000000E-01 
ft. OO0OOO0E-O1 
ft.8000000E-01 
7.0000000E-01 
7.S000000E-01 
q.OOOOOOOE-Ol 

H. bOOOOOOE-Ol 
q.OOOOOOOE-Ol 
9.5000000E-01 
1  .OOOOOOOE  +  OO 

RC  a  8.0000E-1 

XK  a  SOOO.OOO 
C’JWRE'xITIQ 
5.0000000E-02 
\ .C000000E-01 

I . 5000000E-01 
2. OOOOOOOE-O  l 
2. 500 0000 E- 01 
3. OOOOOOOE-O 1 
3. 500000  OF- 01 

4.0000000E-01 
4.5000000F-01 
5 . O00000  OF-O 1 
5.5000000E-01 
ft  •  n  0  0  0  00  OF-O 1 
ft.5000000E-01 
7. 0000000  E->)1 
7 , 50000  QOE-0 ] 
3.0000000E-01 
5. 5000000 E- 01 
0 , 00  OOOQOE-U 1 
R.5000000E-01 
1  .0000000E  +  00 


10  RISE  TIME 

CURRENTIR 
4.900O000E-02  1 

9.9000000E-02  1 

1 . 4qoO000E-0 1  1 

1 .9900COOE-01  1 

2.4900000E-01  1 

2 • 9900000E-0 1  1 

3.4900000E-01  1 

3.9O00000E-01  1 

4.4900000E-01  1 

4.9900000E-O1  1 

S.49000O0E-01  1 

8.9900000E-01  1 

ft.4900000E-01  1 

5.9900000F-01  1 

7 . 49000  OOE-0 1  1 

7.9900000E-O1  1 

3.4900000E-01  1 

A.9900000E-01  1 

9.4900000E-01  1 

9.9900000E-0 1  1 


0  RISE  TIME  : 

CURRENTIR 
4.9800000E-02  2 

9.9POOOOOE-02  2 
1.4980000E-01  2 

1 .9980000E-01  2 

2.49R0000E-01  2 

2.998O000E-O1  2 

3.49BOOOOE-01  2 

3.99R0000E-01  2 

4.4980000E-01  2 

4 . 99R0000E-O 1  2 

S.4980000E-01  2 

5.9980000E-01  2 

ft. 4980 00 OE ”01  2 

ft . 9980000E-O 1  2 

7.498000OE-O1  2 

7 • 9990000E-0 1  2 

R.4980000E-01  2 

5.99R0000E-01  2 

9.4980000E-01  2 

9.9980000E-01  2 


a  4.0000E-07 
DFcRP 

• 0090000E-0 1 
.OOOOOOOE-Ol 
.OOOOOOOE-Ol 
,00 000 OOE-O 1 
,OOOOOOOE-Oi 
,OOOOOOOF-Ol 
.OOOOOOOF-Ol 
.OOOOOOOE-Ol 
,OOOOQOOE-Ol 
.OOOOOOOF-Ol 
,OOOOOOOE-01 
, OOOOOOOF-Ol 
, OOOCOOOE-O 1 
.OOOOOOOF-Ol 

.OOOOQOOE-Ol 

.OOOOOOOF-Ol 

.OOOOOOOF-Ol 

.OOOOOOOE-Ol 

.OOOOOOOF-Ol 

, ooooqoof-oi 


-  4«oOOOE-07 
OFcRP 

.OOOOOOOE-02 
.OOOOOOOE-02 
,0000000E-02 
,  OOOOOOOE-0? 
.OOOOOOOF-02 
.OOOOOOOE-02 
.0000000F-02 
.OOOOOOOE-O? 
,OOOOOOOE-02 
.OOOOOOOE-02 
, OOOOOOOF-02 
.OOOOOOOE-02 
.OOOOOOOE-02 
,0000000E-02 
,0000000E-02 
.OOOOOOOE-02 
.OOOOOOOF-O? 
.OOOOOOOE-02 
.OOOOOOOE-02 
,  OOOOOOOE-02 


36 


TABLE  I  (continued) 


VALUE  OF 

VALUE  oF 
TIME  ImorEmEV'T 
2.000n000c-0B 
4.0000000E-08 
5,000o00‘jF-o8 
8. OOOoOOOF-08 
1.0000000F-07 
1.2000000E-07 
1.4000009F-07 
1.6000000E-n7 
1.8000000F-07 
2,OOOnOOOF>o7 
2,2000«OOE-07 
2.400OP00F-O7 
2.6000000F-07 
2.B00n00»F-o7 
3,oooooooe-o7 
3.2000000E-07 
3.400O00QE-O7 
3.600!}000E-n7 
3.8000000E-07 
4.000000PF-07 


PC  =  4, OOOOE-1 1 

=  loooo.ofno 
cmprfntiq 

S.OOOOOOOF-O? 

1  .AO00000E-01 

i  .soooooof-oi 

2.O000000E-01 
?.SOOOOOOE-Ol 
l.OOOOOOOE-Ol 
S.'jOOOOODE-Ol 
OOOOOOOE-O  1 
4.S000000E-01 
S.OOOOOOOE-ni 
S.500000OF-01 
fy,  300000nE-Ol 
■’-..‘lOOOOOfiE-Ot 
7.ooooooof-oi 
7.C>OOOOOOE-'U 
a.ooooooof-oi 

^.SOOOOOOE-Ol 

o.OOOOOOOE-Ol 

7.5000000E-01 

l.OOOOOOOEOQ 


PISE  TIME 
CUPREMTIR 
4.990O000E-"2  1 

9.9900000E-f'2  1 

1  .4990000E-01  1 

!  .9Q900O0E-01  1 

2.4R90000E-01  1 

2.9990000E-0]  1 

3.4nq0000E-Ol  1 

3.9Q90000E-01  1 

4.4990000E-P1  1 

4.99R0000E-f'l  1 

5.4990000E-O1  1 

F.9990000E-"1  1 

8.4990000E-01  1 

4.9990000E-01  I 
7.4990000E-nl  1 

7,9990000F-M  1 

».4990000E-ol  1 

4.999000PE-ni  1 

9.4990000E-ni  1 

9.9990000E-0)  1 


=  4.0000E-07 

OFcRP 

.ooooooof-o? 

.OOOOOOOE-O? 

.ooooooce-o? 

.OOOOOOOE-O? 

. OOOOOOO F" 02 
.OOOOoOOE-02 
.OOOOOOOF-O? 
.OOOOOOOF-O? 
.OOOOOOOF-O? 
.OOOOOOOF-O? 
.OOOOOOOF-C? 
.OOOOOOOF-O? 
.OOOOOOOF-02 
.OOOOOOOF-O? 
.0000000F-02 
.OOOOOOOE-O? 
.OOOOOOOF-O? 
.O0OOO0OF-O2 
.OOOOOOOF-02 
.OOOOOOOF-O? 
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TABLE  II. 


Results  for  Haversine  Input  to  Ideal  Circuit 


VALUE  OF  Rf!  =  1.0000E-07 


value  of 
t I '■'F'  INrPE^ENiE 
S.OOOOOOOE-<)7 

1 , ncoonooE-Oft 
1  .BnpoOOOE-oft 
2.0000000E-0*> 
?.  SO00000E--16 
3.  rmoonooE-t^ 
3.  SOOOOOOE-''*' 
4.oop0000E-oft 
4. 5noooooE-iift 
5. OOOOOOOE-Oft 
5.c>O0OO00E-o5 
ft.OOOOOOOE-nft 
ft.ftooooooE- vs 
7.0700000E-06 
7.5"0OO()0E-OS 

B.onoOOOOE-(if» 

8.  5P00000E-0ft 

9.  O000000E-0ft 
g.BononooE-ofe 
1 . O00O000E-O5 


X  K  s  100. 000 '• 
C-IPPFMTTO 
2.449l2«3E-02 
9.558584! F-02 
2..'62ftOP4E-nl 

3.45732O4E-01 
5,oo3lM2E-nl 
b.54fl60?3E-ol 
7.9425r»56E-ol 
9,  ',4fln559E-0l 
9.  7‘,7n379E-ol 
9.^9999*0E-o 1 
9. 75312B9E-01 
9. n416?0^E-Ol 
7 • g3??7?8E-0 1 

6  •  L>3ftftft4  3E-0  1 

4.99o5lft3E-ol 
3.4452974E-01 
2.pS?3»SIE-o1 
9.4023B28E-02 
2.^1  oic>o3E-o2 
1  ,t,9«93?3E«06 


RISE  TIMC  = 


C  JRREMTIR 
1.S&71B69E-02 

7.375043ftE-0? 

1 ,02O83ft3E-Ol 
3 .  lft5bl5bE-0 ! 

4.5901000E-0] 
ft, 244943BE-01 
7, S778266E-0 1 
p, 34B3757E-01 

O. 54191R2E-01 

P. 9807145E-01 
9.4315745E-0] 
9.20910P7E-01 
P , 1 7429ft9E“0  J 
ft,R285H8E-0l 

5.303ft019E-Cl 

3.7409467E-O1 
2.3168511E-01 
1.1476105E-01 
1.5576971E-02 
1 a990l590E-03 


l.OOOOE-Os 

DFCPP 
7.8194138E-01 
1  .ftP15405E*0o 
2.41771&1E*00 
2.9170487L^OO 
3.1306123F*00 
3.0374850F*00 
2.ft4ft7902E*00 
1.99b802lE*00 
1.1S11971E*00 
1.9281494E-01 
-7.8445&40E-01 
-1,68407B*E*OO 
•2.4202410E*00 
-?.9185044E*00 
-3.13085ftlE*00 
-3.0364931E*00 
-2.644ft597E*00 
-1.99374lftE*00 
-1.1475067E+0O 
-1.8885ft00F-01 


VALUF  OF  R 

VALUE  of 
T I  •  F  increment 
5. POOOOOOE-O  7 
1 . OOOOOOOE -06 
1 .8OOOOOOE-O6 
?. OOOOOOOE -Ob 
a.ft^oooooE-ob 
3. nooooooE -oft 
3.  v^onnoot-oft 
4. OOOOOOOE-Oft 
4  .BOOOOOOE-'16 
5.  OOOOOOOE-Oft 
5.5000000t-oft 
ft.  OPOOOOOE-^ft 
ft.  bp 00000E-fi6 

7.PO000O0E-06 

7. ^O00000t-o6 

8.  nPOOOOOE-Oft 
R.=>O00000t-0*> 

9. nP00000E-0ft 
o.SOOOOOOE-'ib 
1.OO0O000E-05 


=  2.0000E-0P 

K  =  500.000 

C'  'RPF mT  I Q 
2.4491 7P3E-02 
9»  S5ft5p4l E -02 
2. pft26Q04E-O  1 
3.4873?04£-Ol 
S.  ■>031ftl2E-ol 
ft.F48ft9?3E-0 1 
7.-U  pcn5ftE-0 1 
9.  '|4p0559F-0l 
9.7570379E-"1 
9,99995ftPE-ol 
9.7531?59E-ol 
9.040ft?o8E-0 1 
7.-»  32?7?8E-p  1 
ft. 53ftftft43E-P 1 
4.^905]ft3E-0l 

3.445?a74F-nl 
2.P5238F1E-01 
9.4R?3ft?ftE-0? 
2.M  0 1  9p3E-0? 
1 ,59P99?3E-0ft 


RISE  Til* 
CJRRENiTIR 
2.P623580E-02 
9. 19343Q5E-0? 
2.012214IE-O1 
3.3977B31E-01 
4.9403I35E-O1 
6 . 488ft9 1 1 E-0 1 
I.3912289E-01 

9.0105270E-01 
9, 7369334E-01 
9,9992854E-<!l 
9.7710010E-O1 
9.0770001E-01 
7,9827l4ftE-f)3 
ft,596227lE-03 
5.0533649E-01 
3.50S2749E-01 
?,1036158E-01 
9.3570136E-0? 
?.ftl04823E-02 
ft, 473  2694E-05 


»  =  1.0000t-05 

oecrp 

1.8677035E-01 
3.6314451E-01 
5.0394317E-01 
5,9537297£-01 
6.2847697E-01 
6.0001215E-01 
5.1P76706E-01 
3.752B868E-01 
2.0l0^509E-01 
7.1060898E-03 
-I.8752906E-01 
-3.6379290E-01 
-5.0441771E-01 
-5.9562718E-01 
-ft , 2848596F-0 1 
-5.9977502E-01 
-5.1230705E-CI 
-3.7465085E-01 
-2.0029192E-01 
-6.3ll37ftlE-03 
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TABLE  II 


(continued) 


VALUE  OF 

VALUE  o F 
TI ,JF  Increment 
OrtOOOOOF -0  7 

l.monnooE-06 

1.5O00000E-O6 

2.0000000E-06 

2.5O00O00E-0& 

3»of>oonoot-c6 

3.Sn0OO00£-06 

4. nn00O00E-o6 
4.SO0O000E-06 

5. nooonooE-n6 
5. 6O00000F-06 
ft.OPOOOOOf -06 

6.  c>000000t-(i6 
7.0100000E-06 
7.5000000E-06 
R.OOOOOOOE-nfc 
8.S00O000E-06 
9.OO00000E-O6 

g.^ooooooe-oft 

) .noOOOOOE-OB 


RC  =  l.OOOOE-OR 

XK.  -  lOOO.OOOr, 
C'lPRpMTrQ 
2.4491 2b3E-p? 
9.5565B41E-0? 
2.0626 0R4E-0I 
3.4573?p4E-nl 
5.P031612E-01 
6.5486923E-01 
7.9425o56e-M 
9.04«r.S59E-nl 
9.757p?79E-nl 
9.R999Q60E-01 
9.75312P9£-o1 
9.f)406?o8E-cl 
7 .932P7?PE-n 1 
6.5366643E-01 
4.9905163E-O1 
3.44S2974E-rl 
2.''523R5lE-ol 
9.4P236,?8E-o2 
2.4 1  0190 3E -02 
1.59R9323E-06 


RISE  TIM- 
CJRPENTIR 
2.3538534E-C? 
9.3733924E-0? 
2.0372922E-01 
3.4274873E-01 
4.971733BE-01 
6, 51 87494E-0 1 
7.9169804E-01 
9.0294492E-O1 
f1.747l724E-01 
R.99R83P2E-0] 
9.7626942E-01 
9.05R9723E-O1 
7.9576126E-01 
6,  56651 00E-0 1 
5.0219440E-03 
3.47522P3E-01 
2.0778R70E-01 
9, 5681 0r/RE-02 
2.50P4680E-O2 
1 .3405295E-05 


1 .OOOOt-OP 

decrp 

°.b?74913E-02 

1.8319164E-01 

2.5316197E-01 

2.9833126F-01 

3.142744RE-01 

P.9942977E-01 

2.5525137E-01 

1.8606725E-01 

9.8655017E-02 

1.5780355E-03 

-9.5653539E-02 

■1.8351440E-01 

'2.5339725F-01 

2.9845601E-01 

3.1427648E-01 

2.99308P1E-01 

2.5501932E-01 

1.8574683E-01 

9,B2776c3F-0? 

1.1806363E-03 


VALUE  OF  RC  *  2. 0000E-09 


VALUE  OF  XK  = 
TImf  rMfprMFMT 


5.0000000E-<)7 

1  .OOOOOOOE-O6 

1 .8OOOOOOE-06 

?.  OO00O00E-06 

2.9000000^-06 

3.0000000L-06 

3.9000000E-U6 

4. 0000000 E -06 

4.5000000t -06 

5.O000000E-O6 

5.SO00000E-O6 

6.0000000E-06 

6.5000000E-06 

7.0000000E-O6 

7 ,9000000E-o6 

R.o000000E-o6 

8.P000000E-06 

9.0O00000E-O6 

9.5000000E-06 

1  .OOOOOOOE-05 


5000. dCOn 
C'lppFMJIQ 
2.449i2R3E-n2 
9»5565B4lF -02 
2.r.626np4F-(ii 
3.4573204E-01 
rJ.r’03l6i2E-nl 
6 . 54869?3E'-r  1 
7.94?5f>56E-ol 
9,04P08R9E-nl 
9.7S70379E-r.l 
9.99990^oE-o 1 
9.7531PR9E-01 
9.O4O620BE-01 
7.93?27?8E-nl 
f>.5366643F-ol 
4.9905163E-O1 
3.4452974E-01 
2.0R23R51E-01 
9.4R236?RE-o2 
2.410  19o3E-('2 
1.59R93?3E-o6 


Rise  tim 
CURREMTIR 
2. 4?977?0E-02 
9.5196BR7E-0? 
2.0575264E-01 
7.4513438E-01 
4.99687R5E-01 
6. 5427133E-01 
7 , 9374 1 QoE -0 1 
04436O0E-O  1 
9.755094RE-01 
9.9999960F-01 

9. 75507?] E -P] 
04431 6RE-01 
7.9373595E-01 
6. 5426434E-Q1 
4. 9968020E-01 
3. 4512740E-O 1 
?• 0574670 E -01 
9.5192574E-0? 

2, 4295458E-0? 

7 . 994931 2E -0  7 


-  *  1.000 OE-OR 

DFCRP 

1 .9356304E-02 
3.68953B5E-02 
5.0820006E-02 
5.9766037E-02 
6.28570R1F-02 
5.9790322E-02 
5.0R66197E-02 
3.69589R8F-0? 

1 .9431030F-02 
4.6181547E-07 
1.9431908E-02 
-3.6959706E-02 
-5.0R66740E-02 
-5.9790607E-02 
-6.2B570B0E-02 

•5.9765751 E"02 
-5.0Rl9462r-02 
-3.6R9463HE-02 
•1.9355425e~02 
7.99439]4E-05 
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TABLE  II  (continued) 


VAUJF  Of 

VALUE,  of 
TI'<f  I  Mr  P£M£  ,'jT 
5.o"OooonE-'’7 

1 .0O00000E-'l#’ 

1. snoooooE-nb 

2.  nooooooF-o& 
?.sooonooE-06 

3.  OO00000E-(i6 
3.sooooooE-o6 
4.oooooooe-o6 

4.6n00000e-o6 
S.OOOOOOOE'-Ob 
5.SOOOOOOE-06 
S.OOOOOOOE-Ob 
6.S000O00E-06 
7.OO00000F -0  6 

7. ci000000t“06 
8.0"00000E-O6 

8 .  c>O00O00t«06 

9.  pnoooooE-r 6 
9.t>i'00000E-06 
1  .onooOOOE-05 


«C  =  1.0000E-09 

Hk  =  loono.oooo 

CtipPFNTlO 
2.44912P3E-0? 
9.656684] E-o2 

2  a  5626t)P4E— 0  1 

3.4573?p4E-ol 
S.0031M2E-01 
6.54flf)Qp3E-rl 
7.94?SoF6E-ol 
9.  04805^9E-fll 
9. 7S70179E-PI 
9,OQ94Qh0E-0i 
9. 7B312R9E-01 
9*n406?OHC-Ol 
7.93??7?HE-0l 
6.S366643E-C 1 

4.  99051  b3E-f.l 
3.44B2O74E-01 
2*0523flSlE-M 
9.48?36?8t-o2 
2.4101Q03E-C2 
1.S9893?3E-(i6 


WISE  TIM 
CJRPENJTIR 
2.4394314E-0? 

5381 204E-02 
2.0600663E-O1 
3. 464331 5E-01 
6.0000184E-0I 
b, 54S7034E-01 
7.9399634E-01 
9.0462095E-O1 
9.75606R2E-01 
9, 99999rqE-01 
9, 7541 023E-0 1 
9.0424704E-01 
7.934R173E-01 
6.5396645E-01 
4.9936592E-01 
3.4482851E-01 
2.0649249E-01 


r  =  l.OOOOE-OS 
decrm 

9.6969509F-03 
1.8463694F-02 
2.5421839E-02 
2.9889149E-02 
3.14r8565E-02 
2.9889078E-02 
2.5421503E-02 
1.84635Q7E-02 
9.6967313E-03 
-1 .99859R7E-05 
-9.7347453E-03 
-1.8495839E-02 
-2.5444985E-02 
-2.99014HE-O2 
-3.1428540E-02 
-2.9e767f,8E-02 
-2.5398115E-0? 


9.5C07941E-02.-1.8431332F-02 
2.4198493E-0?  -9.6589213E-03 
1.0016619E-06  5.9727Q35E-05 


VAL'JF  OF 
VALUE  o£ 

T I MF  IMrREMFMT 
6. o000000E*o7 

1 .on0nn00E-O6 

1 .5000000? -06 

2. o^onooOE-06 
?.c;^nooonE-o6 
3.  oooopooE-i'6 
3.S000000E-06 

4. OOOOOOOE-06 

4.60000001-06 
5. OOOOOOOE-06 
5.6000000E-06 

6. OOQOOOOE-06 
6.SO0O000E-O6 

7.  nOOOOOOE -06 
7.6000000E-C6 

8.  nnOOOOOE-06 

8 .6000000E -06 

9.  OOOOOOOE-06 

9.6O00000E-06 
1  .OO00000E-O5 


RC  «  2.0000E-1 

XK  r  6OO00.O00 
C"frf;htiq 

2.449] ?P3E-o2 
9.6666841 E— 02 
2. n626084E-nl 

3«4573?P4E-0l 
10316i2E-ol 
6  *  6486923E-0 1 
7.942S066E-nl 
9. n4R0669E-n 1 
9.7670379E-01 
9.99990f,oE-ol 
9.7531?69E-ol 

9 . 0  4  n  6  2  0  8  E  -  0  1 

7.93227?8E-ol 
6.S366643E-P  1 
4.9906] 63E-01 
3.4462974E-01 
2  •  0923851 E-o 1 
9.482362ME-02 
2.4101903E-02 
1 .69P97P3E-06 


0  RI6E  T I 

fJRREMTIR 

2.4471859E-02 

9.5528R88E-02 

2.062099HE-01 

3.4567225E-0J 

6.00253?  7E -0 1 
6.54«0946E-ol 
7.9419973E-01 
9.0476869E-O1 
9.7568442E-01 
9,99999f,7E-01 
9.7533239E-03 
9.0409910E-01 
7. 9327819E-01 
6.5372625E-01 
4.9911449E-01 
3.4458948E-01 
2.0528929E-O1 
O.4860465E-02 
2,4121 191 E-0? 
1 .44787O1E-06 


•'  *  1.0000E-05 

f)ECRP 

1.9423969E-03 

3.6952969E-03 

5.0861866E-03 

5.978H073E-03 

6.2P57134E-03 
5.9768388E-03 
5.0824424E-03 
3.6901438E-03 
1.9363398E-03 
-7.1580128E-06 
- 1 » 949954 6E- 03 
-3.7Q17236E-03 
-5.0908527F-03 
-5.9812558E-03 
-6.2857044E-03 
-5.9743731E-03 
-5.0777616E-03 
-3.6837064E-03 
-1.9287765E-03 
1.51Q6218E-05 
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SECTION  VI 


SUMMARY  DF  EXPERIMENTAL  WORK  ON  PROJECT 

The  work  done  on  the  crystal  phase  of  the  contract  has 
included  both  analytical  and  experimental  studies  involving 
both  static  and  dynamic  tests.  This  section  of  this  report 
describes  the  work  done,  and  gives  data,  results,  and  con¬ 
clusions  where  applicable. 

STATIC  TESTS 

Static  load  tests  were  conducted  on  several  different 
materials  in  order  to  measure  the  piezoelectric  charge  co¬ 
efficient,  d^.  Information  on  these  tests  has  been  sum¬ 
marized  in  earlier  quarterly  progress  reports,  and  this  will 
be  repeated  where  necessary  and  supplemented. 

The  loads  were  applied  to  the  crystals  using  tv/o  dif¬ 
ferent  machines.  Loads  less  than  5000  lb.  were  applied  in 
a  Dillon  tensile  testing  machine  (using  a  compression  cage) , 
and  loads  greater  than  this  were  applied  with  a  Reihle  uni¬ 
versal  testing  machine.  Figure  6  shows  the  hardware  and 
method  used  to  apply  loads  to  the  crystals. 

Recording  of  the  data  during  loading  was  done  with  an 
X-Y  recorder  for  convenience.  Since  the  piezoelectric 
charge  coefficient  has  dimensions  of  charge  per  unit  of 
force,  it  was  necessary  to  develop  measuring  techniques  and 
circuitry  to  record  these  quantities . 

In  the  low  force  ranges,  measurement  of  force  was  ac¬ 
complished  with  a  Kistler  Corporation  load  cell  and  charge 
amplifier.  In  the  higher  force  ranges,  measurement  was  ac¬ 
complished  with  a  special  load  cell  designed  and  constructed 
at  The  University  of  New  Mexico  (UNM) .  Figure  7  describes 
this  transducer,  and  Figure  8  gives  the  calibration  data  and 
curves  for  it.  In  obtaining  this  calibration  curve,  the  gages 
on  the  transducer  were  connected  in  the  customary  bridge  circuit. 
The  bridge  current  was  adjusted  to  give  the  output  as  indicated 
by  the  specified  calibrated  resistor  and  the  corresponding  de¬ 
flections  . 
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CROSS  SECTION  DETAIL  OF  SPECIMEN  HOLDER 


Figure  6.  Hardware  for  Static  Tests  on  Ferroelectric  Crystal 


Impact  Piace 
0.285  lb. 


Eight  strain  Gages 
on  this  Section 
Type  A- Ik 
R=500-3  ohm 
F=2.05t2% 


MATERIAL:  MILD  STEEL 
MODULUS = 2 9xlOfo  psi 
SPRING  CONSTANT=22xlO^  lb/in 
NATURAL  FREQUENCY=20, 500  Hz 


Figure  7.  High  Range  Force  Transducer 


Figure  9  gives  a  general  description  of  the  recording 
technique  for  these  tests.  The  circuitry  for  recording  the 
charge  developed  is  shown  here.  The  series  capacitor,  Cc, 
charged  with  an  appropriate  voltage  from  the  power  supply  was 
used  to  calibrate  the  circuit.  The  shunt  capacitor,  C„,  was 
used  to  load  the  circuit  in  order  that  the  cathode  follower 
would  not  be  overdriven,  and  as  an  adjustment  on  the  signal 
into  the  X-Y  recorder  to  have  adequate  indicator  deflection 
over  a  wide  range  of  charge. 

Figures  10,  11,  and  12  give  the  results  of  the  calibra¬ 
tion  tests  on  this  circuitry.  Three  curves  are  given  to  cover 
different  ranges  of  charge.  The  curves  on  Figure  12  illus¬ 
trate  two  points:  1)  Repeatability  of  the  calibration  curve 
after  a  time  lapse  of  several  days,  and  2)  the  fact  that  small 
deviation  in  Cc  has  a  negligible  effect  on  the  calibration 

curve  as  long  as  C  >>>C  . 

s  c 

The  test  equipment  was  used  initially  with  a  quartz  cry¬ 
stal.  The  properties  of  quartz  are  well  known,  and  the  re¬ 
sults  of  th_s  test  were  used  as  a  check  on  the  measuring  sys¬ 
tem.  The  1/2-inch  diameter  quartz  disc  was  loaded  tc  4300  lb. 
and  unloaded.  The  charge-force  relation  was  plotted  on  the 
X-Y  recorder.  The  results  showed  a  linear  relationship  with 
essentially  no  hysteresis.  The  slope  of  the  curve  was  com¬ 
puted  to  be  1.93  pcb/N,  which  compares  quite  favorably  with 
the  d^  constant  for  quartz  shown  on  page  4.  It  was  con¬ 
cluded  that  the  test  equipment  v;as  adequate  for  making  the 
desired  measurements,  and  two  samples  each  of  HST-41,  PZT-4 
and  PST- 5  were  tested. 

One  series  of  tests  on  HST-41  were  conducted  primarily 
on  one  sample.  Table  III  describes  these  tests  and  gives 
experimental  values  for  Runs  1  through  4  indicate  that 

the  d32  constant  under  low  stress  values  is  repeatable,  but 
an  upward  shift  is  indicated  even  for  these  low  values  of 
stress.  A  typical  curve  from  one  of  these  tests  is  shown  in 
Figure  13.  This  curve  shows  a  small  hysteresis  loop  and  some 
nonlinearity  during  the  loading  phase  of  the  curve . 
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Figure  9.  Circuitry  for  Static  Tests  on  Ferroelectric  Crystals 


S3iOA-riTW  /^nd^no 
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TABLE  III 


RESULTS  OF  STATIC  LOAD  TESTS  ON  HST-41 


Test  Purpose 

No. 


Peak  force 


1 


lb. 


pcb/N 


1 

2 

3 

4 

Repeatability  of 

for  low  force. 

19.3 

20.6 

18.3 

20.9 

2245 

2261 

2368 

2610 

5 

Change  of  charge- 

1218 

5670 

6 

force  behavior  when 

19.4 

1034 

7 

cycled  to  alternating 

1220 

2945 

8 

high  and  low  force. 

20.2 

1487 

9 

2215 

5260 

10 

24.0 

1249 

11 

2215 

3930 

12 

through 

22 

23 

Change  of  charge- 

force  behavior  when 

cycled  to  the  same 

force. 

2000 

24 

Check  for  recovery  of 

2222 

2430 

25 

crystal  after  30  days 

without  testing. 

2262 

2122 

26 

10,000 

1690 

27 

9,150 

660 

28 

9,100 

450 

29 

9,050 

398 

1  2 
The  crystal  tested  was  1/2  in.  dia . ,  for  an  area  of  .196  in.  . 

Therefore,  approximate  stress  values  can  be  computed  by 

multiplying  the  force  value  by  5. 

2 

The  values  given  here  were  computed  from  the  final  readings 
of  the  charge  and  force,  and  as  such  do  not  represent  any 
slope  associated  with  the  curve  where  nonlinearities  exist. 


50 


50,000 


Stress,  psi 


Runs  5  through  11  were  conducted  to  study  the  behavior  of 
the  crystal  when  cycled  with  alternating  high  and  low  force 
values .  Figure  14  shows  the  curve  obtained  with  the  first  higjn 
force  value.  The  behavior  for  the  low  force  tests  was  similar 
to  that  shown  in  Figure  13,  except  the  hysteresis  loop  was 
slightly  larger.  This  series  of  tests  covered  a  time  interval 
of  approximately  1  hour.  Repeated  tests  to  this  stress  level 
do  show  a  decrease  in  the  charge  developed. 

Runs  12  through  22  were  conducted  one  right  after  another, 
each  to  the  same  load.  Figure  15  shows  the  results  of  this 
series;  as  before,  a  degradation  in  the  charge  developed  was 
noted.  Run  23  was  made  after  the  crystal  "rested"  for  30  min¬ 
utes.  A  small  amount  of  recovery  is  noted. 

The  crystal  was  then  allowed  to  remain  unloaded  for  30 
days.  Re-testing  then  indicated  a  further  recovery  of  charge 
producing  capability,  but  not  to  the  original  level. 

Tests  26  through  29  were  conducted  to  very  high  stress 
levels  —  approximately  50,000  psi.  Figure  16  shows  the  re¬ 
sults  of  test  26.  In  it  noted  that  at  the  high  stress  levels, 
very  little  additional  charge  is  developed  as  the  force  in¬ 
creases.  Repeated  tests  to  the  same  stress  level  show  a  signi¬ 
ficant  decrease  in  "d^-j." 

Table  IV  shows  the  results  of  the  tests  on  PZT-4.  Figure 
17  shows  the  results  of  test  2  on  sample  2  and  test  4  on  sam¬ 
ple  3.  Figure  18  shows  results  of  tests  at  stress  levels  of 
approximately  5000  psi.,  and  Figure  19  shows  results  of  several 
tests  at  the  50,000  psi.  level. 

A  review  of  the  above  mentioned  table  and  figures  gives 
the  following  indications  regarding  the  behavior  of  PZT-4. 

1)  There  may  be  considerable  variation  in  d^  among 
samples  of  this  material. 

2)  The  force-charge  relation  is  closer  to  being  linear 
than  for  HST-41  up  through  5000  psi,  and  there  is 
less  hysteresis. 

3)  At  stress  levels  of  50,000  psi,  the  curve  is  be¬ 
ginning  to  flatten  out. 
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Figure  14.  Charge- 


Sample 
Test  26 


qpd  ' s6 jpq3 
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Charge-Force  Curve  for  HST-41  to  Compressive  Stress  of  50,000  psi 


TABLE  IV 


RESULTS  ON  PZT-4 

Test  No.  SAMPLE  2  SAMPLE  3 


Stress 

psi 

^33 1 
pcb/lb 

Stress 

psi 

d33 

pcb/lb 

1 

180 

1042 

— 

— 

2 

174 

988 

10  3 

1052 

3 

188 

1046 

159 

1701 

4 

195 

1045 

185 

1456 

5 

— 

— 

184 

1487 

6 

5100 

1060 

186 

1392 

7 

5100 

980 

5370 

1694 

8 

5340 

1575 

9 

10,300 

2215 

10 

10,200 

2120 

11 

11,300 

2620 

12 

51,000 

2990 

13 

46,900 

3100 

14 

46,900 

2450 

15 

46,900 

1795 

1The  values  given  here  were  computed  from  the  final  readings 
of  the  charge  and  fores,  and  as  such  do  not  represent  any 
slope  associated  with  the  curve  whsre  nonlinearities  exist. 

2 

There  was  a  30-day  wait  between  test  10  and  11  on  crystal  3. 


56 


50000 


qod  'efijBqo 


57 


20  40  60  60  100  120  140  160 

Stress,  psi 

Figure  17.  Charge-Force  Relations  for  PZT-4  at  Low-Stress  Levels 


Charge-Force  Relations  for  PZT-4  at  Stress  Levels 


A  summary  of  the  results  of  the  tests  on  PZT-5  is  shown 
in  Table  V.  Figures  20,  21,  and  22  show  typical  results  from 
these  tests. 

A  second  series  of  tests  was  conducted  to  determine  a 
basis  for  obtaining  a  nonlinear  and  uniaxial  stress  charge 
coefficient  for  HST-41  at  high  stress  levels.  Fcr  this  series 
a  new  crystal  was  used  for  each  test.  In  order  to  insure  a 
state  of  uniaxial  stress,  a  slightly  different  loading  tech¬ 
nique  was  used.  The  crystals  were  loaded  through  two  bars 
having  the  same  diameter  as  the  crystal,  as  shown  in  Figure  23. 
This  method  was  used  since  the  boundary  conditions  on  the  cry¬ 
stals  in  the  static  tests  would  approximate  those  for  the  split 
Hopkinson  bar  tests  which  would  be  conducted  in  the  dynamic 
tests.  Two  different  materials,  steel  and  aluminum,  were  used 
for  bars  though  which  the  load  was  applied,  and  two  different 
thicknesses  of  crystal  were  used,  —  0.50  inch  and  .100  inch. 
Again  this  was  done  to  duplicate  the  later  dynamic  tests.  The 
various  combinations  of  bar  materials  and  crystals  are  shown 
below. 


Bar  Material 

Crystal  Thickness 
(inches) 

Test  1 

Aluminum 

.100 

Test  2 

Aluminum 

.100 

Test  3 

Aluminum 

.100 

Test  4 

Steel 

.100 

Test  5 

Steel 

.100 

Test  6 

Steel 

.050 

Tests  were  conducted  on  a  Riehle  Universal  Testing  machine 
using  the  same  Kistler  load  cell  in  the  earlier  test  series  to 
monitor  the  load.  Figure  23  also  shows  the  load  washer  and 
test  specimen.  Recording  was  done  as  before,  using  a  large 
shunt  capacitor  to  prevent  the  test  crystal  from  overdriving 
the  Kistler  charge  amplifiers.  A  second  Kistler  charge  ampli¬ 
fier  was  used  in  conjunction  with  the  load  cell.  Both  charge 
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TABLE  V 

RESULTS  ON  PZT-5 

1 


Test  No. 

Stress 

psi 

d33 

pcb/lb 

1 

79 

1995 

2 

100 

2300 

3 

113 

2228 

4 

126 

1859 

5 

5810 

4100 

6 

6140 

7590 

7 

5810 

6060 

8 

5850 

2080 

9 

412 

780 

10 

475 

481 

11 

441 

839 

12 

371 

1210 

132 

2180 

4959 

14 

50,900 

997 

16 

35,203 

348 

The  values  given  here  were  computed  from  the  final  readings 
of  the  charge  and  force,  and  as  such  do  not  represent  any 
slope  associated  with  the  curve  where  nonlinearities  exist. 

2 

Crystal  given  a  30 -day  rest  prior  to  this  test. 

^Crystal  fractured  at  this  load. 
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Figure  20.  Charge-Force  Curve  fo] 


Applied  Load 


Figure  23.  Test  Fixture  Used  for  Determining  the 

Static  Charge-Force  Relation  for  HST-41 
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Charge,  microcoulombs 


amplifiers  then  drove  the  Moseley  X-Y  plotter.  Calibration 
was  done  directly  on  the  recording  circuits  by  using  a  DC  power 
supply  and  a  series  capacitor  to  simulate  a  charge  on  the  cry¬ 
stal. 

The  crystals  were  loaded  through  the  aluminum  bars  to 
about  20,000  psi,  and  to  about  40,000  psi  through  the  steel 
bars.  Data  were  taken  on  the  initial  loading  only.  The 
curves  for  all  of  the  six  tests  fell  within  a  band  which  was 
quite  narrow,  possibly  within  the  range  of  variation  of  speci- 
ment  properties.  Therefore,  only  one  experimentally -determined 
curve  is  shown  for  this  test  searies.  Test  5  was  taken  as 
representative  of  all  data  and  is  shown  in  Figure  24.  It 
was  used  for  determining  the  charge  stress  functions  as  needed 
in  the  numerical  work. 

Two  factors  concerning  these  tests  should  be  noted.  One 
is  that  the  zero  load  charge  coefficient  determined  experi¬ 
mentally  in  this  test  series  agrees  very  closely  with  that 
given  by  the  manufacturer.  However,  at  a  load  of  about 
16,000  psi,  the  charge  coefficient  has  increased  to  four  times 
that  at  zero  load.  At  very  high  loads  then  the  charge  coef¬ 
ficient  drops  off  to  less  than  10%  of  the  zero  load  value. 

This  is  expected  since,  as  the  load  increases,  the  charge  gene¬ 
rated  approaches  the  maximum  available  charge  given  by  Equa¬ 
tion  12.  For  these  tests  the  maximum  charge  approaches 
about  32  ycb.  This  charge  coefficient-average  stress  rela¬ 
tion  is  later  used  in  the  reduction  of  the  data  from  the 
Hopkinson  bar  tests. 

The  second  factor  is  that  the  output  of  the  new  crystals 
was  substantially  higher  than  for  the  older  crystals.  Hence 
it  is  suspected  that  the  ferroelectric  materials  do  lose  some 
of  their  charge  potential  with  age.  It  is  not  believed  that 
the  nonlinear  charge  coefficient  function  determined  for  the 
new  crystals  would  be  valid  for  the  older  crystals. 

DYNAMIC  TESTS— SHOCK  TUBE 

A  considerable  amount  of  work  has  been  done  in  earlier 
investigations  on  response  of  crystals  to  shock  wave  inputs. 
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This  type  of  input  requires  that  the  initial  contact  with 
the  face  of  the  crystal  be  plane,  and  achieving  this  condition 
experimentally  can  be  something  of  a  problem.  The  air  driven 
shock  tube  was  used  as  the  initial  test  vehicle  on  the  crystals 
since  the  loading  condition  satisfied  this  requirement,  the 
tests  could  be  run  quickly,  and  the  tests  were  nondestructive 
(in  the  sense  the  crystal  was  not  broken  during  each  test) . 
These  factors  were  of  importance  in  the  first  series  of  tests. 
The  tests  were  undertaken  to  accomplish  the  following:  (1)  to 
check  the  theory  associated  with  use  of  the  piezoelectric  ef¬ 
fect  as  an  energy  generator,  as  presented  in  Section  III;  (2) 
as  an  initial  probe  into  the  problems  involved  in  the  circuitry 
and  measuring  techniques;  and  (3)  in  order  to  determine  the 
effects  of  the  measuring  instruments  in  crystal  circuits. 
Figures  25  and  26  describe  the  hardware  and  recording  techni¬ 
ques  used  in  these  tests . 

Materials  tested  were  quartz  and  HST-41.  The  recording 
circuit  used  is  shown  in  Figure  27.  The  termination  resistor 
was  50  ohm  on  all  tests.  For  one  series  of  tests,  which  was 
run  only  on  HST-41,  the  effect  of  increasing  the  shock  tube 
driving  pressure  was  studied.  The  load  resistance,  RL,  was 
zero  for  these  tests.  Figure  28  shows  typical  data  for  this 
series  and  a  sketch  of  the  data  for  all  of  the  tests  conducted. 
It  is  observed  that  the  voltage  drop  across  RTERM  increases 
with  driving  pressure,  as  expected.  It  should  be  pointed  out 
that  the  driving  pressure  does  not  represent  the  pressure  ap¬ 
plied  to  the  crystal,  but  is  the  charge  pressure  in  the  shock 
tube. 

Another  series  of  tests  run  on  HST-41  only  involved  a 
study  of  the  effect  of  crystal  thickness.  Tests  were  run  on 
crystals  varying  in  thickness  from  0.02  to  0.25  inch.  The  re¬ 
sults  of  these  tests  are  shown  in  Figure  29.  Theory  indicates 
that  the  current  generated  will  vary  inversely  with  thickness 
(Eq.  4,  p.  3).  Figure  29  shows  that  while  the  voltage  drop 
across  R-TERM  is  dependent  on  the  thickness.  Equation  4  is  not 
applicable  for  the  inputs  of  these  tests.  It  is  believed 
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Figure  25.  Hardware  for  Shock  Tube  Tests 
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Figure  26 .  Recording  Technique  Circuitry  for  Shock  Tube  Tests  on  Crystals 
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Figure  27.  Recording  Circuit  Diagram 
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Figure  28.  Typical  Data  and  Results  for  Shock  Tube  Tests 

on  HST-41  Crystal  (.50  in.  dia.  x  0.10  in.  thick) 
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Figure  29.  Effect  of  Crystal  Thickness  on  Output  for  HST-41 
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that  the  deviation  from  the  theory  is  the  result  of  the  in¬ 
herent  lack  of  planeness  of  the  force  input  to  the  crystal. 

A  study  of  the  effect  of  load  resistance,  R^,  on  current 
and  energy  generated  was  made  for  quartz  and  HST-41.  The  load 
resistance  was  varied  from  zero  to  100,000  ohm  and  driving  pres¬ 
sure  was  50  psig  for  all  tests.  Figures  30  and  31  show  typical 
data  on  quartz  and  HST-41,  respectively.  Figure  32  summarizes 
the  pulse  signatures  for  all  runs  for  quartz,  and  similarly. 
Figure  33  for  HST-41. 

The  data  obtained  for  the  range  of  load  resistors  has  been 
used  in  further  studies  related  to  the  problem.  These  are  as 
follows  s 

1.  Study  of  the  effect  of  the  impedance  of  measur¬ 
ing  instruments  in  the  circuit  used  with  the 
crystals . 

2.  Study  of  the  effect  of  input  pulse  on  the  re¬ 
sponse  of  the  circuit. 

3.  Study  of  current  generation  -  specifically  to 
predict  currents  generated,  and  effect  of 

on  current. 

4.  Investigation  into  the  amounts  of  energy  gene¬ 
rated  for  these  tests. 

Details  of  these  studies  will  be  discussed  in  the  above  order. 

Items  1  and  2  above  were  investigated  with  the  same  ap¬ 
proach.  The  circuit  diagram  of  Figure  27  includes  essentially 
all  the  elements  involved.  Analysis  of  this  circuit  would  be 
quite  difficult.  This  circuit  was  first  reduced  to  the  one 
labelled  "real"  in  Figure  34.  In  doing  this,  the  following 
assumptions  were  made: 

I.  Ry  »>  Rl 

2*  ^IST  <<:<  *1, 

3*  CSCOPE  >>  Cp  +  CDIST 

^TOTAL  p  DIST 
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Driving  Pressure  =  50  psig,  all  tests. 
Total  Capacitance  =  150  pfd,  all  tests. 


.2 


ysec/cm 


Figure  30.  Examples  of  Data  from  Shock  Tube  Tests 

on  Quartz  (0.50  in  dia.  x  0.10  in.  thick) 
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Driving  Pressure  =  50  psig,  all  tests. 
Total  Capacitance  =834  pfd,  all  tests. 
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Figure  31.  Examples  of  Data  from  Shock  Tube  Tests 

on  HST-41  (0.50  in.  dia.  x  0.10  in.  thick) 
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load  resistor,  R. ,  for  that 
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Figure  32.  Summary  of  Data  for  Quartz  Shock  Tube  Tests 
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IDEAL  CIRCUIT 


Figure  34.  Circuits  Used 


This  circuit  retains  the  principal  characteristics  of  the 
actual  circuit,  i.e.,  the  impedance  of  the  measuring  instru¬ 
ment,  and  will  permit  investigation  of  its  effect  on  the  mea¬ 
sured  crystal  output.  The  applicable  circuit  when  the  mea¬ 
suring  instrument  is  not  being  used  is  shown  as  the  "ideal 
circuit"  in  Figure  34. 

The  values  of  the  parameters  used  were  those  from  the 
test  series,  and  computations  were  made  using  a  piezoelectric 
charge  coefficient  for  both  quartz  and  HST-41.  Solutions  for 
the  required  currents  for  the  ideal  circuit  could  be  made 
using  analytical  techniques.  However,  an  analytical  solution 
to  the  differential  equations  for  the  real  circuit  would  be 
quite  messy,  and  a  numerical  technique  was  adopted. 

This  numerical  technique  made  use  of  the  computer  program 
"ECAP" ,  which  was  obtained  from  IBM.  Use  of  the  program  re¬ 
quires  insertion  of  the  desired  impedance  values  between  the 
appropriate  node  points  in  the  circuit.  The  program  then,  in 
essence,  sets  up  the  differential  equations  and  determines 
the  desired  currents  and  voltages. 

These  two  circuits  were  solved  for  both  a  step  input  and 
a  ramp  input.  Figures  35  and  36  show  the  results  for  the 
quartz  and  HST-41,  respectively.  A  study  of  these  figures  per¬ 
mits  the  following  conclusions  to  be  made: 

1.  For  the  value  of  used  in  this  study,  the  presence 
of  the  measuring  instrument  in  the  circuit  reduces 
the  current  through  the  load  resistor  a  noticeable 
amount.  The  amount  of  this  reduction  is  apparently 
dependent  upon  the  parallel  capacitance  in  the  cir¬ 
cuit,  since  the  percent  reduction  is  greater  for 
the  HST-41  than  the  quartz. 

2 .  The  ECAP  program  does  not  give  results  in  complete 
agreement  with  analytical  solutions  to  very  simple 
circuits,  as  evidenced  by  the  lack  of  agreement  be¬ 
tween  solutions  of  the  same  problem  by  each  techni¬ 
que.  This  difference  in  solutions  is  greater  for 
the  step  input  than  for  the  ramp.  However,  this 
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Figure  35.  Results  of  Circuit  Studies,  Quartz  Crystal 
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Figure  36.  Results  of  Circuit  Studies,  HST-41  Crystal 
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apparent  error  in  the  ECAP  method  is  considered  small 
enough  for  the  results  to  still  be  useful. 

3.  The  crystal  is  being  subjected  to  an  approximate 
ramp  input  by  the  shock  tube  having  a  rise  time  of 
about  0.4  ysec.  This  is  indicated  since  the  charac¬ 
teristics  of  the  recorded  signal  are  similar  to 
those  of  the  predicted  currents  for  the  ramp  input 
(see  Fig.  28) . 

The  theory  presented  earlier  specifies  the  current  gene¬ 
rated  when  a  quartz  crystal  is  struck  by  a  shock  wave.  The 
data  presented  gives  measured  currents  for  quartz  and  HZT-41. 
The  work  discussed  in  the  following  pages  is  that  of  computat- 
tion  of  current  generated  using  theory  and  comparison  with 
some  of  the  measured  values. 

These  computations  were  made  based  upon  the  piezoelectric 
theory  presented  in  Section  III.  The  parameters  used  in  the 
computations  were  as  follows: 


Quartz  HST-41 


Driving  Pressure,  P  ,  psig 

50 

50 

Pressure  on  Crystal, 

Pr,  psia 

37.2 

37.2 

Pressure  transmitted 

.495  P 

.66  P 

through  crystal 

r 

r 

Pulse  shape 

Ramp 

Ramp 

Rise  time,  ysec 

.38 

.38 

d33 

10.2  pcb/lb. 

1440  pcb/lb. 

PTERM'  ohmS 

50 

50 

V  Pfd 

130 

777 

CTotal'  P£d 

150 

797 

The  "ideal"  circuit  was  used  in 

the  computations 

£or  V 

Figure  37  shows  the  results 

of  these  computations  for 

quartz  and  Figure  38  for  HST-41. 

Points  from  the 

experiment- 

ally  determined  curves  presented  earlier  have  been  shown  on 
each  of  these  figures. 
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Figure  37.  Effect  of  Size  of  Load  Resistor  on  Current  Generated  by  Quartz 


Figure  38.  Effect  of  Size  of  Load  Resistor  on  Current  Generated  by  HST 


A  study  of  these  figures  indicates  the  following: 

1.  Theory  and  experiment  are  in  general  agreement 
for  low  stresses  used  in  this  phase  of  the  work. 

2.  The  agreement  between  the  computed  and  measured 
values  of  current  improves  as  the  size  of  the  load 
resistor  increases.  This  is  indicative  of  improper 
assumptions  and/or  measurement  regarding  circuit 
impedance.  As  the  resistance  becomes  large,  its 
impedance  in  the  circuit  then  overrides  the  other 
effects . 

The  computed  current  time  histories  presented  above  were 
used  to  determine  the  energy  dissipated  in  the  load  and  termi¬ 
nation  resistor.  The  significant  results  from  these  computa¬ 
tions  are  presented  in  Figures  39  and  40  for  quartz  and  HST-41, 
respectively.  These  figures  give  representative  values  of 
energy  generated  for  this  type  of  excitation  and  give  some 
indication  of  the  effect  of  circuit  resistance.  These  fig¬ 
ures  show  that  there  will  be  some  optimum  value  of  total  load 
resistance  for  maximum  energy  generation. 

DYNAMIC  TESTS— VACUUM-GUN  FACILITY 

Tests  on  the  crystals  utilizing  the  shock  tube  give  very 
low  applied  stresses.  The  vacuum-gun  facility  was  constructed 
for  which  the  function  was  to  test  crystals  under  conditions 
of  shock  excitation  similar  to  the  shock  tube,  but  to  higher 
stress  levels.  Figure  41  shows  a  schematic  of  the  gun.  The 
gun  is  termed  a  vacuum  gun  since  the  chamber  between  the  im¬ 
pact  plate  and  projectile  is  evacuated  before  each  shot. 

This  serves  two  purposes:  (1)  it  eliminates  the  air  cushion 
between  the  projectile  and  the  target,  thus  decreasing  pulse 
rise  time,  and  (2)  it  permits  the  atmospheric  pressure  to  be 
used  for  the  driving  force. 

Preliminary  tests  have  been  run  with  the  gun  and  impact 
velocities  up  to  50  fps  have  been  obtained.  The  most  effec¬ 
tive  way  of  controlling  the  impact  velocity  is  with  the  initial 
position  of  the  projectile  in  the  barrel,  rather  than  the 
amount  of  vacuum  drawn. 
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Figure  40.  Effect  of  Size  of  Load  Resistor  on  Energy  Generated  by  HST-41 


a  o  -h 


The  initial  tests  also  showed  that  the  fit  of  the  pro¬ 
jectile  in  the  barrel  was  not  sufficiently  close  to  give  the 
desired  planeness  of  impact.  An  indication  of  the  importance 
of  planeness  of  impact  is  given  by  this  expression  from  Ref. 
10: 

tan4>  =  tan0 
o 

where  0  =  angle  between  impacting  surfaces 
<J>  =  angle  of  wavefront  in  crystal 
UQ  =  impact  speed. 

A  study  of  this  expression  shows  that  the  tilt  angle  of  the 
impacting  surfaces  will  be  amplified  by  the  ratio  c/Uo.  For 
low  velocity  impact  of  the  magnitude  being  considered  for 
this  study,  this  amplification  will  be  large,  and  0  must  be 
minimized.  Considerable  time  was  spent  grinding  and  polish¬ 
ing  the  base,  and  machining  the  projectiles  to  obtain  a  0 
small  enough.  However,  with  the  equipment  available,  it  was 
not  possible  to  obtain  a  sufficiently  small  0  for  the  wave- 
front  developed  in  the  crystal  to  be  considered  plane. 

DYNAMIC  TESTS — SPLIT  HOPKINSON  PRESSURE  BAR 

Work  was  also  done  on  impulsive  loading  of  crystals  using 
a  split  Hopkinson  pressure  bar.  The  desirable  characteristics 
of  this  method  are  that  the  input  pulse  may  be  controlled  to 
a  greater  degree  than  with  the  other  methods  discussed,  infor¬ 
mation  on  the  pulse  signature  may  be  obtained  from  measure¬ 
ments  rather  than  from  the  crystal,  and  geometry  is  one  which 
could  possible  be  applied  to  a  proposed  system.  The  geometry 
used  is  shown  in  Figure  42.  The  impulse  to  the  input  rod  was 
applied  with  a  projectile  from  the  gas-gun  developed  for  this 
contract,  and  strain  gages  were  used  to  monitor  the  input 
pulse  to  the  crystal.  With  this  geometry,  impact  on  the  in¬ 
put  rod  does  not  have  to  be  as  nearly  plane  as  for  the  case 
when  the  projectile  strikes  the  crystal  directly.  When  non- 
planar  impact  occurs,  the  wavefront  approaches  the  planar 
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condition  as  it  moves  down  the  rod.  Since  this  geometry  is 
less  sensitive  to  the  alignment  of  the  impacting  surfaces  at 
impact,  it  is  somewhat  more  feasible  for  use  under  conditions 
not  found  outside  the  test  laboratory. 

Since  the  purpose  of  this  work  was  to  determine  energy 
levels  which  could  be  obtained  from  crystals  under  impact, 
the  work  for  this  type  of  loading  was  planned  as  follows : 

1.  Develop  through  the  use  of  appropriate  theory,  an 
analytical  method  for  predicting  current  and  energy 
levels  put  out  by  the  crystals  under  the  test  con¬ 
ditions  . 

2.  Conduct  several  experiments  to  determine  the  system 
input  and  output  during  actual  tests. 

3.  Verify  the  analytical  model  through  comparison  of 
measured  outputs  with  computed  values  for  the  same 
system  and  system  input. 

Once  the  analytical  model  had  been  verified,  it  would  have  been 
possible  to  make  studies  with  computer  simulation  in  order  to 
maximize  the  energy  output. 

The  work  done  in  pursuit  of  the  above  investigation  will 
be  presented  in  three  subsections.  The  first  subsection  will 
be  a  discussion  of  the  appropriate  theory  and  computational  pro¬ 
cedures  used  in  data  reduction.  The  second  subsection  will  de¬ 
scribe  the  experimental  phase  of  the  work,  and  the  last  sub¬ 
section  will  present  the  results. 

THEORY  AND  COMPUTATIONAL  PROCEDURES 

When  the  input  rod  is  struck  on  the  end  by  the  projectile, 
the  disturbance  initiated  propagates  down  the  rod  and  passes 
through  the  crystal.  Therefore,  the  theory  and  analytical  pro¬ 
cedures  required  are  wave  propagation,  current  generation,  and 
determination  of  energy  dissipated  by  current  through  a  resistor. 
This  subsection  is  divided  into  the  above  three  subsections. 

1)  Wave  Propagation 

The  governing  equations  for  the  propagation  of  the 
disturbance  in  the  rod  initiated  by  projectile  contact  may  be 
written  in  terms  of  the  field  equations 
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and  the  constitutive  relation 
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(14) 


(15) 


which  is  the  assumed  stress-deformation  law  for  the  materials 
in  which  the  wave  is  propagating.  Equation  (13)  results  from  the 
application  of  momentum  priniciples  to  an  element  of  the  rod, 
and  Equation  (14)  is  the  continuity  requirement. 

The  above  field  equations  were  derived  under  the  assump¬ 
tion  that  (1)  plane  cross  section  of  the  bar  initially  perpendi¬ 
cular  to  the  axis  remain  plane  and  perpendicular  to  the  axis 
at  all  times;  (2)  the  axial  stress  component  is  the  only  non¬ 
zero  stress;  (3)  the  stress  is  uniformly  distributed  over  the 
cross  section;  and  (4)  the  longitudinal  strain  is  given  by 
the  engineering  strain,  3u/3x,  and  it  is  the  only  non-zero  com¬ 
ponent  of  strain.  These  equations  are  therefore  based  upon  a 
one-dimensional  theory. 

The  solutions  to  the  field  equation  depends  upon  the  stress- 
strain  relation  used.  If  the  material  is  linearly  elastic,  we 
have 


f(°,e)  =  |  (15a) 

o 

and  if  it  is  an  elastic-plastic  solid, 

f(0'e)  =  FTiT  (15b) 

In  the  above  two  equations,  f(a,e)  is  the  function  in  Equation 
(15)  and  F(e)  is  the  static  stress-strain  function. 
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Equations  (13) ,  (14)  ,  and  (15) ,  when  solved  with  the  appro¬ 
priate  initial  condition  and  boundary  conditions  for  the  problem 
of  interest,  would  give  stress  and/or  strain  as  a  function  of 
time  at  any  cross  section  along  the  test  specimen.  The  stress 
information  at  the  cross  sections  located  at  the  two  faces  of 
the  crystal  permits  computation  of  the  current  through  the  cir¬ 
cuit  associated  with  the  crystal,  with  the  use  of  Equation  (11) . 

This  system  of  equations  is  hyperbolic  and  therefore  is 
soluble  by  graphical  integration  along  characteristics  in  the 
(x,t)  plane  even  when  the  constitutive  relation  is  of  a  form 
complex  enough  to  prevent  obtaining  a  close  form  analytical 
solution.  This  me the  1 ,  as  it  applies  to  wave  propagation  prob¬ 
lems,  has  been  discussed  in  detail  in  numerous  papers  (33,34), 
and  it  will  not  be  discussed  in  detail  here.  However,  a  gen¬ 
eral  description  of  the  numerical  procedure  used  will  be  given. 

The  significant  relations  upon  which  the  numerical  solu¬ 
tion  is  based  are  the  characteristics 


and  the  characteristic  relations 


(16) 


(17) 
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where 
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In  the  above  four  conditions,  the  superscript  identifies  the 
element  of  the  test  specimen  referred  to:  1  is  the  input  rod, 

2  is  the  crystal,  and  3  is  the  backup  rod.  Conditions  21  and 
22  are  the  result  of  the  crystal  being  securely  bonded  to  the 
input  and  backup  rod  during  the  test  duration  of  interest.  Con¬ 
ditions  21a  and  22a  are  the  result  of  a  force  balance  at  the 
interfaces  between  the  rods  and  the  crystals,  and  the  fact  that 
the  areas  of  all  three  elements  are  equal. 

Since  the  solution  in  closed  analytical  form  is  seldom 
known  for  quasi-linear  systems,  such  as  this  one,  a  finite  dif¬ 
ference  method  was  used  to  obtain  numerical  solutions  to  the 
problem.  Existence,  uniqueness  and  stability  of  the  solution 
of  the  initial  value  problem  for  the  hyperbolic  system  of  'n' 
quasi-linear  first-order  partial  differential  equations  have 
been  established  by  various  authors  in  the  literature. 

It  is  the  direct  replacement  of  the  derivative  along  the 
characteristics  that  form  the  basis  of  the  finite  difference 
method.  The  following  procedure  was  used  for  computing  the 
solution  of  the  initial  value  problem  of  Equations  (17)  and  (18) . 
The  respresentative  points  in  the  (x,t)  net  are  illustrated 
in  Figure  43.  Initially,  only  the  conditions  along  oig  char¬ 
acteristic  are  known.  This  characteristic  has  a  slope  of  c' 

2  3  1.° 

in  the  medium  1,  Cg  in  the  medium  2  and  Cg  (=Cg)  in  the  medium 

3.  Points  of  distance  At  apart  are  marked  off  along  the  t 

1  2 

axis,  and  points  of  distance  CgAt  in  medium  1,  CgAt  in  medium 
2  and  Cg  in  medium  3  are  marked  off  on  the  x  axis  as  shown  in 
Figure  43.  The  solution  for  the  characteristics  and  the  field 
variables  can  then  be  made  by  approximation  from  time  t  to 
time  t+At.  As  an  example,  the  procedure  for  finding  the  values 
of  a,  e  and  v  and  therefore  the  slopes  of  a  and  3  characteris¬ 
tics  along  the  CD  line  will  be  described,  assuming  that  the 
variables  on  the  line  AB  are  known  from  earlier  computations . 

From  the  values  of  a,  e  and  v  on  AB,  the  slope  of  the  a 
and  3  characteristics  along  AB  are  known  to  be  ±  /l/  f (a ,z) , 
and  this  is  shown  in  Figure  44 .  These  slopes  are  used  to  draw 
straight-line  approximations  to  the  characteristics  throughout 
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the  region  ABCD.  Through  each  of  the  points  Pq,  P^,  P •••» 

P  .  ,  P  on  CD,  the  characteristics  are  drawn  by  means  of  the 
n-l  n 

slopes  along  AB.  For  determining  the  field  variables  along  CD 
from  the  known  values  along  AB,  the  following  relations  derived 
from  Equations  (17)  and  (18)  are  used: 


Elastic  case: 


ep  "  £b  ‘  E  (ap-Jb)  =  0 


v  -  v, - —  (a  -a,)  =  0 

p  k  pco  p  k 


v  -  v  -  (a  -a  )  =  0 

p  s  pcQ  p  s' 


Plastic  case: 


e  -  e,  -  f,  (a  -a,  )  =  0 
p  b  b  p  b 


v  -  v  „  +  (a  -a  )  =  0 

p  £  pc£  p  £ 


v  -  v - (a  -a  )  =  0 

P  r  cr  P  r 


(23) 


(23a) 


From  the  known  values  of  o,  e  and  v  at  the  mesh  points  k 
and  b,  the  values  of  a and  v^  are  known  by  linear  inter¬ 
polation.  Similarly  from  the  known  values  of  a,  e  and  v  at  the 
mesh  points  b  and  s,  the  values  of  a^,  er  and  v  are  known  by 
linear  interpolation.  From  the  interpolated  values  a,  e  and  v 
at  £  and  r,  the  slopes  of  the  characteristics  through  £  and  r; 
i.e.,  c^  and  c are  known  from  Equation  (19). 

Equations  (23)  and/or  (23a)  can  then  be  solved  simultaneously, 

for  determination  of  the  values  of  a  ,  e  and  v  .  The  Scftne 

P  P  P 
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process  is  used  to  find  the  values  of  o,  £  and  v  at  all  mesh 
points  P  ,  P^,  . ..,  Pn_^/  pn  along  CD.  By  linear  interpola¬ 
tion,  the  values  of  a,  e  and  v  at  any  arbitrary  location  on 
CD  can  be  found  from  the  computed  values  of  a,  e  and  v  at  the 
neighboring  left  and  right  mesh  points  on  CD.  In  Equation  (23) 
and  (23a),  the  superscripts  identifying  the  medium  in  which  the 
point  lies  have  been  omitted  from  p,  c,  E,  a,  and  e,  for  gene¬ 
rality.  When  specific  computations  are  made,  the  correct 
superscript  must  be  used. 

Successive  applications  of  the  scheme  will  enable  the 
solution  to  be  advanced  in  time  as  long  as  necessary. 

2)  Current  Generation. 

The  short  circuit  current  from  a  crystal  loaded  with  a 
split  Hopkinson  pressure  bar  is  given  by  Equation  (5) ,  repeated 
here  for  convenience: 


i 


sc 


(t) 


d  Ac  (t)  Ac 
em  m 

- K - 


(5) 


The  difference  in  stress  on  the  front  and  rear  face  of  the 
crystal,  Ao^,  occurs  in  the  sp^it  Hopkinson  pressure  bar  test 
when  the  propagating  pulse  is  in  the  rise  or  decay  portion. 

As  previously  mentioned,  d  in  Equation  (5)  is  the  coef- 

em 

ficient  relating  charge  generated  to  force  on  the  crystal,  and 
as  such,  implies  the  relation  is  linear.  This  is  known  to  be 
the  case  for  quartz.*  However,  the  static  tests  have  shown 
this  not  to  be  true  for  the  ferroelectric  materials. 

The  current  generated  by  the  ferroelectric  materials  was 
determined  by  utilizing  a  nonlinear  relation  between  charge 
and  force.  A  cubic  piecewise  polynomial  was  fit  to  the  experi¬ 
mental  data  using  the  method  of  natural  cub  c  interpolating 
splines  (35) .  This  curve  fitting  technique  was  chosen  because 


*The  values  of  dem  and  c  used  in  this  equation  for  quartz  were 
plane  strain  values,  since  the  first  passage  of  a  strain  wave 
in  a  thin  disc  of  the  type  to  be  used  is  essentially  plane. 
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it  yields  smoother  approximations  of  the  first  derivative, 
which  is  the  quantity  of  interest  here.  The  curves  obtained 
fit  the  experimental  data  quite  well;  there  was,  in  fact,  no 
visible  deviation  from  the  original  data  curve.  See  Figure 
24,  which  shows  points  on  the  experimental  curve  which  were 
determined  from  the  curve  fit  routine.  Being  of  polynomial 
form  it  is  quite  easy  to  obtain  the  required  derivative  of  the 
approximating  curve.  Appendix  D  gives  a  listing  of  the  com¬ 
puter  program  used  for  this  curve  fit. 

Application  of  this  program  to  Equation  (5)  required  deter¬ 
mination  of  a  piezoelectric  charge  coefficient,  which  is  the 
slope  of  the  chaige-force  function,  at  a  desired  stress  value. 

At  a  given  instant  of  time  in  the  computational  procedure,  the 
stress  values  used  when  determining  the  charge  coefficient  was 
the  average  of  the  stress  on  the  faces  of  the  crystal  at  that 
time.  This  technique  would  be  correct  only  when  the  slope  of 
the  charge-force  function  is  essentially  constant  within  the 
stress  range  on  the  two  faces  of  the  crystal. 

3)  Energy  Release 

The  usable  energy  released  by  the  crystal  is  assumed  to  be 
that  which  is  dissipated  in  the  termination  and  load  resistors. 
This  may  be  computed  as 

t 

e  -  /  iR<t>  <Rl  +  RrEKM)dt  (24) 

This  equation  was  evaluated  with  the  use  of  a  digital  computer, 
and  the  program  used  is  given  as  Appendix  E. 

In  this  equation,  i  (t)  is  the  current  through  the  re- 
sistors.  This  current  will  be  approximately  the  short  circuit 
current;  the  difference  between  iD  and  i  being  due  to  the 
capacitance  loading  of  the  circuit.  This  affect  has  been  in¬ 
vestigated  and  the  results  discussed  in  Section  V. 
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4 )  Summary 

The  loading  of  the  piezoelectric-ferroelectric  crystals 
with  a  split  Hopkinson  pressure  bar  occurs  as  a  result  of  the 
passage  of  the  strain  wave  through  the  crystal.  The  stress¬ 
time  history  on  the  front  and  rear  surface  of  the  crystal  may 
be  determined  from  a  known  input  to  the  bar  and  a  one-dimen¬ 
sional  wave  propagation  theory .  The  appropriate  crystal  theory 
(with  either  linear  or  nonlinear  force-charge  relation)  may 
then  be  used  to  compute  the  short  circuit  current  from  the  cry¬ 
stal.  This  current  is  used  to  determine  the  electrical  energy 
available  for  release  in  a  resistor. 

The  discussion  just  completed  has  presented  the  appropriate 
theory  and  described  the  computational  procedures  used.  These 
computational  procedures  were  programmed  for  solution  on  a 
digital  computer,  and  listings  of  these  programs  are  given  as 
Appendices. 

Appendix  B  gives  the  program  which  computes  the  short  cir¬ 
cuit  current  and  energy  release  through  a  resistor  from  a  known 
strain-time  input  to  the  split  Hopkinson  pressure  bar.  This 
program  is  usable  for  a  crystal  which  has  either  a  linear  or 
nonlinear  force-charge  relation.  However,  it  is  restricted  to 
the  case  for  which  the  intensity  of  impact  is  such  that  only 
elastic  waves  propagate  in  the  pressure  bar.  Appendix  C  gives 
a  listing  of  the  program  which  is  applicable  when  the  impact 
on  the  bar  is  great  enough  to  exceed  the  elastic  limit  of  the 
material.  In  programming  for  this  case,  the  strain-rate  inde¬ 
pendent  theory  of  plastic  wave  propagation  was  used.  All  other 
aspects  of  this  program  are  similar  to  the  one  in  Appendix  B. 

EXPERIMENTAL  WORK 

The  actual  test  specimens  used  in  this  study  were  similar 
to  that  shown  in  Figure  42.  Tests  were  conducted  on  specimens 
for  which  the  bars  were  made  of  AISI  1095  steel  and  6061-T6 
aluminum.  The  input  rod  was  6  inches  long,  and  the  backup  rod 
was  4  inches  long. 

The  two  different  piezoelectric  materials,  x-cut  quartz 
and  HST-41  were  used  for  the  tests.  The  disks  were  1/2-inch 
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in  diameter,  and  two  thicknesses  of  each  material  were  used, 

.050  in.  and  .100  in.  The  x-cut  quartz  was  chosen  as  one  of 
the  materials  because  its  properties  are  well  known  and  there 
is  considerable  documentation  of  its  previous  usage  in  transient 
measurements  in  the  literature.  In  addition,  it  has  a  nearly 
constant  charge  coefficient  and  shows  little  sensitivity  to 
strain  in  other  than  the  x-direction,  hence  the  effect  of  vio¬ 
lating  the  assumptions  of  uniaxial  stress  or  strain  is  limited 
to  the  difference  in  mechanical  stiffness  between  the  two  cases 
about  10%)  which  is  well  understood.  It  was  felt  that  the  quartz 
tests  would  serve  as  a  rough  check  on  the  validity  of  the  com¬ 
puter  solution  of  the  wave  equations,  the  crystal  current  out¬ 
put  equation,  and  experimental  technique. 

The  HST-41  was  chosen  as  a  representative  ferroelectric. 

It  has  a  nonlinear  charge  coefficient,  and  is  a  readily  avail¬ 
able  material,  but  it  also  has  a  high  cross-sensitivity  (to 
both  manual  and  shear  strain) ,  and  a  high  electro-mechanical 
coupling  coefficient  (about  .3  to  .6).  The  nonlinearity  of  the 
charge  coefficient  made  it  necessary  to  use  the  nonlinear  charge 
force  relation  developed  for  this  purpose  from  static  tests. 

The  piezoelectric  discs  were  bonded  to  the  appropriate 
surfaces  of  the  bar  with  an  epoxy  adhesive .  The  incident  wave 
through  the  crystal  was  measured  with  foil  resistance  strain 
gages  mounted  in  opposed  pairs  on  the  surface  of  the  rod.  Two 
sets  of  gages  were  used.  One  set  was  mounted  2  in.  in  front 
of  the  tests  disc  to  serve  as  trigger  gages  for  the  recording 
oscilloscopes.  The  incident  strain  pulse  profile  was  measured 
by  a  second  pair  1  in.  in  front  of  the  disc.  For  recording, 
the  gages  were  wired  in  opposite  arms  of  a  4-arm  Wheatstone 
bridge  to  cancel  any  bending  strains. 

Recording  of  both  the  strain  gage  signals  and  the  crystal 
signal  was  done  on  Type  555  Tektronix  oscilloscopes  equipped 
with  Type  1A5  plug-ins.  This  system  had  a  bandwidth  from  DC 
to  30  MHz  (at  maximum  plug-in  gain) ,  which  was  more  than  ade¬ 
quate  to  faithfully  display  the  signals  from  the  sensors.  The 
signal  from  the  crystals  was  terminated  using  a  Tektronix 
50- ohm  terminator. 
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Tests  on  these  specimens  were  conducted  at  two  locations. 
The  low  speed  impact  tests  were  conducted  at  The  University  of 
New  Mexico.  The  projectile  was  driven  by  the  air  gun  con¬ 
structed  for  use  on  the  studies  supported  by  this  contract. 

The  high  speed  tests  were  conducted  at  Arnold  Engineering  De¬ 
velopment  Center,  Arnold  Air  Force  Station,  Tennessee. 

The  impact  speeds  for  the  tests  conducted  at  the  UNK 
facility  were  selected  to  cause  only  elastic  strains  in  the 
input  rod.  The  approximate  impact  speed  used  for  all  tests  was 
50  fps.  This  yielded  stress  levels  in  the  steel  input  rod  of 
about  40,000  psi,  and  about  20,000  psi  in  the  aluminum  roa. 

The  projectile  used  on  all  tests  was  a  1/2  in.  diameter  steel 
rod,  about  6  in.  long,  supported  in  the  air  gun  barrel  on  each 
end  by  nylo  sealing  rings.  The  sealing  rings  also  served  to 
keep  the  rod  aligned  in  the  gun  barrel.  This  type  of  projectile 
was  chosen  for  two  reasons.  First,  it  allowed  the  impact  prob¬ 
lem  to  be  simplified  to  that  of  impact  between  two  bars  of  the 
same  diameter,  and  second,  it  minimized  momentum  transfer  to 
the  test  rod,  thereby  reducing  the  damage  done  to  the  bar  and 
instrumentation  when  its  motion  was  arrested  after  impact. 

The  first  two  series  of  tests  were  conducted  with  the  end 
of  the  projectile  ground  flat.  However,  it  was  decided  that 
any  misalignment  between  the  test  rod  and  the  projectile  would 
produce  bending  stress  waves  in  the  rod.  While  this  would  not 
effect  the  output  of  the  strain  circuitry,  it  would  effect  the 
output  of  the  crystal.  Consequently,  subsequent  tests  were 
made  with  the  impact  end  of  the  projectile  ground  slightly 
spherical.  There  was,  however,  no  observable  difference  in 
either  case. 

One  entry  was  made  on  the  "G"  Range  at  AEDC,  to  determine 
crystal  output  for  very  high  impact  speeds.  The  "G"  range  at 
AEDC  is  a  two-state  powder-nitrogen  gun  approximately  100  ft. 
long. 

The  projectile  used  in  these  tests  was  a  2-1/2  inch  diam¬ 
eter  aluminum  slug,  2-1/2  inches  long.  The  recording  equip¬ 
ment  used  during  these  tests  was  essentially  the  same  as  that 
used  for  the  tests  at  UNM. 
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Two  types  of  specimens  were  used;  one  identical  to  that 
shown  in  Figure  42,  and  the  other  a  similar  geometry  in  which 
the  crystal  was  omitted.  The  rod  was  continuous,  and  strain 
gages  were  mounted  on  the  rod  at  the  crystal  station.  The  pur¬ 
pose  of  the  tests  on  the  second  type  of  specimen  was  to  deter¬ 
mine  the  amount  of  attenuation  of  the  propagating  wave  from  the 
monitoring  station  to  the  crystal  location. 

Six  tests  were  conducted  at  AEDC  on  these  types  of  speci¬ 
mens;  three  without  a  crystal  and  three  with.  Due  to  the  dif¬ 
ficulties  encountered  with  noise,  triggering,  and  gain  settings, 
usable  data  were  obtained  only  on  one  specimen  without  a  cry¬ 
stal  and  two  with  the  crystal.  Data  from  these  tests  is  shown 
in  Figures  45,  46,  47  and  48. 

Four  series  of  crystal  tests  were  conducted  on  the  air  gun 
at  UNM.  By  a  series  of  tests,  it  is  meant  one  entry  on  the  gun 
with  a  group  of  prepared  test  specimens.  This  procedure  was 
used  since  the  input  rod  and  back-up  rod  could  be  reused  after 
each  test.  So  preparing  a  specimen,  after  the  rods  were  once 
made,  consisted  of  cleaning  and  bonding  new  crystals  in  place. 

In  the  first  three  series  of  tests  the  ferroelectric  cry¬ 
stals  used  were  some  furnished  earlier  by  AFWL.  It  was  later 
discovered  that  their  exact  origin  and  age  could  not  be  ac¬ 
curately  determined,  and  data  had  indicated  this  to  be  a  factor. 
Consequently,  for  the  fourth  series  of  tests,  a  new  quantity 
of  HST-41  crystals  was  ordered.  The  results  of  the  tests  on 
these  crystals  is  considered  the  most  valid.  The  crystals  used 
in  the  static  determination  of  the  nonlinear  charge  coefficient 
were  also  from  this  new  set  of  HST-41  crystals. 

The  x-cut  quartz  specimens  used  in  all  tests  were,  however, 
of  known  origin.  Since  the  piezoelectric  properties  of  quartz 
are  natural,  they  are  not  affected  by  age  and  so  all  of  the 
tests  in  the  5  series  will  be  used. 

Table  VI  summarizes  the  tests  run  at  UNM.  Figures  49 
through  65  show  the  data  obtained  for  these  tests. 
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I. 


Test  Specimen 


TRIGGER  GAGE 


1/2"  DIAMETER  ,6061 -T  6 
/  ALUMINUM 


1 

'3 

4 

< - 

3  in. 

1  in 

3  in. 

4  in. 

All  measuring  gages  are  type  HE-111  (Budd  Co.),  1/16  in. 
gage  length  (F  =2.04,  R  =  120  ohm). 

II.  Projectile 

2h  inch  diameter  x  2%  inch  long  alu'.inum. 

Impact  speed:  8190  fps 

III.  Circuit  Information 

Each  measuring  gage  was  the  only  active  gage  in  a 
bridge  circuit. 


Gage 

1: 

Rcal‘ 

12,000 

ohms 

Scale 

factor: 

4950 

ye/cm. 

Gage 

3: 

Rcal  ’ 

20,000 

ohms 

Scale 

factor: 

2985 

ye/cm. 

Gage 

4: 

Rcal  ’ 

20,000 

ohms 

Scale 

factor: 

2985 

ye /cm. 

IV.  Scope  Data 


Upper  beam: 
Lower  beam: 
Sweep: 
Delay: 


Gage  1 

Radiometer (trigger) 
5ysec/cm — > 

5ysec 


Upper  beam:  Gage  3 
Lower  beam:  Gage  4 

Sweep:  5ysec/cm— » 
Delay:  25ysec 


Figure  45.  Data  for  Rod  Test  Without  Crystal  AEDC  Test  No. 

1204  (Data  from  other  scopes  on  the  above  gages 
were  obtained,  but  those  shown  were  the  best.) 
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All  measuring  gages  are  type  HE-111  (Budd  Co.),  1/16 
gage  length  (F  =2.04,  R  =  120 

II.  Projectile:  2h  inch  diameter  x  2h  inch  long  aluminum 
Impact  speed:  8380  fps 

III.  Circuit  Information 


Crystal  Circuit:  C1 

=  1210  pfd 

CII 

=  611  pfd 

See  Figure  46 

CP 

=  587  pfd 

Gage  circuit.  Gages  3 

and  4 :  R  . : 

cal 

20,000  ohm 

Scale  factor:  2985  yc/cm 


Scope  2  Upper  Beam:  Gage  3  Scope  3 
Lower  Beam:  Gage  4 

Sweep:  5usec/cm — > 
Delay:  lOysec 


Upper  Beam:  2v/cm 
Lowe:r  Beam:  5v/cm 

Sweep:  2ysec/cm 
Delay:  29ysec 


Figure  47.  Data  for  Rod  Test  with  Crystal.  AEDC  Test  No. 

1206  (Data  from  other  scopes  were  obtained,  but 
those  shown  were  used  in  data  reduction.) 
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I. 


Test  Specimen 


1/2"  DIAMETER, 606I-T6 


All  measuring  gages  are  type  HE-111  (Budd  Co.),  1/16  in. 
gage  length  (F  =  2.04,  R  =  120  ohm). 

II.  Projectile:  2%  in.  diameter  x  2%  in.  long  aluminum 
Impact  speed:  8150  fps 

III.  Circuit  Information 


Crystal  circuit: 


=  1175  pfd 
=  615  pfd 
=  600  pfd 


Gage  circuit:  Gages  3  and  4:  Rcal:  20,000  ohm 


IV.  Scope  Data 


Scale  factor:  2985  ye/cm 
Crystal  data 


Scope  2  Upper  beam:  Gage  3  Scope  3  Upper  beam:  2v/cm 

Lower  beam:  Gage  4  Lower  beam:  5v/cm 

Sweep:  5ysec/cm— »  Sweep:  2ysec/cm-* 

Delay:  lOysec  Delay:  29ysec 


Figure  48.  Data  for  Rod  Test  with  Crystal.  AEDC  Test  No. 

1207  (Data  from  other  scopes  were  obtained,  but 
those  shown  were  used  in  data  reduction.) 
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TABLE  VI 


SUMMARY  OF  TESTS  RUN  AT  UNM 


Series 

Test 

Crystal 
Mat 1 1 

Thick¬ 
ness  , 
Inch 

Rod 
Mat '1 

Data 

Obtaim 

1 

1 

PZT-4 

.100 

Steel 

Yes 

2 

Quartz 

.050 

Steel 

No 

3 

Quartz 

.050 

Alum 

Yes 

4 

Quartz 

.100 

Alum 

Yes 

5 

Quartz 

.100 

Steel 

Yes 

2 

1 

Quartz 

.050 

Steel 

No 

2 

Quartz 

.050 

Alum 

Yes 

3 

Quartz 

.100 

Steel 

Yes 

4 

Quartz 

.100 

Alum 

Yes 

5 

PZT-4 

.100 

Steel 

Yes 

6 

Pzt-4 

.100 

Steel 

Yes 

3 

1 

Quartz 

.050 

Steel 

No 

2 

HST-41 

.050 

Steel 

No 

3 

HST-41 

.050 

Steel 

Yes 

5 

Quartz 

.100 

Alum 

Yes 

6 

Quartz 

.050 

Alum 

Yes 

4 

1 

HST-41 

.100 

Steel 

Yes 

2 

HST-41 

.050 

Steel 

Yes 

3 

Quartz 

.050 

Alum 

Yes 

4 

HST-41 

.050 

Steel 

Yes 

5 

HST-41 

.100 

Steel 

Yes 

Comments 


gain  set  wrong 


bad  connection 


no  trigger 
late  trigger 


peak  current  lost 
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I.  Test  Specimen 


Test  Crystal 
.10  in.  thick  X  1/2  in. 
PZT-4  crystal 


7  in, 


2  in. 


Impact 


Trigger  gage  (2) 


Measuring  gage  (2) 

R  =  120  ohms  F  =  2.0 

9 


diam . 


l/2"inch  diameter  rods  -  1095  Steel 


II.  Projectile  and  Velocity  Information 

Projectile  length  -  6.7  in.  Time  increment  -  14216  usee, 
Firing  pressure  -  3.5  psig.  Impact  Velocity  -  39.?  fps. 

III.  Circuit  Information 

a.  Strain  circuit. 

Number  of  active  arms  -  2 
Bridge  current  -  30  milliamps 
Equivalent  strain  -  750  yin/in/division 
Rcal  =  20  K  ohms/2  divisions 


b.  Current  circuit. 

RL  =  50  ohms 

Rterm  =  50  ohms  (50- ohm  Tektronix  terminating 
resistor  was  used.) 

IV.  Scooe  Data 

Upper  trace  -  Strain 

Sweep  -  5  y sec/di v 

Gain  -  750  yin/in/div 
(equivalent) 

Lower  trace  -  current 

Sweep  -  5  ysec/div 

Gain  -  20  volts/div 

Delay  -  5  ysec  from 
strain  data 


Figure  49. 


Data  for  Crystal  Test  1,  Series  1 
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I.  Test  Specimen 


7  in . 


2  in. 


Test  Crystal 

.050  in.  thick  X  1/2  in.  diam. 
Quartz  crystal 

-+-i  ijM  J _ 


Impact 


Trigger  gage  (2)- 


(  ^Measuring  gage  (2) 

R  =  120  ohms  F  =  2.0 

3 

‘1/2-inch  diameter  rods  -  1095  Steel 


Projectile  and  Velocity  Information 

Projectile  length  -  6.7  in.  Time  increment  -  13697  ysec. 
Firing  pressure  -  3.5  psig.  Impact  Velocity  -  40.7  fps. 

Circuit  Information 

a.  Strain  circuit. 

Number  of  active  arms  -  2 
Bridge  current  -  30  milliamps 
Equivalent  strain  -  750  yin/in/division 
Rcal  =  20  K  ohms/2  divisions 

b.  Current  circuit. 

RL  =  50  ohms 

Rterm  =  50  ohnu~  (50-°hm  Tektronix  terminating 
resistor  was  used.) 

Scope  Data 


Upper  trace  -  Strain 

Sweep  -  5  ysec/div 

Gain  -  750  yin/in/div 
(equivalent) 


Lower  trace  -  current 

Sweep  -  5  ysec/div 

Gain  -  .2  volts/div 

Delay  -  5  ysec  from 
strain  data 


Figure  50.  Data  for  Crystal  Test  3,  Series  1 
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I.  Test  Specimen 


7  in. 


2  in. 


Test  Crystal 

.10  in.  thick  X  1/2  in.  diam. 
Quartz  crystal 

1  1  ^  J  _ 


Impac  t 


Trigger  gage  (2) 


Measuring  gage  (2) 

R  =  120  ohms  F  =  2.0 

g 


l/2-inch  diameter  rods  -  Aluminum 


II.  Projectile  and  Velocity  Information 

Projectile  length  -  6.7  in.  Time  increment  -  13643  ysec. 
Firing  pressure  -  3.5  psig.  Impact  Velocity  -  40.9  f ps . 

III.  Circuit  Information 

a.  Strain  circuit. 

Number  of  active  arms  -  2 
Bridge  current  -  30  milliamps 
Equivalent  strain  -  750  yin/in/division 
Rcal  =  20  K  ohms/2  divisions 

b.  Current  circuit. 

Rl  -  50  ohms 

Rterm  =  50  ohms  (50-ohm  Tektronix  terminating 
resistor  was  used.) 

IV.  Scope  Data 


Upper  trace  -  Strain 

Sweep  -  5  ysec/div 

Gain  -  750  yin/in/div 
(equivalent) 


Lower  trace  -  current 

S  eep  -  5  ysec/div 

Gain  -  .2  volts/div 

Delay  -  5  ysec  from 
strain  data 


Figure  51. 


Data  for  Crystal  Test  4,  Series  1 
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I. 


Test  Specimen 


7  in . 


2  in. 


r 


Test  Crystal 

•10  in.  thick  X  1/2  in.  diam. 
Quartz  crystal 

+1  i!M  7 


Impact 


Trigger  gage  [?) 


] 


Measuring  gage  (2) 

R  =  120  ohms  F  =2.0 
9 

1/2-inch  diameter  rods  -  1095  Steel 


II.  Projectile  and  Velocity  Information 

Projectile  length  -  6.7  in.  Time  increment  -  13915  psec. 
Firing  pressure  -  3.5  psig.  Impact  Velocity  -  40.1  fps. 

III.  Circuit  Information 

a.  Strain  circuit. 

Number  of  active  arms  -  2 
Bridge  current  -  30  milliamps 
Equivalent  strain  -  750  pin/in/division 
Rcal  =  20  K  ohms/2  divisions 

b.  Current  circuit. 

RL  =  50  ohms 

Rterm  =  ^0  °hms  (50-ohm  Tektronix  terminating 
resistor  was  used.) 

IV.  Scope  Data 

Upper  trace  -  Strain 

Sweep  -  5  psec/div 

Gain  -  750  pin/in/di’’ 
(equivalent) 

Lower  trace  -  current 

Sweep  -  5  pscc/div 

Gain  -  2  volts/d.,v 

Delay  -  5  psec  from 
strain  data 


Figure  52.  Data  for  Crystal  Test  5,  Series  1 
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I.  Test  Specimen 


2  in , 


2  in. 


Impact 


Test  Crystal 

.05  in.  thick  X  1/2  in.  diam. 
Quartz  crystal 

4-i  iJH 


Trigger  gage  (2) 


3 


Measuring  gage  (2) 

R  =  120  ohms  F  =  2.0 
9 


1/2-inch  diameter  rods  -  Aluminum 


II.  Projectile  and  Velocity  Information 

Projectile  length  -  6.7  in.  Time  increment  -  13595  ysec. 
Firing  pressure  -  3.5  psig.  Impact  Velocity  -  41.0  fps. 

III.  Circuit  Information 

a.  Strain  circuit. 

Number  of  active  arms  -  2 
Bridge  current  -  30  milliamps 
Equivalent  strain  -  750  yin/in/division 
Rcal  =  20  K  ohms/2  divisions 

b.  Current  circuit. 

Rl  =  50  ohms 

Rterm  =  50  oh,ns  (50~ohm  Tektronix  terminating 
resistor  was  used.) 

IV.  Scope  Data 


Upper  trace  -  Strain 

Sweep  -  5  ysec/div 

Gain  -  750  yin/in/div 
(equivalent) 


Lower  trace  -  current 

Sweep  -  5  ysec/div 

Gain  -  .2  volts/div 

Delay  -  5  ysec  from 
strain  data 


Figure  53.  Data  for  Crystal  Test  2,  Series  2 
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I.  Test  Specimen 


Test  Crystal 

.10  in.  thick  X  1/2  in.  diam. 
Quartz  crystal 


_ 1  in‘ - 

L— . . . 


Impact 


2  in. 


Trigger  gage  (2) 


Measuring  gage  (2) 

R  =  120  ohms  F  =  2.0 
9 


vl/2-inch  diameter  rods  -  1095  Steel 


II.  Projectile  and  Velocity  Information 

Projectile  length  -  6.7  in.  Time  increment  -  13775  ysec. 
Firing  pressure  -  3.5  psig.  Impact  Velocity  -  40.5  fps. 


III.  Circuit  Information 

a.  Strain  circuit. 

Number  of  active  arms  -  2 
Bridge  current  -  30  milliamps 
Equivalent  strain  -  750  yin/in/division 
Rcal  =  20  K  ohms/2  divisions 

b.  Current  circuit. 


Rl  =  50  ohms 

Rterm  ”  ")0  ollms  (50-ohm  Tektronix  terminating 
resistor  was  used.) 

IV.  Scope  Data 


Upper  trace  -  Strain 

Sweep  -  5  ysec/div 

Gain  -  750  yin/in/’div 
(equivalent) 

Lower  trace  -  current 

Sweep  -  5  ysec/div 

Gain  -  .2  volts/div 

Delay  -  5  ysec  from 
strain  data 


Figure  54.  Data  for  Crystal  Test  3,  Series  2 
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I.  Test  Specimen 


2  in. 


2  in . 


Test  Crystal 

.10  in.  thick  i*.  1/2  in.  diam. 
Quartz  crystal 
i  1  in. i  / 


L^1  J 


Impact 


Trigger  gage  (2'- 


(  ''-Measuring  gage  (2) 

R  =  120  ohms  F  =  2.0 

g 

'1/2-inch  diameter  rods  -  Aluminum 


II.  Projectile  and  Velocity  Information 

Projectile  length  -  <3.7  in.  Time  increment  -  13555  ysec. 
Firing  pressure  -  3.5  psig.  Impact  Velocity  -  41.1  fps. 

III.  Circuit  Information 

a.  Strain  circuit. 

Number  of  active  arms  -  2 
Bridge  current  -  30  milliamps 
Equivalent  strain  -  750  yin/ir>/division 
Rcal  =  20  K  ohms/2  divisions 

b.  Current  circuit. 

RL  =  50  ohms 

Rterm  =  ^  °^ins  (50-ohm  Tektronix  terminating 
resistor  was  used.) 


Scope  Data 


Upper  trace  -  Strain 

Sweep  -  5  ysec/div 

Gain  -  750  yin/in/div 
(equivalent) 


Lower  trace  -  current 

Sweep  -  5  ysec/div 

Gain  -  .2  volts/div 

Delay  -  5  ysec  from 
strain  data 


Figure  55.  Data  for  Crystal  Test  4,  Series  2 
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I.  Test  Specimen 


Test  Crystal 

.10  in.  thick  X  1/2  in.  diam. 


PZT-4  crystaj 


II „  Projectile  and  Velocity  Information 

Projectile  length  -  6.7  in.  Time  increment  -  13770  ysec. 
Firing  pressure  -  3.5  psig.  Impact  Velocity  -  40.5  fps. 


III. 


Circuit  Information 
a.  Strain  circuit. 


Number  of  active  arms  -  2 
Bridge  current  -  30  milliamps 
Equivalent  strain  -  750  yin/in/division 


R 


cal 


=  20  K  ohms/2  divisions 


b.  Current  circuit. 


Rl  =  50  ohms 

Rterm  =  ohms  (50 'ohm  Tektronix  terminating 
resistor  was  used.) 


Scope  Data 


Upper  trace  -  Strain 

Sweep  -  5  ysec/div 

Gain  -  750  yin/in/div 
(equivalent) 

Lower  trace  -  current 

Sweep  -  5  ysec/div 

Gain  -  20  volts/div 

Delay  -  5  ysec  from 
strain  data 


Figure  56.  Data  for  Crystal  Test  5,  Series  2 
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I.  lest  Specimen 


7  in, 


2  in. 


Test  Crystal 

.10  in.  thick  X  1/2  in.  diam. 
PZT-4  crystal 

4-i  i!M 


Impact 


Trigger  gage  (2) 


Z 


Measuring  gage  (2) 

R  =  120  ohms  F  =  2.0 

•3 


1/2-inch  diameter  rods  -  1095  Steel 


II.  Projectile  and  Velocity  Information 

Projectile  length  -  6.7  in.  Time  increment  -  13640  psec. 
Firing  pressure  -  3.5  psig.  Impact  Velocity  -  40.9  fps. 

III.  Circuit  Information 

a.  Strain  circuit. 

Number  of  active  arms  -  2 
Bridge  current  -  30  milliamp  ■ 

Equivalent  strain  -  750  pin/in/di vision 
Rcal  =  20  K  ohms/2  divisions 

b.  Current  circuit. 

F.l  =  50  ohms 

Rterm  =  50  ohms  (50- ohm  Tektronix  terminating 
risistor  was  used.) 

IV.  Scope  Data 


Upper  trace  -  Strain 

Sweep  -  5  psec/div 

Gain  -  750  pin/in/div 
(equivalent) 


Lower  trace  -  current 

Sweep  -  5  psec/div 

Gain  -  20  volts/div 

Delay  -  5  psec  from 
strain  data 


Figure  57.  Data  for  Crystal  Test  6,  Series  2 
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I.  Test  Specimen 


7  in . 


2  in , 


Test  Crystal 

.05  in.  thick  X  1/2  in.  diam. 
HST-41  crystal 
i  1  in. i  / 


Impact 


Trigger  aage  (2)- 


(  ^Measuring  gage  (2) 

R  =  1.2C  ohms  F  =  2.0 
9 

'l/2-inch  diameter  rods  -  1095  Steel 


Projectile  and  Velocity  Information 

Projectile  length  -  6.7  in.  Time  increment  -  no  data 
Firing  pressure  -  3.5  psig.  Impact  Velocity  -  approx  40  fps. 

Circuit  Information 

a.  Strain  circuit. 

Number  of  active  arms  -  2 
Bridge  current  -  30  milliamps 
Equivalent  strain  -  750  yin/in/division 
Rcal  =  20  K  ohms/2  divisions 

b.  Current  circuit. 

Rt  =50  ohms 

xj 

Rterm  =  ^0  ohms  (50-ohm  Tektronix  terminating 
resistor  was  used.) 

Scope  Data 


Upper  trace  -  Strain 

Sweep  -  5  ysec/div 

Gain  -  750  yin/in/div 
(equivalent) 


Lower  trace  -  current 

Sweep  -  5  ysec/div 

Gain  -  20  voits/div 

Delay  -  5  ysec  from 
strair  data 


Figure  58.  Data  for  Crystal  Test  3,  Series  3 
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I.  Test  Specimen 


2  in. 


2  in. 


Impact 


Trigger  gage  (2 


Test  Crystal 

.10  in.  thick  X  1/2  in.  diam. 
Quartz  crystal 

U  J  _ 


Measuring  gage  (2) 

R  =  120  ohms  F  =  2.0 
9 

"1/2-inch  diameter  rods  -  Aluminum 


II.  Projectile  and  Velocity  Information 

Projectile  length  -  6.7  in.  Time  increment  -  12964  ysec. 
Firing  pressure  -  3.5  psig.  Impact  Velocity  -  43.0  fps. 

III.  Circuit  Information 

a.  Strain  circuit. 

Number  of  active  arms  -  2 
Bridge  current  -  30  milliamps 
Equivalent  strain  -  750  yin/in/division 
Rcal  =  20  K  ohms/2  divisions 

b.  Current  circuit. 

RL  =  50  ohms 

Rterm  “  50  ohms  (50- ohm  Tektronix  terminating 
resistor  was  used.) 

IV.  Scope  Data 


Upper  trace  -  Strain 

Sweep  -  5  ysec/div 

Gain  -  750  yin/in/div 
(equivalent) 


Lower  trace  -  current 

Sweep  -  5  ysec/div 

Gain  -  .2  volts/div 

Delay  -  5  ysec  from 
strain  data 


Figure  59.  Data  for  Crystal  Test  5,  Series  3 
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I.  Test  Specimen 


2  in. 


2  in . 


Test  Crystal 

.05  in.  thick  X  1/2  in.  diam. 
Quartz  crystal 

1  1  ^  J _ 


Impact 


Trigger  gage  (2) 


Measuring  gage  (2) 

R  =  120  ohms  F  =  2.0 

9 


l/2_inch  diameter  rods  -  Aluminum 


II.  Projectile  and  Velocity  Information 

Projectile  length  -  6.7  in.  Time  increment  -  13277  ysec. 
Firing  pressure  -  3.5  psig.  Impact  Velocity  -  42.0  fps. 


III. 


IV. 


Circuit  Information 


a.  Strain  circuit. 

Number  of  active  arms  -  2 
Bridge  current  -  30  milliamps 
Equivalent  strain  -  750  yin/in/division 
Rcal  =  20  K  ohms/2  divisions 

b.  Current  circuit. 


Rl  =  50  ohms 

Rtrjrm  =  ^0  ohms  (50- ohm  Tektronix  terminating 
resistor  was  used.) 

Scope  Data 


Upper  trace  -  Strain 

Sweep  -  5  ysec/div 

Gain  -  750  yin/in/div 
(equivalent) 


Lower  trace  -  current 

Sweep  -  5  ysec/div 

Gain  -  .2  volts/div 

Delay  -  5  ysec  from 
strain  data 


Figure  60.  Data  for  Crystal  Test  6,  Series  3 
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I.  Test  Specimen 


Test  Crystal 

.100  in.  thick  X  1/2  in.  diam. 
HST-41  crystal 


II.  Projectile  and  Velocity  Information 

Projectile  length  -  6.7  in.  Time  increment  -  13133  ysec. 
Firing  pressure  -  3.5  psig.  Impact  Velocity  -  42.5  fps. 

III.  Circuit  Information 

a.  Strain  circuit. 

Number  of  active  arms  -  2 
Bridge  current  -  30  milliamps 
Equivalent  strain  -  750  yin/in/division 
Rcal  =  20  K  ohms/2  divisions 

b.  Current  circuit. 

RL  =  50  ohms 

Rterm  =  °hms  (50-ohm  Tektronix  terminating 
resistor  was  used.) 

IV.  Scope  Data 


Upper  trace  -  Strain 

Sweep  -  5  ysec/div 

Gain  -  750  yin/in/div 
(equivalent) 

Lower  trace  -  current 

Sweep  -  5  ysec/div 

Gain  -  20  volts/div 

Delay  -  5  ysec  from 
strain  data 


Figure  61.  Data  for  Crystal  Test  1,  Series  4 
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I.  Test  Specimen 


7  in. 


Test  Crystal 

.05  in.  thick  X  1/2  in.  diam. 
HST-41  crystal 
2  ir^j  J 

- 


Impact 


Trigger  gage  (2) 


Measuring  gage  (2) 

R  =  120  ohms  F  =  2.0 

9 


l/2~inch  diameter  rods  -  1095  Steel 


II.  Projectile  and  Velocity  Information 

Projectile  length  -  6.7  in.  Time  increment  -  13261  ysec. 
Firing  pressure  -  3.5  psig.  Impact  Velocity  -  42.0  fps. 


III.  Circuit  Information 

a.  Strain  circuit. 

Number  of  active  arms  -  2 
Bridge  current  -  30  milliamps 
Equivalent  strain  -  750  yin/in/division 
Rca^  =  20  K  ohms/2  divisions 

b.  Current  circuit. 


=  50  ohms 

R.  =  50  ohms  (50-ohm  Tektronix  terminatinq 
term  .  .  . 

resistor  was  used.) 

IV.  Scope  Data 


Upper  trace  -  Strain 

Sweep  -  5  ysec/div 

Gain  -  750  yin/in/div 
(equivalent) 


Lower  trace  -  current 

Sweep  -  5  ysec/div 

Gain  -  20  volts/div 

Delay  -  5  ysec  from 
strain  data 


Figure  62.  Data  for  Crystal  Test  2,  Series  4 
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I.  Test  Specimen 


2  in , 


2  in. 


Test  Crystal 

.05  in.  thick  X  1/2  in.  diam. 
Quartz  crystal 

'J-  ^  J 


Impact 


Trigger  gage  (2) 


] 


Measuring  gage  (2) 

R  =  120  ohms  F  =  2.0 
9 


1/2-inch  diameter  rods  -  Aluminum 


II.  Projectile  and  Velocity  Information 

Projectile  length  -  6.7  in.  Time  increment  -  13370  ysec. 
Firing  pressure  -  3.5  psig.  Impact  Velocity  -  41.7  fps. 

III.  Circuit  Information 

a.  Strain  circuit. 

Number  of  active  arms  -  2 
Bridge  current  -  30  milliamps 
Equivalent  strain  -  750  yin/ in/division 
Rcal  =  20  K  ohms/2  divisions 

b.  Current  circuit. 

RL  =  50  ohms 

Rterm  =  50  ohms  (50-ohm  Tektronix  terminating 
resistor  was  used.) 

IV.  Scope  Data 

Upper  trace  -  Strain 

Sweep  -  5  ysec/div 

Gain  -  750  yin/in/div 
(equivalent) 

Lower  trace  -  current 

Sweep  -  5  ysec/div 

Gain  -  .2  volts/div 

Delay  -  5  ysec  from 
strain  data 


Figure  63.  Data  for  Crystal  Test  3,  Series  4 
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I.  Test  Specimen 


Test  Crystal 

.05  in.  thick  X  x/2  in.  diam. 
HST-41  crystal 


II.  Projectile  and  Velocity  Information 

Projectile  length  -  6.7  in.  Time  increment  -  13458  psec. 
Firing  pressure  -  3.5  psig.  Impact  Velocity  -  41.4  fps. 

III.  Circuit  Information 

a.  Strain  circuit. 

Number  of  active  arms  -  2 
Bridge  current  -  30  milliamps 
Equivalent  strain  -  750  pin/in/division 
Rcal  =  20  K  ohms/2  divisions 

b.  Current  circuit. 

Rl  =  50  ohms 

Rferm  =  50  ohms  (50- ohm  Tektronix  terminating 
resistor  was  used.) 

IV.  Scope  Data 


Upper  trace  -  Strain 

Sweep  -  5  psec/div 

Gain  -  750  pin/in/div 
(equivalent) 

Lower  trace  -  current 

Sweep  -  5  psec/div 

Gain  -  50  volts/div 

Delay  -  5  psec  from 
strain  data 


Figure  64.  Data  for  Crystal  Test  4,  Series  4 
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I.  Test  Specimen 


7  in. 


2  in. 


Impact 


Trigger  gage  (2) 


Test  Crystal 

.10  in.  thick  X  1/2  in.  diam. 
HST-41  crystal 

M  iJH  J 

- . "1 


Measuring  gage  (2) 

R  =  120  ohms  F  =  2.0 

9 

'l/2“inch  diameter  rods  -  1095  Steel 


II.  Projectile  and  Velocity  Information 

Projectile  length  -  6.7  in.  Time  increment  -  13185  psec. 
Firing  pressure  -  3.5  psig.  Impact  Velocity  -  42.3  fps. 


III.  Circuit  Information 


IV. 


a.  Strain  circuit. 

Number  of  active  arms  -  2 
Bridge  current  -  30  milliamps 
Equivalent  strain  -  750  pin/in/division 
Rcal  =  20  K  ohms/2  divisions 

b.  Current  circuit. 


Rl  =  50  ohms 

Rcerm  =  50  ohms  ( 50“ ohm  Tektronix  terminating 
resistor  was  used.) 

Scope  Data 


Upper  trace  -  Strain 

Sweep  -  5  psec/div 

Gain  -  750  pin/in/div 
(equivalent) 


Lower  trace  -  current 

Sweep  -  5  psec/div 

Gain  -  50  volts/div 

Delay  -  5  psec  from 
strain  data 


Figure  65.  Data  for  Crystal  Test  5,  Series  4 
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RESULTS 


The  data  from  the  tests  were  reduced  in  two  steps.  One 
involved  using  the  strain  gage  output  from  the  measuring  gages; 
i.e.,  the  gages  at  the  axial  position  on  the  input  rod  nearest 
the  crystal  as  the  input  to  the  program  of  Appendix  B  or  C. 

The  output  from  this  program  included  tabulated  values  of 
current-time  and  energy-time.  These  values  are  called  the 
theoretical  values  in  subsequent  discussion. 

The  second  step  in  the  data  reduction  was  the  use  of  the 
measured  current-time  function  to  compute  the  electrical  energy 
dissipated  in  a  load  resistor.  This  was  done  by  programming 
for  computation  on  a  digital  computer.  Appendix  E  gives  a 
listing  of  the  program  used  for  this.  This  program  is  basic¬ 
ally  the  last  section  of  the  programs  in  Appendix  B  and  C. 

The  only  difference  is  that  the  program  of  Appendix  E  uses 
experimental  values  for  current,  while  those  of  Appendix  B 
and  C  use  theoretical  values  for  current. 

Earlier  work  had  shown  that  for  the  values  of  capacitance 
and  resistive  loading  of  the  circuit  (with  cables  and  scope 
plug-in)  encountered  with  the  circuits  used  in  the  experi¬ 
mental  work,  the  current  out  of  the  crystal  is,  for  practical 
purposes,  equal  to  that  which  passes  through  the  loading  re¬ 
sistors.  Consequently,  no  correction  was  made  for  the  effect 
of  the  measuring  instruments  on  the  current  available. 

The  data  from  the  trace  on  the  oscilloscope  picture  was 
digitized  for  use  with  the  computer  with  a  Telereader. 

Considering  first  the  results  of  the  tests  on  the  x-cut 
quartz  crystals,  it  was  anticipated  that  the  linearized  equa¬ 
tion  for  the  current  output  (Equation  5)  would  accurately  pre¬ 
dict  the  experimentally  determined  current.  This  is  because 
the  majority  of  the  assumptions  made  in  deriving  this  equation 
are  essentially  satisfied.  The  only  major  assumption  which 
is  not  satisfied  is  the  one  on  the  state  of  strain  in  the  cry¬ 
stal.  In  the  test  situation,  the  crystal  is  loaded  somewhere  be¬ 
tween  plane  stress  and  plane  strain,  and  it  was  assumed  that  the 
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crystal  was  in  plane  strain  for  data  reduction.  The  difference 
in  the  charge  coefficient  between  the  plane  stress  and  strain 
states  is  about  10%.  This  is  due  to  the  difference  in  mechani¬ 
cal  stiffness  between  the  two  states. 

The  data  of  three  runs  were  selected  for  reduction,  and 
the  resulting  plots  are  shown  in  Figures  66,  67  and  68.  These 
graphs  show  that  the  combination  of  the  simple  linear  charge- 
force  relation  for  quartz,  the  one-dimensional  elastic  wave 
propagation  equation  and  its  computer  solution  gave  excellent 
results  for  predicting  the  current  output  and  energy  release 
for  two  cases,  with  results  in  agreement  only  to  the  same 
order  of  magnitude  for  the  third  case. 

The  test  data  on  quartz  which  was  reduced  was  selected  to 
permit  comparison  of  the  effect  of  crystal  thickness  and  the 
effect  of  rod  material.  In  Figures  66  and  67,  the  only  differ¬ 
ent  pat  tmeter  in  the  tests  was  the  bar  material .  There  was 
insufficient  time  to  determine  if  this  was  the  cause  for  the 
disagreement  between  the  numerical  and  experimental  values  for 
Figure  67.  However,  the  results  of  a  duplicate  test  of  a  dif¬ 
ferent  series;  i.e.,  run  5  of  series  1,  was  in  agreement  with 
Figure  67. 

In  Figures  66  and  68,  the  only  different  test  parameter 
was  the  crystal  thickness.  A  study  of  these  graphs  indicates 
the  crystal  output  is  not  a  strong  function  of  this  parameter. 

A  study  of  the  oscilloscope  trace  for  all  the  tests  in¬ 
dicates  that  the  crystal  output  was  the  same  order  of  magnitude 
for  all  tests,  and  that  the  shape  of  the  input  pulse  may  be 
the  most  significant  parade cer  which  affects  the  output  for  a 
particular  crystal  material.  Even  minor  variations  in  the 
shape  of  the  strain/pulse  were  noted  to  have  an  effect  on  the 
output.  It  is  believed  that  the  best  way  to  investigate  the 
effect  of  such  things  as  crystal  thickness,  material  properties 
of  input  and  back-up  rod,  and  shape  of  input  pulse  would  be  to 
use  computer  simulation.  This  is  possible  due  to  the  apparent 
validity  of  the  computer  program. 
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Figure  67.  Results  of  Test  on  Quartz  Crystal  0.1”  Thick  with  Steel  Pressure  Bar 

(Series  2,  Test  3) 


Figure  68.  Results  of  Test  on  Quartz  Crystal  0.05"  Thick  with  Aluminum  Pressure 

(Series  1,  Test  3) 


The  results  of  the  tests  on  the  "new"  HST-41  crystals; 

i.e.,  from  series  4,  for  which  the  data  was  reduced  is  shown 
in  Figures  69  and  70.  The  pressure  bar  was  always  in  the 
elastic  range,  so  the  computer  program  in  Appendix  B  was  used 
in  the  data  reduction.  The  high  impact  speed  tests  which  were 
run  at  AEDC  were  reduced  utilizing  the  plastic  wave  propagation 
program  (Appendix  C) .  The  results  of  these  tests  are  shown  in 
Figures  71  and  72.  The  crystals  used  in  these  tests  were  from 
the  old  batch,  and  this  would  introduce  some  discrepancies  in 
the  results.  The  numerical  solutions  for  all  tests  were  ob¬ 
tained  by  using  the  nonlinear  charge  coefficient  determined 
from  the  static  tests. 

Figures  69  and  70  show  that  the  numerical  technique  pre¬ 
sented  for  computing  current  and  energy  from  the  crystal  gives 
reasonably  good  quantitative  and  qualitative  agreement  with 
the  experimental  results  and  can  be  considered  adequate  for  a 
first  order  approximation.  Since  the  work  on  quartz  has  shown 
that  the  numerical  technique  can  give  excellent  agreement,  the 
differences  in  the  numerical  and  experimental  results  for  the 
HST-41  is  attributed  to  deficiencies  in  the  charge-force  rela¬ 
tion,  as  determined,  and/or  as  considered  in  the  theory.  The 
most  significant  reasons  for  these  differences  are  believed 
to  be 

1.  The  effect  of  current  generated  by  radial  and 
shear  waves  is  not  negligible  for  HST-41,  and 
must  be  included  for  improved  agreement.  The 
one-dimensional  wave  propagation  theory  used 
precludes  easy  incorporation  of  this  change 
into  the  program. 

2.  The  effect  of  the  electromechanical  coupling 
has  also  been  neglected  in  this  one-dimensional 
theory . 

3.  Any  effect  of  strain  rate  on  the  force-charge 
relation  (or  current  generated)  for  the  crystal 
has  been  neglected.  The  function  determined 
under  quasi-static  loading  has  been  used  for 
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Figure  70.  Results  of  Test  on  HST-41  Crystal  0.1"  Thick  with  Steel  Pressure  Bar 

(Series  4,  Test  5) 
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the  dynamic  case,  where  loading  occurs  in 
about  10  psec. 

Two  additional  items  of  importance  are  shown  in  Figure 
69.  One  is  the  current  curve  determined  from  theory  utilizing 
a  linear  charge-force  relation  for  the  crystal.  (The  charge 
coefficient  used  was  one  determined  from  the  end  points  in 
Figure  24) .  When  this  curve  is  compared  to  the  one  deter¬ 
mined  with  the  nonlinear  relation,  it  is  noted  that  the  quan¬ 
titative  agreement  is  much  better  when  the  nonlinear  relation 
is  used. 

Also  shown  in  Figure  69  are  the  plotting  data  from  which 
the  numerical  curve  was  determined.  This  plotting  data  indi¬ 
cates  that  there  was  considerable  noise  in  the  numerical  solu¬ 
tions  for  the  current  output.  This  occurred  for  both  the 
quartz  and  the  HST-41  crystals,  and  the  scatter  in  plotting 
dots  shown  here  is  typical  of  that. for  all  computed  results  in 
Figures  66  through  70.  This  scatter  was  caused  by  errors  in 
the  data  reduction  technique  used  and  the  subsequent  amplifica¬ 
tion  of  these  errors  by  the  finite  difference  solution  of  the 
wave  propagation  problem.  All  of  the  experimental  data  was  re¬ 
corded  by  photographing  the  trace  on  the  oscilloscopes.  The 
traces  on  the  film  were  then  digitized  using  a  Telereader. 

This  required  approximating  the  vertical  center  of  the  trace 
on  the  film.  There  are  errors  in  this  approximation.  Where 
the  slope  of  the  curve  is  the  greatest,  the  errors  in  approxi¬ 
mation  were  the  greatest.  These  errors  were  further  amplified 
by  the  finite  difference  technique  used.  Unfortunately,  the 
peak  current  also  occurs  where  the  slope  of  the  strain  curve 
was  the  greatest  and  hence  the  peak  values  of  current  are  the 
least  accurately  defined,  values. 

The  noise  could  be  greatly  reduced  and  more  accurate  re¬ 
sults  obtained  by  fitting  polynomials  piecewise  to  the  strain 
data  in  the  same  manner  as  was  done  to  obtain  the  nonlinear 
charge-force  relation;  i.e.,  with  natural  cubic  interpolatory 
splines.  This  would  require  some  additional  computer  work  but 
would  yield  much  smoother  results . 
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It  should  be  noted  that  although  there  is  noise  in  the 
numerical  current  solution,  the  energy  computed  was  relatively 
smooth.  This  is  the  result  of  the  averaging  effect  of  the 
intergration  process  required  to  get  energy  from  current. 

The  results  shown  in  Figure  71  and  72  were  determined  for 
impact  speeds  of  about  8000  fps,  whereas  the  results  in  Figures 
69  and  70  were  for  impact  speeds  of  about  50  fps.  When  the 
current  and  energy  for  these  two  curves  are  compared,  it  is 
noted  that  they  are  of  the  same  order  of  magnitude.  It  should 
be  noted  that  essentially  all  of  the  energy  from  the  crystal 
is  included  in  the  curves  of  Figures  69  and  70,  but  such  is 
not  the  case  in  Figures  71  and  72.  For  these  tests,  the  energy 
given  is  that  dissipated  in  the  1.33-ohm  resistor  across  which 
the  voltage  drop  was  measured,  but  the  total  energy  generated 
by  the  crystal  is  dissipated  in  the  1.33-ohm  and  52.14-ohm  re¬ 
sistor  which  are  in  series.  When  a  correction  is  made  for  this, 

the  energy  available  from  the  crystal  for  the  high-speed  shots 

-3  -3 

was  about  7.8x10  joules,  as  compared  to  about  1.5x10  joules 

for  the  low- speed  shots.  (Values  given  are  those  from  measured 

currents . ) 

It  is  of  interest  to  note  that  the  energy  from  the  quartz 

■"8 

crystal  was  about  2x10  joules,  as  compared  to  1.5x10  ~  joules 
from  the  HST-41  crystal  for  identical  test  conditions. 
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SECTION  VII 


CONCLUSIONS  AND  RECOMMENDATIONS 


Conclusions  based  upon  the  static  tests  conducted  on  the 
crystals  used  in  this  study  are  as  follows: 

1.  The  test  set-up  developed  is  adequate  for  deter¬ 
mining  the  force-charge  relation  for  piezoelectric 
or  ferroelectric  materials,  as  indicated  by  the  re¬ 
sults  obtained  with  quartz. 

2.  The  variation  of  the  charge-force  relation  with 
previous  load  history  is  large. 

3.  There  is  an  indication  of  a  variation  in  the 
charge-force  relation  with  age  for  ferroelectric 
materials . 

4.  The  charge-force  relation  is  highly  nonlinear  for 
the  ferroelectr ' cs,  but  at  low-stress  amplitudes, 
it  may  be  approximated  as  linear.  The  stress  level 
to  which  the  linear  approximation  may  be  used  is 
dependent  on  the  material. 

5.  When  loaded  to  high  stresses,  i.e.,  50,0 tO  psi, 
the  increase  in  charge  developed  per  increment  of 
force  becomes  very  small . 

6.  The  variation  between  the  d^  values  for  the  three 
materials  tested  is  small;  generally,  smaller  than 
variations  in  values  occuring  between  different 
crystals  of  a  given  material. 

The  most  significant  conclusions  based  upon  the  results 
of  the  dynamic  tests  on  the  crystals  are: 

1.  The  theory  and  computer  programs  developed  for 
prediction  of  the  current  and  energy  from  either 
quartz  or  ferroelectric  materials,  when  loaded  in 
a  split  Hopkinson  pressure  bar,  gives  agreement 
with  experimental  results  which  is  at  least  within 
the  same  order  of  magnitude.  The  qualitative 
agreement  was  excellent  for  quartz  and  good  for 
the  ferrorelectric  materials. 
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2.  The  lack  of  good  agreement  between  computed  and 
experimental  results  for  the  ferroelectric  mate¬ 
rials  is  attributed  to  the  omission  of  phenomenon 
such  as  electromethanical  coupling,  the  use  of  a 
one-dimensional  wave  propagation  theory,  and  the 
fact  that  the  equation  for  the  short  circuit  cur¬ 
rent  from  a  crystal  may  be  only  an  approximation 
for  ferroelectric  materials. 

3.  With  the  test  geometry  used,  energy  levels  up  to 
about  2x10  ^  joules  were  obtained  from  the  ferro¬ 
electric  materials. 

4.  The  amount  of  energy  obtainable  from  the  test  ge¬ 
ometry  used  in  this  work,  for  a  given  material, 
could  be  increased  by 

a.  Using  the  optimum  size  load  resistor  in  series 
with  the  crystal. 

b.  Moving  the  crystal  very  near  the  impact  end 
of  the  Hopkinson  presure  bar.  (This  would 
tend  to  decrease  the  rise  time  of  the  input 
pulse.) 

c.  Selection  of  the  material  for  the  input  and 
back-up  rods  to  maximize  the  stress  difference. 

5.  The  computer  program  developed  for  predicting  cur¬ 
rent  and  energy  release  would  be  of  value  in  maxi¬ 
mizing  energy  release  by  ferroelectric  materials 
through  variation  and  selection  of  the  values  of 
appropriate  material  parameters  of  the  pressure 
bar  to  accomplish  this.  Further  experimental  work 
is  required  to  determine  the  optimum  load  resistor 
size.  A  two-dimensional  wave  propagation  Lheory 
would  be  required  to  properly  simulate  the  phe¬ 
nomenon  if  the  crystal  were  very  near  the  impact 
end.  (This  would  be  the  minimum  refinement  in 
theory.  Further  refinements  in  other  applicable 
theories  might  also  be  required.) 

The  recommendations  which  could  be  made  for  further  work 
in  this  area  are  numerous,  since  the  work  completed  has  been 
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primarily  in  the  exploratory  stage.  It  is  believed  the  most 
significant  areas  for  further  study  would  be: 

1.  Studies  of  the  effect  on  the  crystal  behavior 
of  such  phenomena  as  load  history,  storage  time, 
radiation,  temperature,  and/or  other  environment 
to  which  the  crystal  will  be  subjected  in  an  ap¬ 
plication. 

2.  Research  to  improve  the  relation  for  the  short 
circuit  current  generated  by  ferroelectric  mate¬ 
rial  under  shock  or  dynamic  loading. 

3.  Studies  to  maximize  the  energy  release  by  ferro¬ 
electric  materials  to  a  load  resistor,  for  the 
specific  test  geometry  used  in  this  work;  i.e., 
the  split  Hopkinson  pressure  bar. 

4.  Studies  to  optimize  the  method  of  loading  of  the 
ferroelectric  materials  in  the  field  (i.e.,  not 
under  laboratory  conditions)  to  give  the  maximum 
energy  release.  For  example,  it  is  possible  that 
crystals  in  series  in  a  pressure  bar  could  in¬ 
crease  the  energy  level,  or  it  is  possible  that 

a  shear  deformation  could  result  in  a  greater 
energy  release. 
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APPENDIX  A 


ABSTRACTS  FROM  LITERATURE  REVIEW 


NOTE:  The  numbering  system  in  Appendix  A  corresponds 
to  the  References  (p.  182).  Those  references  that  do 
not  have  an  abstract  either  cover  introductory  material 
or  have  a  title  that  is  self-explanatory. 


APPENDIX  A 


LITERATURE  REVIEW 

6.  F.  W.  Nielson,  "Ferromagnetic  and  Ferroelectric  One-Shot 
Explosive-Electric  Transducers," 

"This  memorandum  presents  the  theories  of  operation  and 
the  initial  experimental  verification  of  these  theories 
for  ferromagnetic  and  ferroelectric  one-shot  transducers 
which  convert  a  portion  of  the  energy  of  a  high  explo¬ 
sive  to  electrical  energy  with  suitable  power,  current, 
and  voltage. 

The  transducers  are  shown  to  be  extremely  small,  simple, 
and  versatile." 

9.  R.  A.  Graham,  F.  W.  Neilson,  W.  B.  Benedick,  "Piezoelec¬ 
tric  Current  from  Shock-Loaded  Quartz  -  A  Submicrosecond 
Stress  Gauge," 

"Current  from  X-cut  quartz  disks  may  be  used  to  detect 
stress-time  profiles  induced  by  shock  loading.  The  cur¬ 
rent  amplitude  and  its  time  dependence  are  functions  of 
the  dielectric,  piezoelectric,  and  mechanical  properties 
of  quartz  under  shock- loading  conditions.  The  results 
of  an  extensive  experimental  study  of  the  current  from 
shock-loaded  quartz  disks  are  reported  for  shock  stress 
up  to  50  kbar.  The  experiment  is  performed  by  impacting 
precisely  aligned  X-cut  quartz  disks  upon  each  other  at 
various  measured  velocities  and  observing  the  current  in 
one  of  the  disks  during  the  first  wave  transit.  Within 
the  low  signal  range,  the  piezoelectric  stress  constant 
e1!  is  found  to  be  0.174  C-m-2 .  The  coefficient  relat¬ 
ing  current  jump  to  stress  jump  in  one-dimensional  strain 
is  found  to  be  2. 04x10“ 8  C-cm“2-kbar-1  up  to  6  kbar  and 
2. 15x10“ 8  C-cm“2-kbar“ 1  from  9  to  18  kbar.  The  wave  velo¬ 
city  was  determined  to  be  constant  to  25  kbar.  The  ob¬ 
served  current  waveform  could  be  fully  interpreted  in 
terms  of  rate-independent  properties .  Determinations  of 
distortions  to  the  current  from  apparently  minor  devia¬ 
tions  from  one-dimensional  conditions  were  also  made." 

10.  W.  J.  Halpin,  "Current  from  a  Shock-Loaded  Short-Circuited 
Ferroelectric  Ceramic  Disk," 

"Thin  circular  disks  of  polarized  ferroelectric  ceramics 
were  equipped  with  electrodes  on  their  faces,  and  the 
electrodes  were  connected  by  a  short  circuit.  Each  disk 
was  traversed  by  a  stress  wave,  and  the  resulting  current 
pulse  through  the  short  circuit  was  measured.  The  stress 
wave  was  gem-rated  by  impact  and  made  to  propagate  along 
the  disk  axis  in  a  direction  opposite  to  that  of  polari¬ 
zation.  Current-pulse  data  are  presented  from  tests  with 
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several  tens  of  kiiobars.  Over  this  range  of  stress  both 
the  shape  of  the  current  pulse  and  its  time  integral  are 


found  to  be  strongly  stress  dependent.  It  is  concluded 
tnat  this  behavior  results  from  the  combined  effects  of 
stress  on  (1)  remanent  polarization,  (2)  polarizability, 
and  (3)  conductivity  of  the  ferroelectric  material  en¬ 
veloped  by  the  wave.  In  a  narrow  range  of  stress  where 
remanent  polarization  is  reduced  to  near  zero  and  where 
conductivity  is  negligible,  the  observed  behavior  can  be 
adequately  described  by  means  of  a  mathematical  model 
which  is  developed.  The  model  includes  the  effect  of 
stress  wavefront  tilt  which  occurs  in  impact  experiments 
to  an  extent  that  influences  significantly  the  character 
of  the  current  pulse.  Some  conclusions  are  drawn  about 
permittivity  of  shock-loaded  ferroelectrics  and  the  ex¬ 
tent  to  which  remanent  polarization  is  affected  by  the 
rapidly  applied  large  electric  fields  which  accompany 
the  process  of  current  production." 

11.  C.  W.  Beadle  and  J.  W.  Dally,  "Experimental  Methods  for 
Investigating  Strain-Wave  Propagation  and  Associated 
Charge  Release  in  Ferroelectric  Materials," 

"Study  of  strain  levels  associated  with  electric  charge 
release  in  barium-titanate-base  ferroelectric  ceramic 
which  may  be  utilized  as  single-shot,  high-energy  power 
supply,  correlation  between  charge  release  and  stress 
level  were  made  for  stress  levels  up  to  16  kbar;  in  ad¬ 
dition,  fracture  wave  propagation  was  determined,  and 
time  between  fracture  and  charge  release  was  established." 

12.  0.  M.  Stuetzer,  "Secondary  Stresses  in  a  Stress-Pulse- 
Activated  Piezoelectric  Element," 

"Previous  treatments  of  stress-pulse-activated  piezoelec¬ 
tric  and  ferroelectric  stress  gauges  and  power  generators 
neglect  the  secondary  stresses  produced  by  circuit- 
enforced  fields  in  not-directly-s tressed  parts  of  the 
transducer.  It  is  pointed  out  that  under  short-circuit 
conditions  these  secondary  stresses  generate  strong  ef¬ 
fects  in  materials  with  high  electromechanical  coupling 
factors .  To  demonstrate  this ,  a  thin  plate  activated  by 
a  plane  rectangular  low-amplitude  stress  front  is  con¬ 
sidered.  A  rigorous  theoretical  analysis  is  sketched  for 
this  configuration.  It  is  found  that  short-circuit  cur¬ 
rent  pulses  generated  by  the  stress-front  transit  are  not 
rectangular,  but  vary  exponentially  in  time.  The  ex¬ 
ponent  is  proportional  to  the  square  of  the  electro¬ 
mechanical  coupling  factor  and  depends  on  the  mechanical 
boundary  conditions.  Supporting  experiments  are  reported." 
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13.  H.  Krueger  and  D.  Berlincourt,  "Effects  of  High  Static 
Stress  on  the  Piezoelectric  Properties  of  Transducer 
Materials, " 

"Piezoelectric  ceramic  elements  in  high-power  acoustic 
transducers  are  subjected  to  high  static  as  well  as  dy¬ 
namic  stress.  This  is  particularly  true  of  well-matched 
transducers  operating  in  deep  water,  since  the  static 
stress  in  the  piezoelectric  element  may  be  several  times 
the  water  pressure.  The  present  s.udy  was  undertaken  in 
an  effort  to  determine  the  effects  of  static  compressive 
stress  on  the  piezoelectric  properties  of  two  commercial 
lead  titanate  zirconate  compositions,  PZT-4  and  PZT-5, 
and  of  two  barium- titanate  compositions,  commercial 
Ceramic  B  (a  barium  calcium  titanate) ,  and  the  composi¬ 
tion  88  wt%  barium  titanate,  12  wt%  lead  titanate 
(BaPbl2Ti) .  The  permanent  effects  of  stress  exposure, 
determined  at  zero  stress  after  exposure  to  a  given 
stress,  were  found  to  be  more  severe  with  stress  paral¬ 
lel  to  the  polar  axis  than  with  perpendicular  stres^ ,  as 
expected.  Under  maintained  stress,  however,  the  effects 
of  perpendicular  stress  are  more  severe.  PZT-4  and 
BaPbl2Ti,  generally  better  suited  for  use  as  radiating 
transducers,  show  effects  dependent  upon  exposure  time 
but  independent  of  the  number  of  stress  cycles.  Ceramic 
B  and  PZT-5  show  effects  dependent  upon  the  number  of 
stress  cycles  and  less  dependent  upon  the  total  period 
of  stress  exposure.  Of  the  compositions  tested,  PZT-4 
and  BaPbl2Ti  were  least  affected  by  high  static  stress, 
suffering  relatively  little  from  exposure  to  stress  as 
high  as'  15000  psi .  Of  these  two  compositions,  PZT-4  has 
markedly  higher  coupling  (<33=0.64  compared  to  0.365)  and 
therefore  offers  higher  transducer  bandwidth." 

14.  D.  Berlincourt  and  H.  Krueger,  "Domain  Processes  in  Lead 
Titanate  Zirconate  and  Barium  Titanate  Ceramics," 

"The  amount  of  90°  reorientation  during  poling  was  deter¬ 
mined  from  mechanical  strains  measured  during  the  poling 
process.  With  tetragonal  lead  titanate  zirconate  53%  of 
the  possible  30'’  reoiientation  occurred  during  poling, 
but  this  figure  dropped  to  44%  upon  removal  of  the  poling 
field.  With  barium  titanate  the  figures  are  only  17%  and 
12%,  respectively.  Comparison  of  the  polarization  of 
poled  polycrystalline  barium  titanate  with  that  for 
single  crystals  indicates  that  180°  reorientation  is  vir¬ 
tually  perfect.  Application  of  very  high  compressive 
stress  parallel  to  the  polar  axis  causes  90°  switching 
of  nearly  all  aligned  domains,  and,  therefore,  removes 
virtually  all  polarization.  Curves  of  released  charge 
as  function  of  mechanical  strain  are  nearly  linear,  but 
curves  of  released  charge  as  function  of  stress  are 
strongly  nonlinear.  Application  of  high  compressive 
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stress  perpendicular  to  the  polar  axis  also  causes  90° 
domai:  reorientation  and  a  reduction  in  the  total  polari¬ 
zation  of  the  ceramic.  This  domain  reorientation  may  be 
interpreted  as  a  shift  of  the  polar  axes  of  some  domains 
into  a  position  more  closely  corresponding  to  the  plane 
of  cross  expansion,  and  typically  the  total  electric 
moment  is  reduced  by  less  than  10%.  High  electric  stress 
causes  180°  as  well  as  90°  reorientation.  With  prepoled 
specimens  dc  fields  in  the  same  direction  as  the  poling 
field  cause  90°  switching,  while  reverse  dc  fields  cause 
both  90°  and  180°  reorientation,  with  the  latter  pre¬ 
dominating  . " 

15.  D.  G.  Doran,  "Shock -Wave  Compression  of  Barium  Titanate 
and  95/5  Lead  Zirconate  Titanate," 

"The  shock  compression  of  the  ceramics  BaTi03  (5%CaTiC>3) 
and  Pb  (Zrg .95Tio .05^°  (lwt%NB205)  was  measured  in  the 
ranges  5-200  kbar  and  2-140  kbar,  respectively.  Barium 
titanate  exhibits  a  two-wave  structure  above  30  kbar; 
the  first  wave  has  a  velocity  of  6.27  rom/ysec.  The  cusp 
in  the  Hugoniot  at  30  kbar  is  interpreted  as  a  dynamic 
elastic  limit.  Comparison  of  the  first  wave  velocity 
with  the  measured  longitudinal  sound  speeds  of  the  tetra¬ 
gonal  (ferroelectric)  and  cubic  (paraelectric)  phases 
(5.4  and  6.2  mm/ysec,  respectively)  suggests,  as  does 
other  evidence,  that  the  material  begins  to  transform  to 
the  cubic  phase  in  the  neighborhood  of  7  kbar.  Below 
7  kbar,  subsonic  velocities  are  observed  and  it  is  specu¬ 
lated  that  this  phenomenon  is  associated  with  domain  re¬ 
orientation. 

The  particular  lead  zirconate  composition  studied  has  a 
two-wave  structure  above  about  40  kbar;  the  position  of 
the  cusp  in  the  Hugoniot  depends  sensitively  on  initial 
density.  This  cusp  is  presumably  the  Hugoniot  elastic 
limit.  A  weak  cusp  is  also  observed  at  about  2  kbar. 

The  wave  velocity  is  essentially  sonic  below  2  kbar  and 
subsonic  above  2  kbar,  increasing  to  about  sonic  in  the 
neighborhood  of  40  kbar." 

17.  R.  K.  Linde,  "Depolarization  of  Ferroelectrics  at  High 
Strain  Rates , " 

"Poled  specimens  of  barium  titanate  (BT)  and  lead  zircon¬ 
ate  titanate  (PZT)  ceramics  have  been  successfully  re¬ 
covered  after  shock  loading  antiparallel  to  the  polari¬ 
zation  vector  (short-circuit  configuration)  in  the  3-  to 
23kbar  range.  It  has  been  shown  that  permanent  depcling 
as  a  result  of  phase  transitions  or  of  domain  switching 
can  occur  within  the  microsecond  time  scale  of  a  shock 
experiment.  Results  of  dynamic  charge  release  measure¬ 
ments  made  during  shock  transit  in  PZT  52/48  tend  to  cor¬ 
roborate  the  finding  of  other  investigations.  Techniques 
for  recovery  of  shock- loaded  specimens  were  discussed." 
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18.  C.  E.  Reynolds  and  G.  E.  Seay,  "Two-Wave  Shock  Structures 
in  the  Ferroelectric  Ceramics  Barium  Titanate  and  Lead 
Zirconate  Titanate," 

"Pb  (Zrg  ,52*  Tig  .48) -3  wit*1  1  wt%  Nb2C>5  and  pure  BaTi03, 
were  studied  under  plane  shock  waves  having  pressures 
ranging  from  2  to  175  kilobar  Which  were  produced  by  high- 
explosive  driving  systems;  2 -wave  structures  are  character¬ 
istic  of  elastic  properties  of  materials  and  are  not  due 
to  phase  transitions  from  ferroelectric  to  paraelectric. 

19.  W.  J.  Halpin,  "Resistivity  Estimates  for  Some  Shocked 
Ferroelectrics , " 

"It  has  been  observed  previously  that  the  shape  of  the 
current  pulse  delivered  to  an  external  short  circuit  by 
a  shock-loaded  ferroelectric  ceramic  specimen  undergoes 
an  abrupt  change  when  a  sufficiently  high  shock  stress 
level  is  reached.  This  change  was  explained  qualitatively 
as  being  due  to  the  onset  of  electrical  conduction  in  the 
specimen.  Here  a  quantitative  analysis  of  short-circuit 
current  waveshapes  from  some  typical  ferroelectric  mate¬ 
rials  is  made  on  the  basis  of  a  previously  developed 
model  that  has  been  extended  to  include  the  effect  of  con¬ 
duction.  Finite  resistivity  is  assumed  in  the  stressed 
region  of  a  specimen  and  evaluated  in  terms  of  the  various 
model  parameters  and  experimental  quantities,  particularly 
the  short-circuit  current  itself.  With  the  analysis  an 
estimate  is  made  that  the  resistivity  of  the  ferroelec¬ 
tric  materials  under  shock  stress  in  the  range  of  25  to 
35  kbar  is  of  the  order  of  100  ft-cm  with  accompanying 
fields  greater  than  about  20,000  V/cm." 

20.  R.  A.  Graham,  W.  H.  Halpin,  "Dielectric  Breakdown  and  Re¬ 
covery  of  X-Cut  Quartz  Under  Shock  Wave  Compression, " 

"While  a  shock  wa/e  is  traversing  a  disk  of  X-cut  quartz, 
a  piezoelectric  current  flows  in  an  external  circuit  con¬ 
nected  across  the  faces  of  the  disk.  In  this  paper  mea¬ 
surements  of  this  current  are  used  to  study  dielectric 
breakdown  and  subsequent  recovery  which  occurs  in  quartz. 
Quartz  specimen  disks  were  impacted  at  various  stress 
levels  in  such  a  way  as  to  produce  shock  waves  that  pro¬ 
pagated  along  the  X  axis  either  in  the  direction  of  or 
opposite  to  that  of  the  pressure-induced  polarization. 

In  the  latter  case,  short-circuit  current  measurements 
show  that  breakdown  occurs  at  a  threshold  stress  greater 
than  10  and  less  than  13  kbar.  Since  the  impact  experi¬ 
ment  produced  one-dimensional  electrical  and  mechanical 
conditions  in  the  specimen  disk,  it  was  possible  to 
formulate  a  mathematical  model  that  permitted  solutions 
for  internal  electrical  fields  and  resistivity  in  terms 
of  the  measured  current.  Computations  with  this  model 
show  that  the  field  in  the  stressed  portion  of  the  disk 
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at  breakdown  is  about  7.0x10  V/cm,  which  is  an  order  of 
magnitude  lower  than  the  value  observed  at  atmospheric 
pressure.  Computations  with  the  model  also  show  that 
recovery  from  breakdown  to  essentially  infinite  values 
of  resistivity  occurs  during  the  transit  time  of  the 
shock  wave  when  the  field  in  the  stressed  region  of  the 
disk  is  "quenched"  to  a  value  of  about  1.9x10 5  V/cm. 

This  critical  field  appears  to  be  the  same  for  all  shock 
stress  levels  investigated  from  13  to  35  kbar.  The  de¬ 
pendence  of  the  initiation  of  breakdown  on  the  direction 
of  wave  propagation  relative  to  the  polarization  direc¬ 
tion  indicates  that  the  shock-wave  front  furnishes  a 
source  of  free  electrons . " 

21.  R.  A.  Graham,  "Dielectric  Anomaly  in  Quartz  for  High 
Transit  Stress  and  Field," 

"Neilson  and  Benedick  have  reported  an  anomalous  piezo¬ 
electric  behavior  for  negatively-oriented  synthetic 
alpha-quartz  crystals  when  subjected  to  transient  stress 
of  about  65  kbar.  This  paper  describes  the  piezoelectric 
behavior  of  negatively-oriented  quartz  in  the  stress  re¬ 
gion  of  from  5  to  50  kbar.  The  first  indication  of  anom¬ 
alous  behavior  occurs  at  8  kbar.  Between  8  and  24  kbar 
the  negative  x-current  waveforms  show  evidence  of  partial 
electric  breakdown  in  the  quartz.  Between  25  and  34  kbar 
disruptive  breakdown  occurs.  Above  34  kbar  disruptive 
breakdown  is  followed  by  gross  conduction  with  positive 
currents  being  observed  for  stress  greater  than  about  50 
kbar.  The  fields  associated  with  the  piezoelectric  be¬ 
havior  are  lower  than  the  field  for  steady-state  electric 
breakdown  at  atmospheric  pressure.  It  is  proposed  that 
the  anomaly  is  triggered  by  stress-induced  dislocation 
motion  resulting  in  liberated  electrons  which  are  ac¬ 
celerated  into  the  stressed  region  of  the  specimen  by 
the  high  negative  electric  field." 

22.  R.  A.  Graham,  "Technique  for  Studying  Piezoelectricity 
Under  Transient  High  Stress  Conditions," 

"An  experimental  technique  is  described  which  is  being 
successfully  used  to  study  the  transient  high  stress 
piezoelectric  behavior  of  synthetic  alpha-quartz.  Short, 
flat-faced  cylinders  of  quartz  are  impacted  upon  each 
other  at  high  velocity  to  produce  the  desired  transient 
stress  in  the  range  from  5-70  kilobar.  Precision  is 
maintained  in  the  alignment  of  the  flat  impact  surfaces 
so  that  the  Hugoniot  conservation  of  momentum  and  the  re¬ 
lationship  governing  the  impact  of  flat  cylinders  may  be 
used  to  compute  the  stress  imparted  to  the  target  cylin¬ 
ders  of  quartz.  Electric  charge  release  data  are  taken 
in  conjunction  with  the  initial  passage  of  the  stress 
wave  produced  by  impact." 
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23.  E.  A.  Ripperger  and  D.  Hart,  "Stress-Charge  Release  Re¬ 
lationships  for  Ferroelectric  Ceramics," 

"The  relationship  between  charge  released,  electrical 
load,  and  stress  has  been  studied  for  two  ferroelectric 
materials,  barium  titanate  (ceramic  B)  and  lead  zirconium 
titanate  (52-48)  .  Stresses  which  cover  the  range  from 
zero  to  10  kilobars  are  applied  dynamically,  in  the  axial 
mode,  by  allowing  stress  waves  generated  by  an  impact  to 
pass  through  a  specimen  of  the  material.  The  stress  ap¬ 
plied  is  measured  by  the  split  pressure-bar  technique. 

By  properly  correlating  stress  and  charge  measurements 
which  are  recorded  as  functions  of  time,  stress-charge 
curves  for  a  given  electrical  load  are  constructed. 

These  curves  show  that  the  two  materials  investigated 
have  a  dynamic  stress-charge  relationship.  The  satura¬ 
tion  stress,  or  the  stress  at  which  all  the  available 
charge  is  released  is  between  4  and  5  kilobars  for  both 
materials.  Charge  released  is  virtually  independent  of 
electrical  load  for  pure  ohmic  loads  between  8  and  1000 
ohms.  Electrical  energy  released  increases  almost  di¬ 
rectly  with  electrical  load. 

A  PZT  in  the  "Far  Rhombohedral"  range  (92-8)  was  also 
studied  briefly.  Charge  release  from  this  material 
seems  to  be  dependent,  on  both  rate  of  stress  applica¬ 
tion,  and  stress  amplitude." 

24.  R.  W.  Rhode,  0.  E.  Jones,  "Mechanical  and  Piezoelectric 
Properties  of  Shock-Loaded  X-Cut  Quartz  at  573°K," 

"Properties  of  uniaxially  strained  X-cut  quartz  were 
studied  at  573°K  over  stress  range  of  6  to  21  kbar  by 
monitoring  current  output  from  shock- loaded  quartz  disk; 
in  addition  to  providing  fundamental  knowledge  of  be¬ 
havior  of  quartz  under  these  conditions,  results  yield 
data  necessary  for  application  of  quartz  as  submicro¬ 
second  stress  gage  at  this  temperature;  experimental 
results  on  temperature-stress  dependence  of  piezoelec¬ 
tric  and  mechanical  properties  of  uniaxially  strained 
X-cut  alpha  quartz  ar_  presented." 

25.  ...  A.  Krueger,  "Stress  Sensitivity  of  Piezoelectric 
Ceramics- 2,  3," 

"Measurements  on  two  piezoelectric  ceramics  used  for  high- 
power,  high-stress  applications  for  which  stabilizing 
heat  treatment  has  reduced  changes  in  permittivity; 
stabilizing  effect  shown  to  be  permanent;  several  stress 
cycles  to  10  to  20  kpsi  stabilize  these  ceramics;  mea¬ 
surements  of  changes  in  permittivity,  loss  and  piezo¬ 
electric  constants  for  piezoelectric  ceramic  transducer 
materials  subject  to  high  stress  perpendicular  to  polar 
axis;  difference  in  stress  sensitivity  for  "hard" 
ceramics  and  "soft"  donor-doped  ceramics." 
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26.  E.  Fatuzzo,  W.  J.  Merz,  Selected  Topics  in  Solid  State 
Physics,  Vol.  7,  Ferroelectricity, 

"A  number  of  books  and  review  articles  have  been  published 
on  the  subject,  however,  since  the  advances  made  in  the 
last  few  years,  this  book  offers  an  up  to  date  review  of 
ferroelectricity.  Little  space  is  devoted  to  the  "clas¬ 
sical  aspects"  i.e.  dependence  of  the  dielectric  proper¬ 
ties  on  pressure,  temperature,  strain,  etc.  and  the  ap¬ 
plications  of  ferroelectricity  are  only  briefly  summarized. 
These  intentional  omissions  apart,  the  book  gives  a 
thorough  physical  picture  of  the  ferroelectric  phenomena. 
The  various  topics  include:  microwave  and  infrared 
studies,  fatigue  effects,  thermodynamic  theory,  static  and 
dynamic  properties  of  domains,  lattice  dynamic  theory, 
theory  of  polarization  and  radiation  damage.  A  chapter 
is  devoted  to  recent  developments  in  ferroelectric  re¬ 
search,  these  include  special  effects  like  optical  second- 
harmonic  generation,  semi-conducting  properties  and  tem¬ 
perature  autostabilization  effects  in  ferroelectrics . 

As  shown  above  the  book  deals  with  purely  scientific  prob¬ 
lems  and  would  be  of  value  to  the  research  worker." 

27.  D.  Berlincourt,  "Piezoelectric  and  Ferroelectric  Energy 
Conversion, " 

"The  parameters  of  the  materials  involved  and  the  princi¬ 
ples  of  these  electromechanical  energy  conversion  methods 
are  discussed  in  general  terms.  Use  of  some  of  these 
materials  for  high-voltage  thermoelectric  energy  con¬ 
version  is  discussed  briefly;  power  conversion  efficiency 
appears  limited  to  the  1  percent  range." 

28.  J.  C.  Muirhead,  W.  J.  Fenrick,  "Studies  on  Shock  Wave 
Pressure-Time  Gauges  III:  The  Spontaneous  Charge  Release 
Effect, " 

"This  note  describes  further  studies  on  spontaneous  charge 
release  and  absorption  in  piezoelectric  ceramics.  Charges 
that  are  developed  by  stressing  ceramics  (or  are  impressed 
on  the  ceramic  from  a  power  source)  tend  to  change  the 
polarities  of  individual  crystals  within  the  ceramics, 
which  thereby  store  charges .  Application  of  pressure 
pulses  assists  in  returning  these  modified  crystals  to 
their  normal  state.  This  in  turn  releases  the  stored 
electrical  charges,  which  show  up  on  pressure-time  re¬ 
cords  as  spontaneous  charge  release  or  absorption." 

29.  W.  Ohlwiler,  "High  Temperature  Ferroelectric  Materials," 

"This  report  reviews  published  data  on  high-temperature 
ferroelectric  materials.  Applicable  terms  are  defined 
and  some  characteristics  of  this  type  of  material  are 
given.  High  temperature  ferroelectric  materials  studied 
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were  lead  metaniobate,  lithium  tantalate,  potassium 
niobate,  lead  titanate,  and  lead  metantantalate .  Other 
ferroelectric  materials  with  less  supporting  data  are 
discussed  also." 

30.  H.  Jaffe,  D.  A.  Berlincourt,  "Piezoelectric  Transducer 
Materials , " 

"Properties  of  piezoelectric  crystals  and  ceramics  are 
reviewed  as  they  affect  use  of  such  materials  in  electro 
acoustic  transducers.  Extensive  data  on  lead  titanate- 
zirconate  ceramics  are  presented  to  help  in  selection  of 
the  appropriate  material  for  a  variety  of  applications 
in  the  lower  and  medium  frequency  range,  up  to  several 
megacycles  per  second.  Quartz  and  several  recently  dis¬ 
covered  piezoelectric  crystals  will  share  the  higher 
megacycle  range . " 

31.  Bruce  D.  Wedlock,  "Properties  of  Piezoelectric  Materials 

"The  purpose  of  this  monograph  is  to  provide  a  source  of 
data  on  a  wide  variety  of  piezoelectric  materials  which 
will  be  useful  to  those  conducting  research  in  this  area 
as  well  as  to  the  engineer  designing  ultrasonic  trans¬ 
ducers,  filters,  and  other  piezoelectric  devices.  The 
numerical  data  are  given  in  the  rationalized  MKS  system 
in  accordance  with  the  IRE  Standards  on  Piezoelectricity 
as  are  the  symbols  employed.  Table  I  defines  the  sym¬ 
bols  employed  and  gives  their  MKS  unit.  Since  a  large 
amount  of  literature  still  uses  the  CGS  system.  Table  II 
presents  the  factors  to  convert  CGS  units  to  MKS  units . 
On  the  data  pages,  the  digits  in  the  data  columns  are  to 
be  multiplied  by  the  power  of  ten  in  the  column  immedi¬ 
ately  to  the  right  of  the  symbol  to  obtain  the  value  in 
the  MKS  system.  The  compounds  are  arranged  alphabetic- 
ally  by  chemical  name.  The  pure  compounds  are  followed 
by  the  binary  alloys,  followed  by  ternary  compounds,  etc 
The  multiple  component  compounds  are  arranged  alphabetic 
ally  by  main  constituent,  secondary  constituent,  etc., 
and  by  decreasing  amount  of  main  constituent." 

32.  I.  Lefkowitz,  K.  Kramer,  and  P.  Kroeger,  "Ferroelectrics 
Their  Electrical  Behavior  During,  and  Subsequent  to. 
Ionizing  Radiations," 

"A  study  was  made  on  certain  ferroelectric  materials 
which  are  mong  those  used  extensively  in  military  appli¬ 
cations,  both  in  a  polarized  condition  in  such  items  as 
impact  fuzes  and  in  an  unpoiarized  condition  as  capa¬ 
citors.  The  transient  voltages  across  a  load  produced 
by  charge  generated  on  the  surface  of  various  types  of 
ferroelectric  specimens  during  irradiation  have  been 
measured.  These  outputs  were  found  to  vary  widly  in 
amplitude  and  polarity  even  from  specimens  ostensibly 
alike  (i.e.,  from  the  same  lot  of  a  given  manufacturer). 
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and  even  from  a  single  specimen  pulsed  repeatedly  with 
gamma  and  neutron  radiation.  Pulses  were  found  to  vary 
from  background  levels  to  a  high  of  300  volts  on  a  large 
specimen  with  a  load  of  10  ?  ohms  and  radiation  of  5.4  x 
10 fast  (Pu)  neutrons  and  associated  gammas.  The 
maximum  voltage  possible  is  unknown,  as  are  the  factors 
causing  the  variation.  Both  polarized  and  unpolarized 
specimens  showed  outputs  well  above  background  transients . 
Some  specimens  were  identical  to  those  used  in  impact 
fuze  applications  and  were  in  simulated  housings  which 
reproduce  the  mechanical  environment  of  the  device. 

Some  fuze  circuits  in  present  use  are  discussed  with 
reference  to  their  use  of  ferroelectric  materials  and 
a  calculation  of  the  energy  transferred  by  a  voltage 
pulse  to  a  load  is  presented." 
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APPENDIX  B 

LIST  OF  COMPUTER  PROGRAM 

FOR 

ELASTIC  WAVE  PROPAGATION  IN  A  COMPOSITE  BAR 
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General : 

The  program  "CRYSTAL"  instructs  the  computer  to  perform 

the  following  operations  in  the  order  listed: 

(1)  Read  in  the  structural  data  including  geometry  and 
material  properties  and  circuitry  data. 

(2)  Read  in  the  time  increment. 

(3)  Read  in  the  strain  input  data  from  the  experimental 
scope  recordings. 

(4)  Calcualte  the  initial  values  of  the  field  variables. 

(5)  Calculate  the  one  dimensional  wave  velocity, 

C  -  /E/p  for  the  three  media. 

(6)  Calculate  the  values  of  the  field  variables,  a,  e,  and 
v  for  each  mesh  point,  up  to  the  leading  elastic  wave 
front. 

(7)  Calculate  the  stress  difference  across  the  crystal  faces. 

{8)  Calculate  the  short  circuit  current. 

(9)  Calculate  the  energy  produced  by  the  short  circuit 

current. 

(10)  Print  on  the  field  variables  at  desired  mesh  points. 

Also  print  on  the  stress  difference,  short  circuit 
current  and  energy. 

(11)  Repeat  steps  (6)  through  (10)  for  each  time  increment. 

Data : 

(1)  El,  E2,  E3,  RHOl ,  RII02 ,  RH03  (FORMAT  100) 

(2)  IMAX,  LI,  L2 ,  X2,  DIEC,  DIA  (FORMAT  101) 

(3)  Ml,  M2,  M3,  R,  DT,  FACT  (To  be  given  in  the  body  of 
the  program) . 

(4)  XINPUT  (To  be  given  in  the  body  of  the  program). 

El,  E2,  E3  =  Youngs  modulus  for  medium  1,2,3  resp. 

RHOl,  RH02 ,  RH03  =  Mass  density  for  medium  1,2,3  resp. 

IMAX  =  Total  number  of  time  incremer  '■:» 

LI  =  Spatial  increment  number  to  the  front 
face  of  the  crystal 
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L2  =  Spatial  increment  number  to  the  rear 
face  of  the  crystal 
X2  =  Thickness  of  the  crystal 
DIEC  =  Dielectric  constant  for  the  crystal 
DIA  =  Diameter  of  the  bar 
Ml  =  1 

M2  =  Spatial  increment  number  to  the  middle 
of  the  crystal 

M3  =  Spatial  increment  number  of  teh  output 
gage  in  the  backup  bar 
P  =  Short  circuit  resistance 
DT  =  Time  increment 
FACT  =  Scope  vertical  gain 
XINPUT  =  Vertical  deflection  of  scope  for  input 
gage 

SIG,  STS  =  Stress  at  tiem  =  t  and  (t-DT)  resp. 

EPS,  STR  =  Strain  at  time  =  t  and  (t-DT)  resp. 

VEL,  V  =  Particle  velocity  at  time  =  t  and  (t-DT) 
resp. 

Cl  =  Short  circuit  current 
ENERGY  =  Energy  produced  by  the  short  circuit 
current 
T  =  Time 

Di  =  1.0/pi*Ci;  (i=l,2,3) 

Ci  =  /Ei/pi  ;  (i-1,2,3) 

Running  time: 

If  the  calcualtion  is  done  for  t=o  to  t=30ysec 
DT=0.25ysec,  the  running  time  is  approximately  15  seconds. 
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Figure  73.  Computer  Program  Flow  Chart  for  Elastic  Wave  Propagati 


JFn »60  »  30* 55000. VALATHUR, 5791, UNM- 
RUN ( S ) 

CRYSTAL. 

i 

PROGRAM  CRYSTAL ( INPUT* OUTPUT) 
r*****#t  ***** *******  ***4**************  ********* 

C*ELASTIC  WAVE  PROPAGATION  IN  A  COMPOSITE  BAR  * 
r  ***  ******************************************* 

DIMENSION  SIG(501 ) »EPS(5C1) »VFL( 501) , STS (501) 
1*XINPUT(5C1) *FMfRGY(501 ) »CI (501) ,STR(501) ,V( 501) 

c  strain  input  problem 

C  LENGTH  OFMEDIUMl*Xl* 

C  LENGTH  0FMEDIUM2*X2* 

C  DX1  IS  THC  X  INCREMFNT  IN  MEDIUM1 
C  DX2  IS  THC  X  INCREMENT  IN  MEDIUM? 

C  DX3  IS  Ihc  X  INCREMFNT  IN  V€DIU«3 
C  E  IS  y''UNGS  MODULUS 
r  RHO  is  MASS  DFNSITY 

C  R  IS  RESISTANCE  OF  THE  STRAIN  GAGES  IN  OH^S 

c  dt  is  the  time  increment 

C  L1=1+X1/DX1 

c  L2=li+;;2/dx2 
C  L3=l2+X3/DX1 

C  IMAX  IS  THE  TOTAL  NUM8ER  OF  INCREMENTS  OF  TI^IE 
C  IMAX=TOTAL  TIME  /  DT 

C  N0TE**CH00SF  DT  SUCH  THAT  (L2-L1JIS  AM  INTEGER** 

C  DATA  TO  BF  SUPPLIED  FOR  THE  PROBLEM 
'-******************************************** 

C  *NOT e  on  UNITS***INCHES,POUNDS»SECONDS** 

r****.»*4******ii  ********4********************* 

r  UNITS  OP  DATA,  F=LR/IN**2,  RH0=LR-SEC**2/IN**4 
C  UNITS  OF  THE  DATA.  X2=IN.,  DIEC=AMP-SCC/L8. ,  DIA=IN. 
450  READ  100,E1*E2,-3,RH01,Rh02»RH03 

100  FORMAT ( 6E10.3 ) 

READ  101»IMAX,L1,L2,X2,DIEC»DIA 

101  FORMAT (3I5»F10«3*E10.3»P10.3) 

C  OUTPUT  STRAIN  GAG1*  LOCATIONS 

0  Ml  IS  INPUT  POINT 

M 1  =  1 

C  M2  IS  THE  MIDPOINT  OF  MEDIUM  2 
C  V3  IS  THP  SECOND  GAGE  LOCATION 
AREA=^.1417*(DIA**2. )/a. 

R=53.47 
DT=0. 1 BE-06 

C  DETAILS  OP  STRAIN  INPUT  FROM  EXPERIMENTS 
C  VALUES  OF  FIELD  VARIABLES  AT  TIME=0  SPECIFIED  BY 
C  SUBSCRIPT(l) 

IMAXP1=IMAX+1 
T  =  C  .0 

C  XINPUT  IS  THE  VERTICIL  DEFLECTION  OF  SCOPES  FOR 
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C  STRAIN  GAGES 

DO  80°  J= 1  ♦  56 
T=T+DT 

800  XINPUK  J)=2.6E+05*T 
DO  801  J=57» 112 
T=T+DT 

801  XINPUK J)=2.6+6.0E+04*(T-10.0E-06) 

DO  802  J= 1 13  » I MAXP1 

T=T+OT 

802  XINPUK  J)=3. 2 

C  SCALE  FACTOR  POR  INPUT  IS*FACT* 

FACT=2P24.CE-06/4.0 
C  END  OF  INPUT  STRAIN 
475  CONTINUE 

c  the  speed  of  propagation  is  in  units  of  in. /sec. 

C1=(SQRT(E1/RH01) ) 

C2* (SQRT ( E2/RH02 ) ) 

C3= ( SORT ( E3/RH03 ) ) 

D1=1.0/(RH01*C1) 

02  =  1.0/ ( RH02*C 2  ) 

D3=1.0/(RH03*C3) 

STR ( 1 )  =  XI NPUT ( 1 )  *FACT 
STS ( 1 )  =E1*STR ( 1 ) 

V(1!=D1*STS(1) 

ENERGY ( 1 ) =0 .0 

C  FIELD  VARIABLES  FOR  TIME  T=DT 
EPS(1)=FACT*XINPUT(2) 

SI G( 1 ) =STS ( 1 J+El* ( EPS ( 1 )-STR ( 1 ) ) 
VEL(1)=V(1)-D1*(SIG(1)-STS(1)  ) 

EPS ( 2 ) =STR ( 1 ) 

SIG(2)=STS(1) 

ENERGY(2)=ENERGY(l) 

VEL ( 2 ) =V ( 1 ) 

PRINT  299 

C  THE  NUMBER  AND  DATE  OF  THE  DATA  IS  PRINTED  5Y  FORMAT  169 
501  PRINT  169 

169  FORMAT ( 15X  »23H  RUN  NUMBER  3  7/10/1968,/) 

PRINT  199  ,  IMAX 
PRINT  20  ,X2 

20  FORMAT ( 3X  » 11H  Th! C<NESS= , F 10 .4  ) 

PRINT  308  ,  D I  EC 

3C3  FORMAT ( 3X  »21H  DIELECTRIC  CONSTAN T= » F10.3 ♦ 10HAMP-SEC/LB » 
1//) 

PRINT  300 
PRINT  310 

PRINT  301  ,  El »E2,E3,RH01,RH02,RH03*L1,L2 
pR I  NT  311 
PRINT  312 
PRINT  302 
T  =  DT 
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PRINT  3C3»  T,EPS(Ml)tSIG(Ll) ,EPSBM1.PPS(L1) *SIG(L2) » 
1EPSBM3*EPS(L2) *DFSTR»CI (1) * ENERGY ( 1 ) 

MIX=3 

DC  2CC  IX=3»IVAxP] 

T=T+DT 
IX1=IX-1 
DO  201  IJ=l*IXl 
STS( I J)=SIG( IJ) 

STR( IJ)=FPS( IJ) 

201  V( I J ) =VEL ( IJ) 
f  AREAS  THE  MED I A  ARE  EQUAL 

r  AT  THE  INTERFACE  FORCES  AND  VELOCITIES  ARE  SAME 
S  T  R ( IX ) =STR ( 1X1 ) 

STS(IX)=STS( 1X1) 

V(  IX)=V( 1X1) 

EPM1)=FACT*XINPUT(  IX)  . 

51 G( 1 ) =STS( 1 )+El* ( EPS( 1 )-5TR ( 1 ) ) 

VEL( 1 )=V( 1)-D1*(SIG( 1)-STS( 1 ) ) 

SIG( IX) =STS(  1X1 ) 

VEL ( I X ) =V ( I XI ) 

IF ( IX-L1)212»21?»214 

213  EPS( IX)=SIG( IX)/E2 
GO  TO  215 

214  IF( IX-L2)212*212*215 

215  FPS(  IX )=SIG( IX) /E3 
GO  TO  216 

212  FPS( IX)=STR(IX1) 

216  DO  20?  I J=2» 1X1 
I  JL=I J-l 

I  JR= I J+l 

IF ( IJ-L1)203.204,205 

204  DL=1.0/D1 
QR= 1 «  0/D2 
GO  TO  209 

205  I F ( IJ-L2)206»207»203 

207  DL=1.Q/D2 
f)R= 1 • 0/D3 

209  VEL ( I J )  =  ( STS( IJR)-STS( IJL)+DL*V(  I JL ) +DR*V ( I  JR ) ) / ( DL+DR ) 
S I G( I J ) =  ST S ( IJL)+DL*(VFL( IJ)-V(I  JL) ) 

EPS ( I J ) =STR ( IJ) ■ ( SIG ( I J)-STS( IJ) )/E2 
GO  TO  202 

206  PX=D2 
EX  =  P2 

GO  TO  210 

208  DX=D3 
FX=r3 

GO  TO  210 
203  DX=0 1 
EX  =  r  1 

210  V5L ( IJ)  =  (V( I  JR ) +V ( I JL )+CX* ( STS  II  JR )-Sf S( I JL ) ) )/2.0 
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S I G( I J)=STS< I JL  J+UVEL ( I J ) -V ( I JL ) )/DX 
FPS( I J)=STR( IJ)+(SIG( IJ)-STS( I J) ) /EX 
CONTINUE 

EPS9M1=SIG(L1)/F1 


FPSBM3=SIG(L2) /E3 

C  DFSTR  IS  THE  STRESS  DIFFERENCE  ACROSS  THE  CRYSTAL 
DFSTR=SIG(L1)-SIG(L2) 

C  CONVERSION  FACTOR  FROM  WATT-SFC  TO  JOULES  IS  'CON' 


C0N=1 «0 

C  Cl (IX)  IS  THE  SHORT  CIRCUIT  CURRENT  AS  CALCULATED  BY 
C  STRESS  EQUATION  AND  THE  UNITS  ARE  AMPERES. 

Cl (  I X ) =DFSTR*AREA*C2*DIEC/X2 

C  ENFRGY  IS  THE  POWER  PRODUCED  BY  THE  SHORT  CIRCUIT  CURRENT 
C  ***  ENFRGY=DX/3*R*(CI (N)**2+4*CI (N+1)**2+CI (N+2)**2) 

IF( IX-MIX)600,600,700 
600  CO'-iINUE 

ENERGY ( IX)=DT/3.*R*CON*(CI ( IX-2) **2+4. *CI ( IX-1)**2+ 

1CI  ( IX)**2 )+ENERGY( IX-2) 

211  PRINT  303 ♦  T,EPS(M1),SIG(L1) .EPSBMl.EPS(Ll) »SIG(L2) ♦ 
1EPSBM3 .EPS ( L2 ) .DFSTR.CI ( IX) .ENERGY ( IX) 

GO  TO  200 
700  CONTINUE 

PRINT  303 ♦  T  » EPS (Ml ) »SIG( LI ) »EPSBM1 »EPS(L1) »SIG(L2) » 
1EPSRM^,EPS(L2) .DFSTR.CI (IX) 

MI X=M I X+2 
200  CONTINUE 

299  FORMAT( 1H1.10X.44H  ELASTIC  WAVE  PROPAGATION  IN  A 
1COMDOSITE  BAR/) 

199  cORMAT ( 3X  »6H  IMAX=.I5./) 

300  FORMAT ( 3X , 3H  E1.7X.3H  E2.7X.3H  E3.6X.5H  RH01.5X. 

15H  RH02.5X.5H  RH03.5X.3H  L1.2X.3H  L2 ) 

301  FORMAT (6E10. 3,215//) 

302  FORMAT ( // »2X,5H  TtME.6X.6H  EPSM1.5X.6H  SIGL1.4X, 

1 7H  ^dSRM1,5X,6H  ECSL1.5X.6H  SIGL2.4X.7H  EPSBM3.5X, 

26H  FPSL2.5X.6H  DFSTR.2X.8H  CURRENT ,4X ,7H  ENERGY) 


303 

310 

311 


312 


1000 


FORMAT (E10.3.10F11.4) 

FORMAT  (  1X.27H....3....  (Li3/IN**2) . .4X.29H 

1(LB-SFC**2/IN**4) . . 


FORMAT ( 14H****LCGEND****» / .4X ♦ 35HEPSM1  =  STRAI N  AT  FI 
1RST  GAGE  LOCATION, /,4X,40HSIGL1=STRESS  IN  MEDI 
2UM  1  AT  CRYSTAL  FACE»/,4X.41HEPSBM1=STRAI N  IN  MED 
3  I  UN!  1  AT  CRYSTAL  FACE  »  /  »4X  »40HEPSL1  =  STRAI  N  IN  MED 
4IUM2  AT  CRYSTAL  FACE »/ ,4X ,33HSIGL2=STRESS  AT  CRY 
5STAL  REAR  FACE  »/»4X»4lHEPS3M3=STRAIN  IN  MEDIUM  3  A 
6T  CRYSTAL  FACE./.4X, 36HEPSL2=STR AIN  IN  CRYSTAL  A 
71  P.PAR  FACE) 


FORMAT ( 4X ,38HDFSTR= STRESS  DIFFERENCE  ACROSS 
1AL,/,4X»40HCURRENT=CIRCUIT  current  in  units 
2PS »/.4X.40HENERGY=CRYSTAL  ENERGY  INUNITS  OF 
CONTINUE 
STOP 


CRYST 
OF  AM 
JOULES,/) 


END 
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APPENDIX  C 


LISTING  OF  COMPUTER  PROGRAM 

FOR 


PLASTIC  WAVE  PROPAGATION  IN  A  COMPOSITE  BAR 


General : 

The  program  "WAVES"  instructs  the  computer  to  perform 

the  following  operations  in  the  order  listed: 

(1)  Read  in  the  s  ructural  data  including  geometry  and 
material  properties  and  circuitry  data 

(2)  Read  in  the  time  increment 

(3)  Read  in  the  strain  input  data  from  the  experimental 
scope  recordings 

(4)  Calculate  the  initial  values  of  the  field  variables 

(5)  Calculate  the  one  dimensional  wave  velocity,  C  =  /E/p 
for  the  three  midia 

(6)  Knowing  the  values  of  the  field  variables  at  the  pre¬ 
vious  instant  of  time  ie  t-t-Dt,  the  slopes  of  the 
characteristics  starting  from  the  previous  time  base 
and  passing  through  the  mesh  points  considered  are 
calculated.  This  is  done  by  calling  the  subroutine 
"INTERP".  This  subroutine  calculates  the  field  var¬ 
iables  at  the  starting  point  of  the  characteristic  by 
linear  polation. 

(7)  Calculate  the  values  of  the  field  variables,  a,  e,  and  v 

(8)  Repeat  steps  (6)  and  (7)  for  all  the  mesh  points,  up 
to  the  leading  elastic  wave  front. 

(9)  Calculate  the  stress  difference  across  the  crystal  faces 

(10)  Calculate  the  short  circuit  current 

(11)  Calculate  the  energy  produced  by  the  short  circuit 
current 

(12)  Print  out  the  field  variables  at  desired  mesh  points. 
Also  print  out  the  stress  difference,  short  circuit 
current  and  energy 

(13)  Repeat  steps  (6)  through  (12)  for  each  time  increment 

Data: 

(1)  El,  E2,  E3,  RHOl ,  RH02,  RH03  (FORMAT  100) 

(2)  IMAX,  LI,  L2,  X2,  DIEC,  DIA  (FORMAT  101) 

(3)  XK,  XN  (FORMAT  99) 
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(4)  YIELD,  R,  DT,  FACT,  Ml  (To  be  given  in  the  body 
of  the  program) 

(5)  XINPUT  (To  be  given  in  the  body  of  the  program) 


List  of  symbols: 

El,  E2,  E3 
RH01 ,  RH02 ,  RH03 
IMAX 
LI 

L2 

X2 
DIEC 
DIA 
XK,  XN 
YIELD 
R 
DT 
FACT 
XINPUT 

Ml 

SIG,  STS 
EPS,  STR 
VEL,  V 

F,  S 

Cl 

ENERGY 

T 

Dj 


Youngs  modulus  for  medium  1,2,3  resp. 
Mass  density  for  medium  1,2,3  resp. 
Total  number  of  time  increments 
Spatial  increment  number  to  the  front 
face  of  the  crystal 

Spatial  increment  number  to  the  rear 
face  of  the  crystal 
Thickness  of  the  crystal 
Dielectric  constant  for  the  crystal 
Diameter  of  the  bar 

Material  parameters  for  the  input  bar 
Yield  stress  for  input  bar 
Short  circuit  resistance 
Time  increment 
Scope  vertical  gain 

Vertical  deflection  of  scope  for  input 
gage 
1 

Stress  at  time=t  and  (t-DT)  resp. 
Strain  at  time=t  and  (t-DT)  resp. 
Particle  velocity  at  time=t  and  (t-DT) 
resp. 

Function  variable  at  time-c  and  (t-DT) 
resp. 

Short  circuit  current 

Energy  produced  by  the  short  circuit 

current 

Time 

1 . 0/ p j  *Cj ;  ( j=l ,2,3) 
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Cj  =  /Ej/pj  ;  (j=l,2,3) 

XNl  =  XN-1 
XN12  =  XN1/2.0 

SGI,  SGV,  VI,  VLV  =  Variable  quantities 
Running  time: 

If  the  calculation  is  done  for  t=o  to  t=30ysec  with 
DT  =  O.lysec. ,  the  running  time  is  approximately  100  seconds. 
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Read  _ m  Value  of  T Calculate  the  t  ,  Se t  t  Initial 

Data  Strain  Input  I  *1  Wave  Velocities  Ix=Ix+l  f  Variabl 


Figure  74.  Computer  Program  Flow  Chart  for  Plastic  Wave  Propagation  in  a  Composite  Bar 


JFH  ,6  , 50C» 55000. VALATHUR* 579 1»UNM* 3743. 

RUN ( S » » , ♦ ,*25000) 

WAVES. 

i 

program  waves ( input  *  OUTPUT ) 

COMMON  SGI ,SGV,VI ,VLV»Rh01 ,C1 ,XK ,XN,DT ,DX1,SGL,VL* 
1X£1,CF,RCF,XKN,XN1 
COMMON  YIELD 

DIMENSION  SIG(5C1),EPS(501)»VEL( 501 5 ,STS< 501 ) *STR ( 501 ) 
l»S(50l) *F(5Cl) ,XINPUT(501) *CI ( 501), ENERGY ( 501 ) ,V( 501 ) 
r **************  **********  ************* ********* 

C*PLASTIC  WAVE  PROPAGATION  IN  A  COMPOSITE  BAR  * 

q**** ****************** ft** ********* ************ 

C  MEDIUM  2  IS  IN  ELASTIC  STATE  AND  MEDIA  1AND3  ARE 
C  IN  PLASTIC  STATE 

Q******* *************** ft*** 

c*  strain  input  problem* 

£****•#*#  ************  ******* 

C  LENGTH  OF  MEDI UM2*X2* 

C  DX1  IS  THE  X  INCREMENT  IN  MEDIUM1 

C  DX2  IS  THE  X  INCREMENT  IN  MEDIUM2 

C  DX3  IS  THE  X  INCREMENT  IN  MEDIUM3 

C  STRESS  IS  EXPRESSED  IN  KSI  UNITS 

C  E  IS  YOUNGS  MODULUS 
C  RHO  IS  MASS  DENSITY 

C  XK  AND  XN... .MATERIAL  PARAMETERS ( RATE  INDEPENDENT 

C  THEORY) 

C  *R*  IS  THE  RESISTANCE  OF  THE  CIRCUIT  IN  OHMS 

C  DT  IS  THE  TIME  INCREMENT 
C  L1=1+X1/DX1 
C  L2=L1+X2/DX2 

C  IMAX  IS  THE  TOTAL  NUMBER  OF  INCREMENTS  OF  TIME 
C  ImaXsTOTAL  TIME  /  DT 

C  N0TE«*CH00SE  DT  SUCH  THAT  (L2-L1) IS  AN  INTEGER** 

q* ****•>: )  >****ft»******-»*^** ************ ********* 

C*****NOTE  ON  UNITS. . .KIPS, INCHES » SECONDS****** 

C*** ***********  ********  ************************ 

C  UNITS  OF  DATA...E=KIPS/IN**2»RH0=KIPS-SEC**2/IN**4 

C  UNITS  OF  DATA...X2=INCHES,DIEC=AMP-SEC/KIP,DIA=INCH 

C  DATA  TO  BE  SUPPLIED  FOR  THE  PROBLEM 

READ  IOC  »E1»E2*E3»RH01»RH02»RH03 

100  FORMAT  <  6E10.3 ) 

READ  101»IMAX,L1,L2,X2,DIEC»DIA 

101  FORMAT (3I5»F10.3»E10.3,F10.3) 

READ  99»XK,XN 

99  FORMAT! 2E 12. 4) 

C  YIELD  STRESS  FOR  MATERIAL  1  IS  EGUAL  T030  KSI  (YIELD) 

C 
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n  n  rt  n  n 


n— xt  jul  uj 


Y I ELD=3C • 0 

DETAILS  OF  STRAIN  INPUT  FROM  EXPERIMENTS 
SCALE  FACTOR  FOR  INPUT  IS  #FACT* 

0T=0.1SE-06 
R  =  5  3  •  4  7 
IYAXP1=IMAX+1 
T  =  0 .0 

XINPUT  IS  THE  VERTICAL  DEFLECTION  GF  SCOPES  FOR 

strain  gages 

READ  803. ( XINPUT ( J ) ,J=1 .IMAXP1 ) 

800  FORMAT (  1 5F5 • 2 ) 

C  SCALE  FACTOR  FOR  INPUT  IS  FACT 

FACT=2.985E-03 
C  end  of  READING  TH^  INPUT 

AREA=22.0/7.0*DIA**2/4,C 
C1=SQRT ( El/RHOl ) 

C2=SQRT ( E2/RH02 ) 

C3  =SORT ( E3/RH03 ) 

Dl = 1 • 0/ ( RH01*C 1 ) 

D2=1.0/(RH02*C2) 

D3  =  l • r/ ( RH03*C  3 ) 

01=1. ~/D2 

DX1=C1*0T 

DX2=C2*DT 

DX3=C4*DT 

XE1=1./E1 

XE2=1.0/E2 

XKN=XK*XN 

XN1=XN~1 • 3 

XN12=XNl/2.0 

C  SPECIFY  THE  VALUCS  OF  THE  FIELD  VARIABLES  FOR  ZERO  TIME 
C  BY  SUBSCRIPT ( 1 ) 

STR(1)=XI.\PUT(1)  *fact 
S  T  S ( 1 ) =E1*STR ( 1 ) 

V( 1 )=D1*STS( 1 J 
S(1)-1.0/E1 
ENERGY ( 1 ) =0  #C 

C  FIELD  VAR  I AaL  F  S  FOR  TIME  T  =  0T 
EPStl )=FACT*XI!\PUT(  2  ) 

SIGC 1)=STS(1)+E1*(EPS( 1 )  —  ST  R (  1)) 

STRS2=SIG! 1)*  SIG( 1) 

F( 1 ) =XE1+XKN*STRS2**XN12 

VEL ( 1 ) =V ( 1 ) — Dl * ( S I G C l)-STS(l)  ) 

FPS(2)=STR(1) 

S I G ( 2 ) =ST  S ( 1 ) 

VEL ( 2 ) =V ( 1 ) 

F  (  2  )  —  S  (  1 ) 

FNERGY ( 2 ) =0.0 
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r>  r» 


PRINT  299 
PRINT  169 

169  FORMAT ( 15X»24H  RUN  NUMBER  ARO  8/7/1963./) 
PRINT  808  »  DIEC 

808  FORMAT ( 3X  »21H  DIELECTRIC  CONSTANT=,E10.3. 
112H  AMP-SEC/KIP,//) 

PRINT  300 

PRINT  301  .  E1.E2.E3.RHC1.RH02.RH03.L1.L2 
PRINT  311 
PRINT  312 
PRINT  302 
T=DT 
MI  X=3 

FPM1=1.0/E1 
FPM3=1.0/E3 

output  strain  gage  locations 

Ml  IS  INPUT  POINT 
Ml  =  l 

DO  200  IX  =  3» IMAXPl 
T=T+DT 
IX1=IX-1 
DO  201  IJ-=1.IX1 
STS(IJ)=SIG( IJ) 

STR(IJ)=EPS( I J  ) 

S( I J ) =F ( I J ) 

201  VC  I J)=VEL( I J) 

C  AREAS  OF  THE  MEDIA  ATE  EQUAL 

C  AT  T8E  INTERFACE  FOR3ES  AND  VELOCITIES  ARE  SAME 
STR(IX)=STR{ 1X1) 

STS ( I X ) =STS ( I XI ) 

V( I X ) =V  C I XI ) 

S ( I  X ) =S ( I XI ) 

EPS ( 1 ) =FACT*XI NPUT (IX) 

SIGC 1 ) =STS  < 1 )  +  (EPS(l)-STR(l) ) /SC  1) 

STRS2=SIG( 1 )*  SIGC 1 ) 

F ( 1 )=XE1+XKN*STRS2**XN12 
SGI =STS ( 1 ) 

VI=V(1) 

SGV=STS ( 2 ) -STS ( 1 ) 

VLV=V(2)-V(1) 

CALL  INTERP 

VEL(1)=VL-(SIG(1)-SGL)/RCF 
SIGC IX)»STS( I X 1 ) 

VEL (IX) = V ( I X 1 ) 

I F  (  IX-LD212.213.214 

213  EPS ( I  X ) =S I G ( I X ) / E  2 
F( IX)*1./E2 

GO  TO  216 

214  I F ( IX-L2)2 12 .212.215 

215  EPSC IX)=SIG< IX)/E3 


F( IX)=1.0/E3 
GO  TO  216 

212  FPS( IX.=STR( 1X1) 

FI  I X ) =S { I XI ) 

216  DO  202  I J=2 • IXi 
I JL= I J-l 
I JR= I J+l 

I F ( I J-L1)203»2C4»205 

204  SG I =$TS ( I J ) 

Vl=V( I J) 

SGV=STS( I JL)-STS( I J) 

VLV=V( I JL ) -V ( I J ) 

CALL  INTERP 

VEL ( I J )  =  ( STS ( I  JR  J-SGL+V ( I JS )/D2+RCF*VL ) / ( RCF+DI ? 

S I G ( I J)=STS( IJR)-(VEL( I J)~V(I JR) >/D2 
EP5( I J)=STR( IJ)+(SIG( IJ)-STS( IJ) )/E2 
F ( I J ) =XE2 

EPSMl  =  EP8tfl+(SIG(IJ)-STS( I J  ) ) *FPM1 
STRS2=SIG( IJ)*SIG( I J) 

FPM1=XE1+XKN*STRS2**XN12 
GO  TO  202 

205  I F ( IJ-L2)206*207*208 
207  OGI =STS ( I J ) 

V I =V ( I J) 

SGV=STS( I JR)-STS( I J) 

VLV=V( I  JR) -VC  I J) 

CALL  INTERP 

VEL ( I J ) = ( 3GL-STS ( I JL ) +RCF*VL+V ( I JL ) /D2 ) / ( RCF+DI ) 
5IG(IJ)=5TS( IJL)  +  (VEL( IJ)-V(IJL)  )/D2 
EPS! 1 J)=STR( I J)+{$IG( I J)-STS( I J) )/S2 
F ( I J ) =XE2 

FPB“3=EPe?/.3+iSIGI  I J )  —STS  l  I  J)  )*FPM3 
$TRS2  =  SIG(  U)*SIG(  IJ) 

FPV3  =  XE1  +  XK.\>STR£2^*XN12 
GO  TO  2C2 

206  DX=D2 
EX  =  'r2 

210  VEL  (  I J  )  =  >'  V  (  I  JR  )  +V  (  I  JL  )  +DX#  (  ST5  ( I  JR  )  -ST  S  ( I  JL  )  )  J/2.0 
SIG(  I J  )  =STS  (  IJL)  +  (VELC  IJ)-V(UL)  )/DX 
EPS( IJ)=STR( IJ)  +  (SIG(IJ)-STS( IJ)  )/EX 
F ( I J ) =XE2 
GO  TO  202 
206  CONTINUE 
203  ?GI=STS(IJ) 

VI=V( I J) 

SGV  =  ST5(  I  JR  ) -STS  t  U) 

VLV=V( I  JR ) —V (  I  J) 

CALL  INTERP 

sgr=sgl 


CR  =  CF 
RCR=RCF 

SGV=STS( I JL)-STS( I J) 

VLV=V( I JL)-V( I J) 

CALL  INTERP 
CL  =  CF 
RCL=RCF 

VEL { I J )  =  ( (SGR-SGL)/RH01+CR*VR+CL*VL)/(CL+CR) 

S I G ( I J ) =SGL+RCL* ( VEL (I J ) - VL ) 

EPS( IJ)=STR( IJ )+S{ I J)*(STG( IJ)-STS( I J) ) 
STRS2=SIGUJ)*SIG(  I  J) 

F (  I  J)=XE1+XKN*STRS2**X\'12 
304  FORMAT (I5»3E12.4) 

2C2  CONTINUE 

DFSTR=SIG(L1)-SIG(L2) 

C  CONVERSION  factor  from  watt-sec  to  joules  is  *con* 

CON=1.0 

C  THE  UNITS  OF  Cl  ARE  AMPERES. 

Cl ( IX)=DFSTR*AREA*C2*DIEC/X2 

C  ENERGY  IS  the  POWER  PRODUCED  BY  THF  SHORT  CIRCUIT  CURRENT 
IF ( IX-MIX 1600,600,700 
6C0  CONTINUE 

ENERGY! IX)=DT/3.*R*CON*(CI ( IX-2) **2+4. *C I i IX~1)**2+ 

1CI ( IX)**2)+ENERGY(IX-2) 

211  PRINT  303,  T , EPS (Ml ) ,SIG(L1 ) ,  EPBM1 , EPS (LI) ,SIG(L2), 
1EPBM3,EPS(L2),DFSTR,CI ( IX ) , ENERGY (IX) 

GO  TO  20C 
700  CONTINUE 

PRINT  303,  T,EPS(M1 ) ,SIG(L1) ♦  EPRM1 *EPS( LI ) »SIG( L2 ) , 
1EPBM3 ,EPS ( L2 ) , DFSTR  ,CI (IX) 

MI/=MIX+2 
200  CONTINUE 

299  FORMAT ( 1H1,1CX,44H  PLASTIC  WAVE  PROPAGATION  IN  A  COM 
IPOS  I TE  BAR/) 

300  FORMAT ( 3X ,3H  E1,7X»3H  E2,7X,3H  E3,6a,5H  RH01,5X, 

I5H  RH02,5X,5H  RH03,5X,3H  L1,2X,3H  L2) 

301  FORMAT (6E10.3»2I5) 

302  FORMAT ( 4X ,5H  TIME,7X,6H  EPSM1,6X,6H  SIGLl,5X, 

17H  FPSBMl ,6X ,6H  EPSL1,6X,6H  SIGL2,5X»7H  EPSBM3,6X, 

26H  EPSL2,6X,6H  DFSTR, 2X,8H  CURRENT ,4X ,7H  ENERGY) 

303  FORMAT ( 11E12.4 ) 

311  FORMAT ( 14H****LEGEND****»/»4X,35HEPSM1=STRAIN  AT  FI 
1RST  GAGE  LOCATION ,/ ,4X ,40HS I GL1- STRESS  IN  MEDI 

2UM  1  AT  CRYSTAL  FACE,/ ,4X ,41HEPSBM1=STRAIN  IN  MED 
3IUM  1  AT  CRYSTAL  FACE, /,4X,40HEPSL1=ST RAIN  IN  MED 
4IUM2  AT  CRYSTAL  FACE, / ,4X ,33HSIGL2=STRESS  AT  CRY 
5STAL  REAR  FACE , / ,4X ,4lHEPSBM3=ST RAI N  IN  MEDIUM  3  A 
6T  CRYSTAL  FACE , / ,4X ,36HEPSL2=STR A I N  IN  CRYSTAL  A 
7T  REAR  FACE) 

312  FORMAT ( 4X »38HDFSTR= STRESS  DIFFERENCE  ACROSS  CRYST 
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1AL,/*4X*40HCURRENT=CIRCUIT  current  in  units  of  am 
2PS./*4X.40HENERGY=CRYSTAL  ENERGY  INUNITS  OF  JOULES*/) 
STOP 
END 

SUBROUTINE  interp 

COMMON  SGI.SGV.VI .VLV.RHOl ,C1 »XK ,XN ,DT ,DX 1 »SGL »VL » 
1XE1.CF.RCF.XXN.XN1 
COMMON  YIELD 

C  YIELD  STRESS  FOR  MATERIAL  1  IS  EQUAL  TO  30  KSI (YIELD) 
IFISGI+SGV-  YIELD  )105*105»106 

105  SGL=SGI+SGV 
VL=VI+VLV 
CF=C1 

GO  TO  107 

106  XN12=XNl/2.0 
DDX=DXl/60.0 
SGLV=SGV/60.0 
SGL=SGI 
XMUL=0.0 

100  sgl=sgl+sglv 

XMUL=XMUL+DDX 
STRS2=SGL*SGL 
104  FP=XE1+XKN*STRS2**XN12 
CF=SQRT(1.C/(RHC1*FP) ) 

IF(XMUL/CF-DT) 102,101.101 
102  GO  TO  100 

101  VL=VI+XMUL*VLV/DX1 

107  RCF=RH01*CF 
RETURN 

END 
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APPENDIX  d 


LISTING  OF  COMPUTER  PROGRAM 

FOR 


CURVE  FIT 


rionnnnnn^nn 


T8IS  PROGRAM  COMPUTES  THE  COEFFICIENTS  FOR  A  NATURAL 
CUBIC  INTERPOLATORY  SPLINE  FIT  TO  EQUALLY  SPACED 
DATA  POINTS. 

INPUT  AND  OUTPUT  VARIABLES  ARE  AS  DEFINED  BELOW 
N  -  NUMBER  OF  SEGMENTS  OR  SPLINES 
H  -  WIDTH  OF  EACH  SEGMENT  { NORN AL I ZE  CLOSE  TO  1 
TO  MINIMIZE  COMPUTATIONAL  ERRORS) 

Y  -  DEPENDENT  VARIABLE  INPUT  DATA 

FACT ( I )  -  THESE  ARE  THE  MULTIPLIERS  IN  THE  GENERAL 
SPLINE  FORMULA.  INDEXED  FROM  1TO  N+l 
DIMENSION  EM(  11 )  *L'(  11 )  *F(  11  )  *3(11)  .FACT ( 13  )  »Y{  13  ) 

D1TA  Y/O. 0,1. 8 .5.3 .10. 4. 15. 9. 20. A, 23. 6. 25.8, 27. 4. 
128.5,29.3,29.9.30.1/ 

H=  1 .0 
N=  12 
NP1=N+1 
NM1=N-1 
DO  10  1=1, NM1 
EM ( I ) =H/6.0 
U( I )=2.0*H/3.0 
FI  I )=H/6.0 
10  CONTINUE 
DO  20  J=2,N 
JM1 = J— 1 
JP1=J+1 

B(  JM1 )  =  (Y( JP1)-2.0*Y{ J)+Y( JM1 )  )/H 
20  CONTINUE 

DO  100  1=2, NM1 
IM1=I-1 

EM ( I )=EM( I )/U( I  Ml ) 

U(  I  )=U( I ) -EM { I ) *F { I  Ml ) 

100  CONTINUE 

FACT  1 1 ) =0 .0 
FACT ( NP1 ) =0.0 
FACT ( 2 ) =3  I  1 ) 

DO  200  1=2, NM1 
I Pl= I +1 

FACT ( I  PI )=3( I )-FM{ I  )*FACT{  I  ) 

200  CONTINUE 

FACT (N )=F ACT (N )/U(NMl) 

DO  30C  J=2»NM1 
I =NP 1- J 
IP1=I+1 
I M 1  =  I  -  1 

FACT ( I )  =  ( FACT ( I )-F( I  Ml ) *FACT ( I  Pi ) )/U( I  Ml) 

300  CONTINUE 

PRINT  1 » ( FACT ( I ) ,  1  =  1, NP1) 

1  FORMAT ( 12F10.3 ) 

CALL  EXIT 

end 

/DATA 

/END 
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APPENDIX  E 


LISTING  OF  COMPUTER  PROGRAM  FOR 


DETERMINATION  OF  ENERGY  OUTPUT  FROM  CURRENT 


F 


i 


> 


fe 


General : 

The  program  "ELCENR"  computes  the  energy  produced  from 
the  experimental  voltage  input.  This  program  instructs  the 
computer  to  perform  the  following  operations  in  the  order 
listed: 

(1)  Read  in  the  circuitry  data 

(2)  Read  in  the  time  increment 

(3)  Read  in  the  scope  deflection  data 

(4)  Calculate  voltage,  current  for  each  time  increment 

(5)  Calculate  the  energy  by  Simpson's  rule 

(6)  Print  out  current  and  energy 
Data: 

DT,  IMAX,  SGAIN ,  R,  TG  (FORMAT  100) 

RSC  (To  be  given  in  the  body  of  the  program) 

List  of  symbols: 

DT  =  Time  increment 
IMAX  =  Total  number  of  time  increments 
SGAIN  =  Scope  vertical  gain 

R  =  Resistance  of  the  strain  gage 
TG  =  Trigger  delay  time 
RSC  =  Short  circuit  resistance 
DEF  =  Scope  deflection  data 
V  =  Voltage 

I  =  Short  circuit  current 
ENERGY  =  Energy  developed 
T  =  Time 


1 
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FLOW  CHART 
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Figure  75.  Computer  Program  Flow  Chart  for  Determination  of  Energy 
Output  from  Current 


JFH  »60  .  3C  »  5C00C.  VAL  ATHUR  *  5791  .UNM, 

run ( s ) 

ELCEN\. 

• 

PROGRAM  ELCENR ( INPUT .OUTPUT ) 

C  **  *4*  ********  ******m**  *************■*..*************** 

C  **RcD'JCTION  OF  CRYSTAL  DATA  FOR  ELECTRICAL  ENERGY** 

Q  ***************************  ************************ 

DIMENSION  I ( 501 )»V(5C1) .ENERGY (bci) .DEF ( 501 ) 

REAL  I 
r 

C  READ  IN  TIME  STEP  DT,  tOTAL  NUMBER  OF  TIME  I NCR 5 MEN TS 
C  I  MAX.  SCOPE  VERTICAL  GAIN  SGAIN.  GAGE  RESISTANCE  R. 

C  and  trigger  DELAY  TG 

C  DATA  MUST  BE  GIVEN  FOR  EACH  SET  ACCORDING  TO  FORMATS 
C  100  AND  101 

C  DIMENSIONS.  SGAIN=vOLTS/CFNTIMETER.  R=OHMS .DT=SECONDS ♦ 

C  TG=SECONDS 

READ  100. DT.IMAX. SGAIN. R.TG 
100  FORMAT (E10.3.I10.3E1C.3) 

C 

C  RSC  =  SHORT  CIRCUIT  RESISTANCE 
RS^= 1.33 
C 

IMAXP1=(2*IMAX)+1 

c 

C  READ  IN  SCOPE  DEFLECTIONS  AT  THE  PRESCRIBED  TIME 

C  intervals 

C  PUNCH  SCOPE  DEFLECTION  DATA  AT  PRESCRIBED  TIME  INTERVAL 
C  AS  DETERMINED  BY  THE  INTEGRATION  METHOD 
C  LET  EACH  VALUE  CORRESPOND  TO  1/2  TIME  STEP  (1/2*DT) 

C 

C  INITIALIZE  VALUES  TO  ZERO 
DC  33  J=1 » IMAXP1 
DEF( J)=0.0 
V( J) =0.0 
I  (  J) =0.0 

33  ENERGY! J) =0.0 
C 

READ  10 1 » ( DEF ( J  )  ,  J= 1  , 1 MAXP 1  ) 

101  FORMAT! 16F5. 2) 

DO  1  M=I  ♦  IMAXP1 
C 

C  CALCULATE  VOLTAGE 

'/(M  )=SGAIN*DEF(M) 

C 

C  CALCULATE  CURRENT 
I (M) =V(M) /RSC 
C 

1  CONTINUE 
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n  n  n  n  on 


PRINT  150 
PRINT  151 
450  PRINT  169 

169  FORMAT ( 15X  *  30HSH0T  NO  1206  8/7/66  SCOPE  NO  3  * / ) 

180  PRINT  152*1 MAX  »SGAIN»R»TG 
PRINT  153 
PRINT  154 

CONVERT  UNITS  OF  ENERGY' N)  FROM  WATT-SEC  TO  JOULES 
C0N=1 . 0 

T  =  0.0 

ZERO  TIME  SPECIFIED  2Y  SUBSCRIPT  (1). 

ENERGY ( 1 ) -DT /6.*R* ( I ( 1 ) **2+4.*I ! 2 ) **2+I (3)** 2 )*CON 
PRINT  200»T*V(1)*I(11*ENERGY(1) 

C 

C  CALCULATE  ENERGY  from  CURRENT  USING  SIMPSONS S  RULE  FOR 

c  numerical  integration 

C  ENERGY  =  1/ 3* ( B-A ) /2.\'*R* ! I ( N-l ) **2+4* I ( N ) **2+I !N+1)**2) 

C  INCREMENTS  IN  THAT  INTERVAL.  N  MUST  BE  INTEGER. 

C  example 

DO  2  N  =  3*  IMAXP1 *2 
T=T+DT 

ENERGY! N)=DT/6.*RM I  I ( N-2 I **2+4.* I ( N-l ) **2+It N ) **2 ) * 
lCON+ENERGY(N-2) 

PRINT  20C»T*V(N)*I (N)*ENERGY(N) 

2  CONTINUE 

150  FORMAT! 1H 1,1 CX,34HWAVE  PROPAGATION  IN  COMPOSITE  BAR,/) 

151  FORMAT! 15X,31H  EXPERIMENTAL  ELECTRICAL  ENERGY,//) 

152  FORMAT  (  5X  ,6H  I  MAX ••=  ,  I  5  »/ ,  5X  ,7H  SGAIN=  ,  E10. 3  ♦  /  ,  5X , 

13H  R=,F8.3,5X,7H  DELAY®, E10.3, // ) 

153  FORMAT ( 5X , 5H  TTME,6X,8H  VOLTAGE, 5X , 8H  CURRENT, 5X, 

17H  -NERGY) 

154  FORMAT ( 5X ,6H  (SEC)»5X,8H  ( VOLTS) ,6X , 7H  (AMPS)»4X, 

19H  (JOULES),/) 

200  FORMAT (2Y , El  0. 3 , 3X ,F10.3 * 3X ,F10. 3»2X »E10. 3 ) 

1000  CONTINUE 
C 

C  IMPUT  DATA  ACCORDING  TO  FORMATS  100  AND  101 
STOP 
END 
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EXAMPLE  RUN 


RSC  =  1.33 


53.47 


DT  =  0.183Msec 


TG  =  30pisec 


SGAIN  =5.00 


IMAX  =140 


TIME  CURRENT  CAMPS)  ENERGY ( JOULES ) 


0. 

.0001 

9c22  3E-14 

3. 650E-07 

.0001 

3.505E-13 

7.3Q0E-G7 

•  c  002 

9. 7 63 E-. 13 

1. 095E-C6 

.C  004 

2.681E-12 

1.46GE-C6 

.  0  OG  6 

6. 2P*+E“12 

1.  825E-Q6 

.0  °C 6 

1.319E-11 

2. 190E-06 

.rni 

2.667E-11 

2.  555E-06 

•  i  014 

5.548E-11 

2* 92OE-06 

.U18 

1.0  2  E-10 

3.  285E-C6 

.  u  0  2  3 

1.822E-10 

3. 650E-G6 

.0028 

3.0  HE- 10 

4. C15E-06 

•  0  035 

4.3C6E-10 

4.380E-06 

.0  041 

7 • 358E-1D 

4. 745E-06 

.CO  47 

1. j91E-C9 

5.110E-06 

.0  0  54 

1.556E-09 

5.475E-U6 

.0059 

2.155E-G9 

5. 64GE-06 

•  C  063 

2.833E-09 

6. 205 E- 06 

.0  0  63 

3 • 52 7 E- 09 

6. 57CE-06 

.0069 

4.485E-C9 
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